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Chapter  1 
Introduction 


A  time  domain  version  of  the  uniform  geometrical  theory  of  diffraction  (UTD),  hence¬ 
forth  abbreviated  as  TD-UTD,  is  developed  in  this  study  for  describing  in  relatively 
simple  form  the  transient  radiation  and  scattering  from  a  class  of  perfectly  conduct¬ 
ing  configurations  when  they  are  excited  by  an  electromagnetic  pulse.  In  particular, 
a  TD-UTD  formulation  is  obtained  for  describing  the  scattering  from  an  arbitrary 
curved  wedge,  and  also  the  scattering  from  an  arbitrary  convex  surface.  Due  to  the 
fact  that  the  TD-UTD  is  developed  via  an  analytical  inversion  into  the  time  domain 
(TD)  of  the  corresponding  high  frequency  asymptotic  results  obtained  previously  in 
the  frequency  domain  based  UTD,  the  TD-UTD  therefore  provides  an  “early  time” 
asymptotic  representation  that  is  valid  close  to  the  arrival  times  of  each  of  the  various 
wave  fronts.  More  specifically,  each  term  in  the  TD-UTD  formulation  is  asymptot¬ 
ically  valid  at  and  near  a  particular  arrival  time  associated  with  the  path  length  of 
the  corresponding  ray  (i.e.  incident  ray,  reflected  ray,  diffracted  ray,  etc.)  which  con¬ 
tributes  at  a  given  observation  point.  Due  to  this  “early  time”  asymptotic  nature  of 
the  TD-UTD  representation,  each  TD  ray  field  may  not  have  a  causal  impulse  re¬ 
sponse  since  there  may  be  information  in  the  time  domain  occurring  after  and  before 
the  arrival  time  of  the  ray.  This  is  an  important  observation  since  earlier  attempts 
at  the  interpretation  of  approximate  impulse  responses  of  scattering  targets  did  not 
fully  understand  the  possible  non-causal  nature  of  “early  time”  solutions  [1].  Since 
the  TD-UTD  employs  the  same  rays  as  the  corresponding  frequency  domain  UTD, 
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it  provides  a  clear  physical  picture  for  radiation  and  scattering  in  exactly  the  same 
manner  as  does  the  UTD  on  which  it  is  based. 

The  development  of  the  TD-UTD  is  based  on  an  integral  transform  which  is  called 
the  analytic  time  transform  (ATT)  in  this  report.  The  ATT  is  a  one-sided  transform 
which  transforms  a  frequency  domain  function  (which  is  a  function  of  a  real  frequency 
variable  o>)  into  an  analytic  time  function  which  is  analytic  in  the  upper  half  time 
plane  (Im  (f)  >  0).  By  construction,  the  real  part  of  the  "analytic  time  function 
when  evaluated  for  real  time  (  Im  (£)  =  0)  is  exactly  the  inverse  Fourier  transform 
of  the  frequency  domain  function  whenever  the  frequency  domain  function  obeys  the 
conjugate  symmetry  relationship  (F(—lo)  =  F*(u>))  which  guarantees  that  the  inverse 
Fourier  transform  is  a  real  time  function.  The  imaginary  part  of  the  analytic  time 
function  is  the  Hilbert  transform  of  the  real  part  as  expected  from  complex  function 
theory  since  the  analytic  time  function  is  analytic  in  the  upper  half  time  plane. 

The  ATT  is  developed  because  the  traditional  Fourier  transform  and  Laplace 
transform  are  not  convenient  for  transforming  the  asymptotic  high  frequency  UTD 
results  into  the  time  domain.  The  main  reason  for  the  difficulty  in  using  the  traditional 
Fourier  or  Laplace  transforms  on  the  frequency  domain  UTD  results  lies  in  the  fact 
that  the  resulting  time  functions  may  not  be  one-sided  (i.e.  may  not  be  causal).  Such  a 
difficulty  can  occur  for  even  the  simplest  of  examples,  such  as  an  incident  geometrical 
optics  ray  which  is  reflected  from  a  smooth  concave  surface  and  passes  through  one 
caustic  before  reaching  the  observer.  The  reflected  ray  field  after  the  caustic  has  an 
extra  factor  of  eJ7T^2  in  the  frequency  domain  which  makes  the  use  of  the  conventional 
inverse  Laplace  or  Fourier  transforms  inconvenient.  In  principle,  a  two-sided  Laplace 
transform  or  a  Fourier  transform  can  be  used  to  derive  the  results  in  this  report,  but 
would  require  a  much  greater  effort  and  the  resulting  time  domain  formulas  would 
have  been  much  more  complicated. 

In  addition  to  the  convenience  of  using  the  ATT  on  the  frequency  domain  UTD 
results,  there  is  the  added  benefit  of  the  ability  of  performing  efficient  convolutions 
with  a  broad  class  of  excitation  pulses.  An  algorithm  is  developed  to  convolve  any 
analytic  time  function  obtained  by  using  the  ATT  with  an  excitation  with  a  frequency 
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response  which  can  be  written  as  a  series  of  complex  exponentials.  This  is  a  very 
general  description  of  an  excitation  pulse,  but  it  is  especially  useful  when  the  pulse 
has  a  narrow  temporal  width  and  a  broad  banded  frequency  response.  When  this  is 
the  case,  the  excitation  pulse  can  be  modeled  with  only  a  few  basis  functions  (around 
10  or  15  basis  functions  for  many  practical  cases)  and  the  convolution  is  a  simple 
summation  over  a  few  terms. 

Although  the  term  “UTD”  encompasses  a  number  of  different  asymptotic  solutions 
for  various  special  geometries,  only  two  specific  but  most  often  required  geometries 
are  analyzed  here;  namely,  an  arbitrary  curved  wedge,  and  an  arbitrary  smooth  con¬ 
vex  surface.  These  geometries  can  be  used  to  build  up  the  solution  to  a  more  general 
complex  configuration  (in  terms  of  wedges,  plates,  and  cylinders  or  ellipsoids,  etc.) 
when  the  standard  assumption  of  localization  of  UTD  fields  is  invoked.  A  TD-UTD 
formulation  for  the  transient  scattering  from  a  pulse  excited  general  curved  wedge  is 
accomplished  by  applying  the  ATT  to  the  well  understood  frequency  domain  UTD 
formulas  that  have  been  obtained  previously  by  Kouyoumjian  and  Pathak  [2].  In 
particular,  the  transient  response  of  the  various  ray  mechanisms — such  as  incident, 
reflected  and  diffracted  rays — is  explicitly  obtained  in  closed  form  when  the  excita¬ 
tion  is  a  time  impulsive,  general  astigmatic  wavefront.  There  are  explicit  and  easy 
to  calculate  formulas  for  the  TD-GO  rays,  both  incident  and  reflected,  for  the  first 
order  diffracted  ray  and  also  for  the  slope  diffracted  ray.  The  slope  diffraction  con¬ 
tribution  is  important  when  the  incident  field  at  the  edge  is  rapidly  varying  with 
respect  to  space  (not  time).  Likewise,  the  TD-UTD  solution  for  the  diffraction  by  a 
smooth  convex  surface  is  developed  via  the  ATT,  by  analytical  inversion  of  the  corre¬ 
sponding  frequency  domain  UTD  solution  obtained  previously  by  Pathak  fit  al.  [3,  4]. 
The  analytic  time  functions  obtained  from  this  transformation  can  not  be  obtained 
completely  in  closed  form,  but  efficient  algorithms  are  developed  for  their  calculation 
where  new  special  functions  are  defined  for  the  computation  of  the  TD-UTD  con¬ 
vex  surface  diffraction.  The  scattering  case  where  the  observer  and  the  source  are 
removed  from  the  surface  is  analyzed  here;  on  the  other  hand,  the  case  of  a  source 
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(pulsed  antenna)  and/or  the  observer  located  on  the  surface  can  be  dealt  with  in  a 
similar  manner  in  the  future. 

The  interest  in  transient  electromagnetic  analysis  where  the  excitation  is  a  rel¬ 
atively  short  pulse  has  been  increasing  over  the  years  and  it  constitutes  an  area  of 
research  which  will  continue  to  receive  attention.  This  is  mainly  due  to  the  increased 
use  of  ultra-wide  band  radars  and  the  development  of  impulse  radiating  systems. 
Also,  there  has  been  a  long  time  interest  in  the  effects  of  natural  and  manmade  short 
electromagnetic  pulses  (EMPs)  on  communication  and  radar  systems.  It  is  most  nat¬ 
ural  to  analyze  such  transient  problems  directly  in  the  time  domain.  This  is  especially 
true  when  the  scattering  objects  are  large  in  terms  of  pulse  width. 

To  analyze  the  scattering  or  radiation  from  complex  structures,  one  could  alterna¬ 
tively  perform  a  frequency  domain  UTD  analysis  at  a  series  of  frequencies,  then  simply 
apply  an  appropriate  window  and  use  the  inverse  fast  Fourier  transform  (IFFT)  to 
obtain  the  scattered  field  due  to  a  transient  pulse  excitation.  However,  a  direct  TD- 
UTD  analysis  of  the  same  problem  is  more  efficient  when  the  pulse  is  very  narrow  (for 
example  a  pulse  with  a  frequency  response  with  more  than  100  percent  bandwidth). 
Also,  the  TD-UTD  can  be  used  to  compute  an  approximate  impulse  response  which 
can  usually  provide  more  insight  since  the  impulse  response  is  essentially  independent 
of  the  excitation,  whereas  the  response  to  a  particular  excitation  is  obtained  from  the 
convolution  operation.  The  TD-UTD  provides  a  new  way  of  examining  electromag¬ 
netic  radiating  and  scattering  phenomena  as  compared  to  what  has  been  done  in  the 
past,  and  it  therefore  could  lead  to  new  insights  for  other  applications.  Some  other 
applications  which  may  result  from  this  TD-UTD  development  are  parametric  mod¬ 
els  used  for  radar  data  analysis  [5,  6]  and  also  the  development  of  a  new  hybrid  time 
domain  integral  equation  based  numerical  method. 

There  exist  various  methods  for  the  numerical  analysis  of  transient  electromag¬ 
netic  problems  such  as  the  finite  difference  time  domain  (FDTD)  method  or  the  time 
domain  integral  equation  (TD-IE)  approaches.  These  methods  have  been  found  to 
be  very  useful  for  a  great  many  problems,  but  become  very  cumbersome  or  even  in¬ 
tractable  for  problems  where  the  geometry  is  large  in  terms  of  pulse  width  (or  in  terms 
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of  the  smallest  wavelength  of  interest).  Because  of  this,  the  TD-UTD  should  provide 
a  very  useful  complement  to  these  more  robust  numerical  methods  when  analyzing 
general  problems.  Also,  the  robust  numerical  methods  such  as  the  FDTD  and  the 
TD-IE  are  not  well  suited  for  computing  approximate  impulse  responses. 

There  is  an  abundance  of  research  which  has  been  done  in  the  past  on  transient 
electromagnetic  (or  acoustic)  analysis.  In  fact,  there  are  some  books  dedicated  to 
this  subject  [7,  8,  9].  There  is  also  a  large  amount  of  literature  on  the  subject  of 
numerical  time  domain  methods  which  are  not  discussed  here.  There  is  some  work  on 
the  formulation  of  a  time  domain  version  of  the  UTD  for  wedge  diffraction  [10, 11, 12], 
but  this  previous  work  concentrates  on  the  time  domain  analysis  of  the  diffraction 
from  a  straight  perfectly  conducting  wedge.  Unfortunately,  this  analysis  uses  an 
inverse  Laplace  transform  to  obtain  the  time  domain  formulas  and  is  thus  limited  to 
simple  excitations  and  straight  wedges.  In  contrast,  the  work  in  this  report  can  easily 
handle  general  astigmatic  incident  ray  fields  which  may  have  passed  through  a  caustic 
before  arriving  at  the  edge  of  the  wedge.  Also,  in  this  report,  the  wedge  is  a  general 
curved  wedge  which  can  have  curved  faces  or  a  curved  edge.  Other  research  which 
is  most  similar  to  the  TD-UTD  development  in  this  report  is  the  spectral  theory  of 
transients  (STT)  [13,  14,  15,  16,  17,  18,  19].  In  the  STT,  an  analytic  time  transform 
is  used  which  is  essentially  the  same  as  the  ATT  of  this  report  except  that  their  time 
convention  for  the  frequency  domain  fields  is  e~,ut  whereas  the  time  convention  used  in 
this  report  is  eJU,(.  The  important  difference  between  the  STT  and  the  TD-UTD  is  that 
the  STT  is  used  to  find  the  exact  time  domain  solution  to  canonical  geometries  which 
can  then  be  approximated  for  “early  time”,  whereas  the  TD-UTD  obtains  an  “early 
time”  solution  directly  by  using  the  ATT  on  the  frequency  domain  UTD  formulas. 
In  principle,  the  STT  is  a  very  powerful  and  accurate  method  but  unfortunately 
the  “early  time”  solutions  obtained  by  the  STT  method  are  not  generalized  to  handle 
relatively  arbitrary  geometries  such  as  the  arbitrary  curved  wedge  analyzed  in  Chapter 
3  of  this  report.  Furthermore,  it  appears  to  be  rather  complicated  to  extend  the  STT 
to  deal  with  smooth  convex  surfaces  as  are  analyzed  in  Chapter  4  of  this  report.  The 
STT  has  been  used  to  obtain  some  interesting  solutions  for  straight  wedge  diffraction 
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due  to  a  collimated  pulsed  beam  type  illumination.  It  should  be  mentioned  that  the 
STT  also  provides  inspiration  especially  for  the  use  of  the  analytic  time  transform  in 
the  present  research  which  may  not  have  been  accomplished  otherwise. 

This  report  is  organized  as  follows.  Chapter  2  discusses  the  analytic  time  trans¬ 
form  (ATT).  The  general  properties  of  the  ATT  and  its  relations  to  the  traditional 
Fourier  and  Laplace  transforms  are  given.  The  ATT  of  a  frequency  domain  function 
produces  a  complex  time  function  which  is  analytic  in  the  upper  half  t-plane,  where 
t  is  the  complex  time  variable.  But  there  are  special  properties  of  the  analytic  time 
function  when  it  is  evaluated  on  its  boundary  of  analyticity,  the  real  time  axis.  These 
special  properties  are  explained  in  Chapter  2.  There  are  times  when  the  ATT  of  a 
particular  frequency  domain  function  can  not  be  obtained  in  closed  form  and  therefore 
it  is  important  to  be  able  to  calculate  the  analytic  time  function  for  complex  time 
values.  For  some  functions,  an  early  time  representation  can  be  obtained  from  an 
asymptotic  high  frequency  domain  expansion;  while,  a  late  time  representation  can 
be  obtained  from  a  frequency  domain  power  series  (i.e.  low  frequency)  expansion. 
When  these  two  representations  overlap  for  intermediate  time  values,  then  these  rep¬ 
resentations  can  be  used  to  calculate  the  ATT  of  the  original  function.  This  procedure 
is  explained  in  Chapter  2  and  a  particular  example  is  provided  for  demonstration.  Fi¬ 
nally,  a  method  for  performing  an  efficient  convolution  on  a  broad  class  of  excitation 
functions  with  the  analytic  time  impulse  response  is  presented.  Further,  Appendix 
A  presents  a  simple  method  for  synthesizing  a  given  smooth  transient  pulse  function 
in  terms  of  a  set  of  special  basis  functions  which  provide  the  efficient  closed  form 
convolution. 

Chapter  3  presents  the  derivation  of  the  TD-UTD  analysis  of  the  scattering  from 
a  perfectly  conducting  curved  wedge  where  the  curved  wedge  may  have  curved  faces 
and  also  a  curved  edge.  There  are  four  important  scattering  mechanisms  for  this 
geometry;  namely,  the  geometrical  optics  (GO)  ray  fields  (an  incident  ray  field  and 
a  reflected  ray  field),  the  first  order  edge  diffracted  field,  and  the  higher  order  slope 
‘  diffracted  field.  The  incident  GO  ray  may  have  any  polarization  including  circular 
polarization  or  in  general  elliptical  polarization  and  the  ray  tube  may  be  an  astigmatic 
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ray  tube  which  has  possibly  passed  through  some  caustics.  The  impulse  response  in 
Chapter  3  is  defined  in  terms  of  the  temporal  behavior  of  the  incident  field  (which  is  an 
analytic  delta  function  for  the  impulse  response  computation);  while,  the  geometrical 
and  polarization  properties  of  the  incident  field  are  completely  general.  Two  versions 
of  the  TD-UTD  slope  diffraction  coefficient  are  derived;  one  is  based  on  the  frequency 
domain  slope  diffraction  coefficient  by  Veruttipong  and  Kouyoumjian  [20];  while,  the 
other  is  based  on  the  frequency  domain  slope  diffraction  coefficient  by  Hwang  and 
Kouyoumjian  [21].  The  details  of  the  derivation  of  these  slope  diffraction  coefficients 
are  presented  in  Appendix  C.  These  two  slope  diffraction  formulations  are  compared 
in  the  numerical  examples  in  Chapter  3.  An  exact  solution  for  the  scattering  from 
a  straight  wedge  when  it  is  illuminated  by  an  electric  or  magnetic  dipole  is  used  as 
a  reference  here.  This  exact  solution  is  based  on  Felsen’s  work  [22]  and  is  summa¬ 
rized  in  Appendix  B.  The  numerical  examples  also  reveal  how  the  impulse  response 
approximated  by  the  asymptotic  ’’early  time”  TD-UTD  may  not  be  causal. 

Chapter  4  presents  the  derivation  of  the  TD-UTD  analysis  for  the  scattering  from 
a  smooth  convex  surface.  The  scattering  case  is  the  only  one  considered  here  where 
the  source  and  observer  are  removed  from  the  boundary.  Smooth  convex  surface 
diffraction  can  be  modeled  by  creeping  wave  modes  when  the  observer  is  in  the  deep 
shadow  region,  or  for  multiple  encirclements  around  closed  bodies.  These  creep¬ 
ing  wave  modes  are  dispersive  in  the  frequency  domain  and  therefore  their  impulse 
response  is  effectively  delayed  in  time.  In  other  words,  the  impulse  response  of  a 
creeping  wave  mode  starts  with  zero  amplitude  and  gradually  builds  up  to  a  peak 
value  occurring  after  the  arrival  time  corresponding  to  a  simple  surface  diffracted 
ray  path.  The  ATT  of  the  frequency  domain  fields  for  smooth  surface  diffraction  is 
not  obtained  in  closed  form,  so  a  numerical  procedure  is  developed  to  compute  the 
ATT  of  the  creeping  wave  mode  with  extreme  efficiency.  This  numerical  algorithm 
is  presented  in  Appendix  D.  Chapter  4  also  includes  the  development  of  a  uniform 
TD-UTD  formulation  which  is  valid  in  the  lit  region,  the  shadow  region,  and  across 
the  shadow  boundary.  To  obtain  this  uniform  result,  the  ATT  of  the  Fock  type 
functions  occurring  in  the  corresponding  frequency  domain  UTD 
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representation  are  needed.  These  transforms  are  not  obtained  in  closed  form,  but 
instead  approximate  expansions  for  early,  late  and  intermediate  time  values  are  ob¬ 
tained,  and  the  ATT  of  the  Fock  type  functions  are  thus  calculated  using  an  algorithm 
which  chooses  the  appropriate  representation  depending  on  the  time  variable.  This 
numerical  algorithm  is  presented  in  Appendix  E 

Finally,  some  conclusions  are  made  in  Chapter  5.  This  chapter  also  discusses 
some  ideas  for  continuing  the  development  of  the  TD-UTD  so  that  even  more  general 
geometries  can  be  analyzed.  The  possible  impact  of  the  present  TD-UTD  development 
along  with  the  use  of  the  ATT  is  also  discussed.  In  particular,  one  of  the  possible 
applications  which  might  benefit  from  the  TD-UTD  is  the  use  of  parametric  models 
in  the  time  domain  for  the  analysis  of  radar  scattering  data. 

An  eJU,<  time  convention  for  the  frequency  domain  fields  is  assumed  and  suppressed 
in  the  following  analysis. 


8 


Chapter  2 

Analytic  Time  Functions 


The  Analytic  Time  Transform  (ATT)  is  fully  explained  in  this  chapter.  Since  many 
texts  define  the  conventional  Fourier  Transform  (FT)  or  the  Laplace  Transform  (LT) 
in  slightly  different  ways,  these  conventional  transforms  are  explicitly  defined  here. 
Also,  the  relationship  between  all  three  transforms  is  presented.  The  properties  of 
the  ATT  are  then  explored  including  the  important  time  convolution  formula.  When 
the  analytic  time  function  obtained  from  the  ATT  is  evaluated  on  the  real  time  axis 
(Im  t  =  0),  some  very  useful  properties  arise  and  these  real  time  properties  are  ex¬ 
plained.  Some  important  examples  of  analytic  time  functions  which  are  useful  in  the 
TD-UTD  are  presented  in  an  easy  to  use  table.  Also,  the  ATT  of  a  generic  asymp¬ 
totic  high  frequency  series  and  a  generic  low  frequency  power  series  are  investigated. 
The  ATT  of  an  asymptotic  high  frequency  series  or  a  low  frequency  power  series  pro¬ 
duces  an  early  time  or  late  time  representation,  respectively.  These  “early”  time  and 
“late”  time  representations  are  sometimes  useful  in  the  calculation  of  special  func¬ 
tions.  A  method  for  efficiently  computing  the  time  convolution  of  two  analytic  time 
functions  is  also  presented.  This  efficient  convolution  is  especially  important  for  the 
implementation  of  the  TD-UTD. 

It  is  appropriate  at  this  juncture  to  describe  some  of  the  notation  used  in  this 
chapter.  Throughout  this  chapter,  time  domain  functions  which  are  real  functions  of 
real  time,  are  written  in  the  usual  way,  for  example 

F(t)  for  real  t. 
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Note  that  these  real  time  functions  (when  used  with  the  Fourier  transform)  may  be 
two  sided  (i.e.  not  causal).  The  topic  of  causality  is  further  discussed  later.  On  the 
other  hand,  the  analytic  time  functions  obtained  from  the  ATT  are  analytic  (complex) 
functions  of  complex  time  and  are  denoted  by  a  ’+’  sign  on  top,  such  as 

F  (t)  for  complex  t. 

The  frequency  (or  Fourier)  domain  and  Laplace  domain  functions  are  denoted  by  a 
on  top.  A  Laplace  domain  function  is  an  analytic  function  of  a  complex  variable 
s,  for  example 

F(s)  for  complex  s  =  cr  +  ju. 

where  a  and  u  are  real.  In  many  instances,  the  frequency  domain  function  F(w)  can 
be  obtained  from  an  evaluation  of  the  Laplace  domain  function  F(s)  at  imaginary 
s  values  (i.e.  a  =  0),  which  motivates  the  use  of  similar  notation  for  the  frequency 
domain  functions,  such  as 

F(u>)  for  real  u>. 

Although  the  notation  for  F(s)  and  F(u>)  appear  the  same,  there  are  many  instances 
when  the  functions  are  different.  For  example,  if  F(s )  =  1/s  then  F(uj)  =  + 

pv  [l/(yu>)]  where  “pv”  denotes  Cauchy’s  principle  value.  Nonetheless,  the  same 
notation  is  used  for  convenience,  and  it  is  obvious  from  the  context  whether  the 
function  is  a  Laplace  domain  function  F(s)  or  a  frequency  (Fourier)  domain  function 
F(u>). 

The  word  “function”  is  used  loosely  here  since  F(t),  or  F(cj),  may  in  general 
be  a  distribution  (or  what  is  sometimes  called  a  generalized  function).  Although 
'  most  mathematicians  will  not  be  happy  with  the  use  of  the  word  “function”,  it  is 
nevertheless  used  in  this  loose  sense  throughout  this  report.  This  is  consistent  with 
much  of  the  engineering  literature  which,  for  example,  calls  the  impulse  distribution 
£(f)  as  an  impulse  function. 

All  of  the  frequency  domain  functions,  such  as  F(uj),  are  assumed  to  be  of  expo¬ 
nential  order  a  where  a  <  0.  This  means  that  a  is  the  largest  constant  such  that 
F(w)ea 1“!  — >  0  when  |u;|  — ►  oo.  Notice  that  there  may  be  times  when  a  is  arbitrarily 


10 


less  than  0  but  not  equal  to  zero.  For  example,  when  F(u>)  =  u>  then  a  =  e  <  0  where 
e  is  a  constant  arbitrarily  close  to  0. 

2.1  Definitions 

The  Fourier  transform  as  used  in  this  report  is  defined  as  follows.  The  Fourier  trans¬ 
form  is  an  integral  transform  that  relates  a  time  domain  function  F(t),  which  is  a 
function  of  a  real  time  variable  t,  to  a  frequency  domain  function  F(u;)  which  is  a 
function  of  a  real  frequency  variable  uj.  This  transform  is  defined  as 

r°° 

F(u>)  =  /  F{t)e~3Wtdu  (2.1) 

J—oo 

and  the  inverse  Fourier  transform  is 

1  r°°  - 

m  =  £  j_x  F(w)e“*it  (2.2) 

The  above  convention  for  the  Fourier  transform  is  what  is  typically  used  in  engineering 
texts.  The  Fourier  transform  may  also  be  written  as  an  operator  T  as  in 

F(W)  =  F|F(«)] 

and  the  inverse  Fourier  transform  denoted  by 

m  =  r-'  [FH] 

or  one  can  write  the  relationship  between  the  pair  of  functions  simply  as 

F(t)  F( u). 

It  should  be  kept  in  mind,  for  the  Fourier  transform,  that  the  time  domain  func¬ 
tion  F(t)  and  the  frequency  domain  function  F(uj)  are  actually  distributions.  Also, 
although  in  general  the  Fourier  transform  is  defined  for  complex  time  functions,  in 
this  report  the  time  domain  functions  (without  the  +  on  top)  is  always  assumed  real. 
This  implies  the  following  conjugate  symmetry  relationship  for  the  frequency  domain 
function 

F(-v)  =  F»  (2.3) 
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where  *  denotes  a  complex  conjugate.  This  conjugate  symmetry  relationship  is  used 
when  finding  the  connection  between  the  ATT  and  the  Fourier  transform.  Finally, 
notice  that  F(t)  may  be  a  two-sided  function;  this  fact  turns  out  to  be  an  important 
property  when  one  transforms  asymptotic  high  frequency  solutions  into  the  time 
domain  which  is  examined  later  in  this  chapter. 

The  Laplace  transform  is  defined  as  follows.  In  this  report,  only  the  one-sided 
Laplace  transform  is  discussed.  Although  there  are  some  benefits  from  using  the  two- 
sided  Laplace  transform  (such  as  the  ability  to  transform  two-sided  time  functions) 
only  the  one-sided  version  is  used  here  for  two  reasons.  The  first  reason  is  that  there  is 
not  a  sufficient  benefit  from  using  the  two-sided  Laplace  transform  on  the  applications 
discussed  in  this  report.  On  the  other  hand,  the  ATT  is  far  more  versatile  especially 
when  dealing  with  asymptotic  high  frequency  results.  The  second  reason  is  that  the 
one-sided  version  of  the  Laplace  transform  is  more  widely  used  and  understood.  The 
Laplace  transform  is  defined  by 

F(s )  =  f  F(t)e~stdt  for  Re  s  >  8  (2-4) 

Jo 

and  the  inverse  Laplace  transform  is 

1  f£+j  oo  _ 

F(t )  =  —  /  F(sy'ia  (2.5) 

ZTVJ  Jt-j oo 

where  e  >  6  and  F(t)  is  of  exponential  order  8  (i.e.  8  is  the  largest  value  such  that 
F(t)esM  — »  0  when  t  — >  oo).  The  Laplace  transform  is  defined  whenever  the  integral 
in  Equation  (2.4)  exists,  and  can  be  extended  to  include  distributions  for  which  the 
integral  will  not,  in  general,  exist  in  the  usual  sense  [23].  The  Laplace  transform  may 
also  be  written  as  an  operator  C  as  in  » 

F(s)  =  C  [F(t)l 

and  the  inverse  Laplace  transform  denoted  by 

F(t)  =  [F(S)] 

or  one  can  write  the  relationship  between  the  pair  of  functions  simply  as 

m  ~ 
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Now  the  analytic  time  transform  (ATT)  is  defined  as  follows.  There  are  two 
versions  of  this  transform.  One  version  transforms  a  frequency  domain  function  F( a>) 
into  a  complex  analytic  time  function  F  (<)>  while  the  other  version  transforms  a  real 
time  domain  function  F(t)  into  a  complex  analytic  time  function  F  (t).  When  F(u>) 
is  the  Fourier  transform  of  F(t)  then  both  versions  produce  the  same  F  (<)•  The  first 
version  is  defined  as  follows 

F  (t)  =  —  /  F(u)u(u)eJwtdu)  for  Im  t  >  a  (2-6) 

7T  J-oo 

When  F(tu)  contains  no  impulses  at  u  =  0,  the  above  integral  can  be  written  as 

F  (t)  =  —  f  F(w)e3Utdu)  for  Im  t  >  a  (2-7) 

ir  Jo  ' 

The  corresponding  inverse  ATT  is 

1  r°°+Je  +  •  . 

F(l>)u(u>)  =  -  F  ( t)e  3U)tdt  for  e  >  a  (2-8) 

2  J  —  oo+je 

where  u(ui)  is  the  Heaviside  unit  step  function  defined  as 

r 

0  u>  <  0 

u{ u>)  =  <  1/2  u  =  0  (2.9) 

1  u  >  0 

s. 

This  first  version  of  the  ATT  may  be  also  written  as  an  operator  Aw  as  in 

F  (0  =  Aw  [#(«)] 

and  the  inverse  ATT  is 

F{u)u( u)  =  AZ1  [F  (*)] 

or  one  can  write  the  relationship  between  the  pair  of  functions  simply  as 

F  (0  F(v). 

The  second  version  of  the  ATT  is  defined  by 

F  (t)  =  —  f  — — ^-dr  for  Im  t  >  0  (2.10) 

7 T  J-oo  t  —  T 
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(where  a  <  0  has  been  assumed).  Equation  (2.10)  is  sometimes  called  the  Cauchy 
Representation  [24].  The  inverse  ATT  for  the  second  version  is 

F(t)=  Re  [F(tj\  for  Im  <  =  0  (2.11) 

where  Re  (z)  denotes  the  real  part  of  a  complex  z.  Notice  that  (2.11)  uses  F  (t) 
at  Im  t  =  0  even  though  in  (2.10)  it  is  only  defined  for  Im  t  >  0.  If  F  ( t )  turns 
out  to  be  analytic  at  Im  t  =  0  then  this  poses  no  difficulty,  on  the  other  hand,  if 
F  ( t )  is  not  analytic  on  the  real  time  axis  ( Im  t  =  0)  then  F  (0  at  Im  t  =  0  is  the 
distribution  defined  by  F  ( t )  (Im  t  >  0)  on  its  boundary  of  analyticity  [24,  25].  One 
may  refer  to  the  book  by  Beltrami  and  Wohlers  for  a  thorough  explanation  of  the 
connection  between  distributions  and  analytic  functions  evaluated  on  their  boundaries 
of  analyticity  [25].  The  second  version  of  the  ATT  in  (2.10)  may  also  be  written  as 
an  operator  At  as  in 

F  (t)  =  A,  [F(l)] 

and  the  inverse  ATT  is 

F(t)  =  Ai1  [>  ( t )] 

or  one  can  write  the  relationship  between  the  pair  of  functions  simply  as 

F  0 1 )  F(t). 

Whenever  it  is  obvious  from  the  context  which  version  of  the  ATT  is  being  used,  the 
subscript  on  the  operator  A  is  dropped. 

Figure  2.1  shows  the  relationship  between  the  Laplace  transform,  Fourier  trans¬ 
form  and  the  ATT.  The  H  operator  in  the  figure  represents  a  Hilbert  transform 
which  is  defined  shortly.  Notice  that  the  transformation  from  F  {t)  (Im  t  >  0)  to 
F{t)  -\-j'HF{t)  ( Im  t  =  0)  is  simply  evaluating  F  ( t )  on  the  real  time  axis  ( Im  t  =  0). 
But,  the  function  F  ( t )  may  not  be  analytic  for  Im  t  =  0.  In  fact,  F  (t)  at  ( Im  t  =  0) 
may  be  a  distribution  which  appears  quite  different  from  the  analytic  function  F  ( t ) 
at  Im  t  >  0.  An  important  example  of  this  is  the  analytic  delta  function  defined  by 


± 

7 r£ 

S(t)  +  p v  £ 


Im  t  >  0 
Im  t  =  0 


(2.12) 
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~  ~  ^  CD 

F(S)  - >  F(  CO)  - 

(Re  s  >  0)  (Im  co  =  0) 


+ 

>  F(t)  - ->F(l)  +  jflF(t) 

(Im  t  >  0)  (Im  t  =  0) 


Figure  2.1:  Diagram  showing  the  relationship  between  the  Laplace  transform  (LT), 
Fourier  transform  (FT)  and  analytic  time  transform  (ATT).  Notice  that  F(t)  can  be 
obtained  directly  from  an  inverse  Laplace  transform  only  if  F(t )  is  one-sided  since 
the  one-sided  Laplace  transform  is  used  here. 


where  “pv”  denotes  that  a  Cauchy  principal  value  is  understood  when  integrating  over 
this  function.  Notice  that  the  frequency  response  corresponding  to  the  analytic  delta 
function  is  F(uj)  —  1.  Also  notice  that  the  inverse  Fourier  transform  of  F( u>)  =  1, 
which  is  the  Dirac  delta  function  8(t),  is  properly  obtained  from  Re  [5  (t)]  when 
Im  t  —  0. 

An  obvious  yet  very  useful  observation  is  now  made.  The  ATT  defined  by  the 
integral  (2.6)  is  proportional  to  a  one-sided  Laplace  transform  of  the  frequency  domain 
function  F(u>)  where  the  frequency  variable  u>  acts  as  the  time  variable  in  (2.4)  and  the 
Laplace  variable  s  in  (2.4)  is  related  to  the  complex  time  variable  t  in  (2.6)  by  s  —  —jt. 
This  is  a  very  useful  observation  because  it  allows  one  to  use  the  extensive  tables  of 
one-sided  Laplace  transforms  that  exist  to  evaluate  the  ATT.  This  relationship  can 
be  written  symbolically  as 

F  (()  =  \C.  [#(*)]  I  (2.13) 

where  the  subscript  x  on  the  Laplace  operator  signifies  that  x  is  the  variable  of 
integration.  This  relationship  is  also  used  to  find  the  ATT  of  distributions  which 
are  not  typically  included  in  Laplace  transform  theory,  but  can  be  found  within  the 
context  of  distribution  theory  [23]. 
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Table  2.1:  Properties  of  the  analytic  time  transform  (ATT) 


i  S-o o  F(u;)u(u;)ejuitdw 

F(u) 

F(t) 

'J(i) 

F(aw),  a  >  0 

F  (t  -  a) 

Im  (t)  >  a  +  Im  (a) 

e-^aF{uj) 

ej0t  F  (t) 

0>  0 

F(u  -  0)u( u  -  0) 

ej0t  [F  (t)  +  l  !%  F^e^dw] 

F(u-/3) 

iF{t) 

juF{tjj) 

&F(t) 

(ju>)nF(w) 

~jt  F  ( t )  -  1F(0) 

in-) 

(-Jtr  F  ( t )  -  i(-it)"-'F( 0) - Ijr<"-»(0) 

Bn-) 

1  p  ((),  h(t)=\  F  (r)  G(t-  r)ir 

Im  t  >  a-f  +  ag 

Im  (£)  —  ag  >  e  >  ay 

F(u>)G(w) 

2.2  Properties  of  the  Analytic  Time  Transform 
(ATT) 

Some  important  properties  of  the  analytic  time  transform  (ATT)  are  presented  in 
this  section.  Also,  some  useful  special  properties  of  the  analytic  time  function  when 
it  is  evaluated  for  real  time  ( Im  t  =  0)  are  explored. 

Table  2.1  summarizes  the  properties  of  the  ATT.  Most  of  these  properties  can  be 
easily  derived  using  a  change  of  variables  as  in  (2.13)  and  a  good  reference  on  Laplace 
transforms  [26,  23]. 

2.2.1  Properties  of  Analytic  Time  Functions  for  Real  Time 

The  properties  of  the  analytic  time  functions  evaluated  on  the  real  time  axis  ( Im  t  = 
0)  are  now  examined.  During  this  discussion,  it  is  assumed  that  F(t)  and  F(ui) 
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Table  2.2:  Properties  of  the  Hilbert  transform.  The  Hilbert  transform  of  F(t )  is 
denoted  by  TiF[t)  or  F(t)  and  similarly  for  G(t). 


are  square  integrable  functions  unless  specified  otherwise.  Nonetheless,  many  of  the 
properties  discussed  here  can  be  easily  generalized  for  the  case  when  F{t )  and  F(u) 
are  distributions  [25,  23]. 

When  the  analytic  time  function  F  (£)  is  evaluated  on  the  real  time  axis,  Im  t  =  0, 
it  can  be  written  as 

F  {t)  =  F[t)  +  jHF(t)  for  Im  £  =  0  (2.14) 

where  F(t)  is  the  inverse  Fourier  transform  of  F(u) 

F(t)  =  F~'[F(u)}  (2.15) 

assuming  that  F(—  w)  —  F*( u>)  where  *  denotes  a  complex  conjugate.  The  H  operator 
denotes  a  Hilbert  transform 


HF{t)  =  T~x  [— jsgn(u;)F'(u;)] 


(2.16) 

(2.17) 


Table  2.2  shows  a  list  of  useful  properties  of  the  Hilbert  transform  [27].  Another 
notation  which  is  used  for  the  Hilbert  transform  is 


F(t)  =  HF{t) 
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The  following  Fourier  transform  pairs 


(2.18) 

(2.19) 

(2.20) 


F  ( t )  2F(u>)u(u) 

F(t)  +£+  F{u) 

F(t)  *-£->  -jsgn(u;)F(u>) 

(and  similarly  for  G  (t),  G(t)  and  G(t))  can  be  used  to  derive  some  additional  inter¬ 
esting  relationships,  such  as 


f°°  F(t)G{t)dt  =  l  Re  r  F  (0  G*  ( t)dt 

J~~  OO  2  J  — oo 


“  oo 

^F(t)*G{t)  =  F(t)*G(t) 
=  F[t)*G(t) 
l  Re  [F  (0*  G  (t)]  =  m  *  G(t) 


(2.21) 

(2.22) 

(2.23) 

(2.24) 


for  real  time,  Im  t  =  0,  and  *  here  denotes  a  convolution  along  the  real  time  axis. 


2.3  Some  Important  Analytic  Time  Functions 

Table  2.4  summarizes  some  useful  ATT  transform  pairs.  Most  of  the  transforms  in 
Table  2.4  can  be  derived  by  using  (2.13)  along  with  a  good  table  of  Laplace  transforms 
[28,  23].  Nonetheless,  some  care  must  be  taken  when  using  this  table,  especially  if  one 
is  applying  the  ATT  to  an  asymptotic  high  frequency  representation.  This  is  discussed 
in  more  detail  later  in  this  chapter.  If  one  is  working  with  a  Laplace  domain  function 
then  Table  2.3  can  be  used  to  find  the  corresponding  Fourier  domain  function  before 
using  Table  2.4.  Note  that  the  “Pf”  in  Table  2.3  stands  for  pseudofunction  as  defined 
by  Zemanian  [23].  When  one  integrates  over  a  pseudofunction,  one  takes  Hadamard’s 
finite  part  of  the  integral  which  may  be  thought  of  qualitatively  as  a  generalization 
of  Cauchy’s  principle  value  for  higher  order  singular  integrands.  For  the  Laplace 
domain  function  F(s)  =  s’1  where  rj  is  not  an  integer,  the  branch  cut  is  chosen  to  lie 
on  the  negative  real  s  axis  (  Re  s  <  0,  Im  s  =  0)  as  usual.  For  the  analytic  time 
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Table  2.3:  Relationship  between  some  Laplace  domain  functions,  F(s)  and  corre¬ 
sponding  frequency  (or  Fourier)  domain  functions,  F(w) 


FM 

s  =  a  +  ju 

(u>  is  real) 

1 

8 

TrS(u)  +  pV'jrj 

1 

3n  + 1 

n=0,l,2,.,. 

1 

P>-1 

Pf(JU))M+l 

^integer 

sn 

(M" 

n=0,l,2,... 

a” 

r(  ^integer 

a>0 

1 

ju+a, 

e— (t0-$)  e-*(*o+$) 

8  8 

Tejut°  sinc(^) 

A1 

o 

function  the  branch  cut  is  on  the  negative  real  axis  in  the  —jt  plane,  which 

corresponds  the  the  negative  imaginary  axis  in  the  t  plane  ( Re  t  =  0,  Im  t  <  0). 

The  notation  for  special  functions  in  Table  2.4,  such  as  Ei(z)  to  denote  the  ex¬ 
ponential  integral  or  erfc(z)  to  denote  the  complementary  error  function  follows  the 
notation  in  the  standard  reference  of  Abramowitz  and  Stegun  [29]. 

2.4  Representations  of  Analytic  Time  Functions 

This  section  is  concerned  with  various  representations  of  analytic  time  functions  which 
may  be  useful  for  both  analytical  and  computational  reasons.  A  late  time  series  in 
terms  of  inverse  powers  of  time,  (1/t),  is  derived  by  the  ATT  from  a  low  frequency 
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Table  2.4:  The  analytic  time  transform  (ATT)  of  some  useful  functions.  In  this  table, 
n  is  a  non-negative  integer,  fi  >  —1  is  not  an  integer,  and  77  is  not  an  integer. 


FM 

F(t) 

Im  t  >  0 

fXO  +  jW^t) 

Im  t  —  0 

1 

jL 

7T  t 

*(0  +  P  vti 

7T%)  +  pv4; 

I  _  larg(-jf) 

+  £M*l  +  7] 

wW  +  *tln  1*1  +  7] 

£  [|-iars  HO] 

£<*(*) 

|<|-^(n  +  l)] 

r(i+i)(1+icot*70  *>0 

jcSCTTH^^j  t<  0 

jnn! 

7r(-jOn+l 

«■«+» 

j”r(t,+i) 

7r(— ji)’)+1 

r(>>+i)  Aci*n  *  >  n 

-jT(r)-H)  #  <  fi 

1 

jw+a 

e  at[l  +  ^ij't(of)],  f  >  0 
— e"at^i(— at),  i<0 

£-,>/(^)erfc(^) 

e-«V(2<r2) 

+ip  vf^'r’dx 

Teju,t°  sinc(!f:) 

u(f  -  to  +  7)  -  «(t  -  <0  -  f ) 

,  ,•  ,  t-io+f 

+*ln 

e-“H 

a  >  0 

_ i _ 

7r(f+ja) 

ir(t2+a2)  +  ■?  ir(t2+a2) 
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power  series.  Then  an  early  time  series  in  terms  of  powers  of  time,  t ,  is  derived 
by  the  ATT  from  a  high  frequency  asymptotic  series.  It  is  shown,  by  an  example, 
that  an  early  time  series  obtained  by  a  term-by-term  application  of  the  ATT  may  be 
incomplete.  The  mathematical  justification  for  these  representations  may  be  found  in 
[26,  30,  31]  where  the  roles  of  time  and  frequency  are  switched  as  compared  to  here. 

2.4.1  Analytic  Time  Transform  of  Low  Frequency  Power 
Series 


The  ATT  of  a  series  which  is  accurate  at  low  frequencies  produces  a  time  represen¬ 
tation  which  is  asymptotically  valid  at  late  times.  The  mathematical  basis  for  this 
fact  is  well  known  and  can  be  found  in  texts  on  Laplace  transform  theory  [30,  26] 
where  the  connection  between  the  ATT  and  the  Laplace  transform  in  (2.13)  is  un¬ 
derstood.  In  particular,  the  general  theorem  on  the  Laplace  transform  of  asymptotic 
representations  is  on  page  6  in  Seigel  [32]. 

Consider  the  following  asymptotic  representation 

N 

F( u )  ~  ^2  An(ju>r  f°r  u  — *  0+  (2.25) 

n=l 

where  fin  >  —  1  and  fin  j  oo  as  n  — ►  oo.  Then  the  “early  time”  analytic  time 


representation  corresponding  to  this  asymptotic  low  frequency  representation  is 


for  M 


(2.26) 


where  Im  t  >  0. 


2.4.2  Analytic  Time  Transform  of  High  Frequency  Asymp¬ 
totic  Series 


When  solving  for  “early  time”  solutions  to  electromagnetic  radiation  or  scattering 
problems,  one  may  solve  the  frequency  domain  (FD)  (i.e.  time  harmonic)  problem 
first,  then  obtain  a  high  frequency  asymptotic  expansion  of  this  FD  solution  and 
finally  perform  an  inverse  Laplace  transform  to  obtain  a  time  domain  power  series 


[33,  34,  9].  Let  the  asymptotic  expansion  of  the  FD  solution  have  the  form 


Re  s  >  0 


(2.27) 


for  |s|  — »  oo  ,  where  s  =  a  +  jw,  Co  and  Cn  are  constant  with  respect  to  frequency, 
t0  is  some  time  delay  dependent  on  the  problem  geometry  and  vn  >  —1.  Note  that 
vn  is  not  necessarily  an  integer.  Also,  notice  that  N  in  (2.27)  is  not  the  same  as  the 
N  in  (2.26).  Typically,  an  inverse  Laplace  transform  is  applied  to  (2.27)  to  find  the 
corresponding  “early  time”  power  series  [33,  34] 

F(t)  ~  C0S(t  -  to)  +  E  r,  -"■:;■(*  -  toYnu{t  - 10)  (2.28) 

for  t  — >  0,  where  £(f)  is  the  Dirac  delta  function,  r(z)  is  the  Gamma  function,  and  u(t) 
is  the  Heaviside  unit  step  function.  This  time  domain  series  provides  a  satisfactory 
solution  if  the  constants  Co  and  { C„ }  are  real,  but  unfortunately  this  is  not  always 
the  case.  Consider  the  extremely  simple  case  of  the  plane  wave  reflection  from  a 
concave  surface.  The  first  term  in  the  high  frequency  asymptotic  expansion  is  the 
geometrical  optics  (GO)  reflected  field  which  may  pass  through  a  caustic  and  will  then 
gain  a  factor  of  eJir/2.  This  factor  makes  the  constant  Co  in  the  FD  asymptotic 
expansion  in  (2.27)  complex  and  then  the  time  domain  series  in  (2.28)  no  longer 
makes  sense  since  the  time  domain  solution  should  be  a  real  function. 

When  the  constants  {Co)Cn}  in  (2.27)  are  complex,  one  can  obtain  a  real  time 
domain  solution  F(t)  by  first  evaluating  jF'(.s)  for  Re  (a)  =  0  (<r  =  0)  to  obtain 
F(u>)  and  then  enforcing  the  relationship  F(-u)  =  F*(w)  for  negative  frequencies 
(u>  <  0)  and  applying  an  inverse  Fourier  transform,  where  the  frequency  variable  u 
and  the  time  variable  t  are  real.  For  example,  the  exp(j7r/2)  factor  discussed  above 
applies  only  for  positive  frequencies  and  more  generally  one  would  use  a  factor  of 
exp (ysgn (u> )7r /2 )  to  ensure  that  F{i)  is  a  real  function.  Notice  that  each  term  in 
(2.27)  is  not  absolutely  integrable  (when  Re  (s)  =  0),  and  therefore,  each  term  must 
be  considered  a  distribution,  see  Table  2.3,  and  one  must  resort  to  using  the  theory 
of  distributions  [23,  24]  or  generalized  functions  [35,  36]  to  find  the  inverse  Fourier 
transform. 

Alternatively,  one  could  obtain  a  real  time  domain  function  F[t)  from  the  Laplace 
domain  function  F(s)  by  enforcing  the  conjugate  symmetry  relationship  F(s*)  = 
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F*(s)  and  applying  an  inverse  two-sided  Laplace  transform.  But  this  can  be  quite 
complicated  even  for  a  simple  algebraic  F(s)  function. 

Instead  of  using  an  inverse  Fourier  transform  or  an  inverse  two-sided  Laplace 
transform  as  just  described,  it  may  be  easier  to  obtain  the  inverse  Fourier  transform 
via  the  analytic  signal  representation  described  earlier.  First  find  the  ATT  of  (2.27) 
to  provide  an  analytic  function  of  time.  Then  evaluate  this  analytic  function  for  real 
time  (Im  t  =  0)  and  take  the  real  part  to  obtain  the  final  real  time  function  which 
equals  the  inverse  Fourier  transform  [25,  37].  Note  that  the  analytic  time  function  is 
analytic  for  Im  t  >  0  and  its  boundary  value  at  Im  t  =  0  is  a  distribution  [25,  37]. 

The  form  of  the  early  time  asymptotic  expansion  obtained  from  the  high  frequency 
asymptotic  expansion  using  the  analytic  time  transform  depends  on  the  form  of  the 
high  frequency  expansion.  Three  particular  cases  are  presented  here.  The  following 
early  time  expansions  are  derived  by  using  the  results  in  Chapter  4  of  Bleistein  and 
Handelsman’s  book  [31]  and  the  relationship  between  Laplace  transforms  and  the 
analytic  time  transform  in  (2.13). 

Consider  the  following  high  frequency  asymptotic  expansion 

F( u)  ~  J  Am(ju>)~rm  for  W-+CO  (2.29) 

m=0 

where  rm  is  not  a  positive  integer  here  and  Re  (rin )  |  oo  as  m  — >  oo.  Then  the  early 
time  expansion  corresponding  to  the  function  in  (2.29)  obtained  from  the  asymptotic 
results  in  Bleistein  and  Handelsman  [31]  and  (2.13)  is 

F  (t)  ~  -  Arnj~Tm T(1  -  rrn)(-jt)Tm~l 

7rm= 0 

1  oo  (_  l'|m 

+-  E  (2.30) 

7r  __n  ml 

m~U 

where  M.[F\m  +  1]  is  the  generalized  Mellin  transform  of  F(ui)  which  is  defined 
shortly.  Notice  that  the  first  summation  in  (2.30)  is  exactly  what  would  be  obtained 
if  the  analytic  time  transform  in  (2.6)  is  applied  to  the  frequency  domain  asymptotic 
expansion  in  (2.29)  term-by-term.  Also,  notice  that  Tables  2.3  and  2.4  are  very 
useful  for  this  term-by-term  transformation.  But  the  second  summation  in  (2.30) 
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is  not  predicted  by  applying  the  analytic  time  transform  term-by-term  and  yet  it  is 
obviously  important.  In  fact,  the  example  in  the  next  section  demonstrates  this  fact. 

The  ordinary  (i.e.  not  generalized)  Mellin  transform,  denoted  by  the  operator 
symbol  M.  is  defined  as 

_  /*oo  _ 

M[F-,z}=  F(u)  u>z-'<Lj  (2.31) 

Jo 

The  generalized  Mellin  transform  is  defined  by 

M[F]z]  =  M[Fi;z]  +  M[Frtz\  (2.32) 

where 

M[Fx\z)  =  /^(wy-’dw  (2.33) 

Jo 

which  can  usually  be  analytically  extended  as  a  meromorphic  function  into  the  entire 
z-plane  and 

M[F2\z]  =  J |  F{u)uz~ldu  (2.34) 

which  can  usually  also  be  analytically  extended  as  a  meromorphic  function  into  the 
entire  z-plane.  So  the  expression  in  (2.32)  should  be  interpreted  as  a  meromorphic 
function  in  the  z-plane  even  though  the  integral  expressions  in  (2.33)  and  (2.34)  may 
not  converge  for  the  same  z  values.  See  Chapter  4  of  [31]  for  more  information  on 
the  generalized  Mellin  transform. 

Now  consider  the  following  possible  high  frequency  asymptotic  expansion 

OO 

F(uj)  ~  Am(ju})~Tm  for  u;  — >  oo  (2.35) 

m=0 

where  Re  (d)  >  0,  v  >  0  and  Re  (rm)  f  oo  as  m  -4  oo.  The  early  time  asymptotic 
expansion  corresponding  to  (2.35)  is 

F  (t)  -  -  E  ™  +  1  (2-36) 

7T  m_n  TO! 
m—U 

where  M.[F\  m+1]  is  the  generalized  Mellin  transform  of  F(w).  If  the  Mellin  transform 
M[F\m  +  1]  exists  in  the  ordinary  sense  (i.e.  the  integral  in  (2.31)  is  convergent)  for 
m  =  0,1,2,...  then  the  expansion  in  (2.36)  is  an  “expansion  by  moments”.  This 
expansion  is  used  in  Appendix  D  to  derive  (D.12). 


24 


Finally,  consider  a  third  possible  high  frequency  expansion  which  has  the  form 

OO 

F( u> )  ~  ^  Am{ju))~Tm  for  uj  — >  oo  (2.37) 

m= 0 

where  Re  (rm)  f  oo  as  m  — *  oo.  This  is  the  same  as  the  high  frequency  expansion  in 
(2.29)  and  for  the  terms  where  rm  is  not  an  integer  one  could  use  the  result  in  (2.30). 
However  if  rm  is  an  integer,  there  exists  a  different  corresponding  term  in  the 
early  time  series.  So,  consider  the  case  when  rm  is  an  integer  (say  rm  =  n-f  1  where  n 
is  an  integer),  then  the  corresponding  term  in  the  early  time  expansion  has  the  form 

-AJ--1  M-j'O  -  *(»  +  i)I 

7l\ 

+  (2.38) 

n\ 

where 

=  (2-39) 

and 

V>(1)  -  -7  = -0.577216...  (2.40) 

"  1 

tj){n  +  1)  =  ~7  +  7  ^or  71  —  1  (2-41) 

i=i  1 

where  7  is  Euler’s  constant.  For  example,  consider  the  following  high  frequency 
expansion 

00 

F(uj)  ~  y)  j4m(ju;)_m_1  for  u>  — ■>  oo  (2.42) 

m= 0 

where  m  is  an  integer,  then  the  corresponding  early  time  expansion  is 

+  1  00  f  (  — 

F(t)  ~  -  £  M- M-j'O  -  *(»»  + 1)) 

W  m=0  l  771  * 

+rm-1(^C7m(_J,r}  (2.43) 

Notice  that  if  one  assumes  that  the  expansion  in  (2.42)  can  be  transformed  term-by- 
term  and  uses  the  results  in  Table  2.4  for  each  term,  then  the  constant  Cm  in  (2.43) 
would  be  Cm  =  j  |  Am .  Further,  notice  that  if  the  one  sided  inverse  Laplace  transform 
is  applied  to  a  function  with  the  following  asymptotic  series,  which  has  a  finite  radius 
of  convergence  [30] , 

F(s)  ~  y)  for  |s|  — >  00  (2-44) 
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where  |arg(s)|  <  7r/2  ,  then  one  obtains  a  corresponding  early  time  expansion 

oo  A 

m  ~  E  “(')  <2-45) 

m=0  m' 

Now,  if  one  uses  Cm  —  j\Am  in  (2.43)  and  takes  the  real  part  of  F  ( t )  than  one 
recovers  the  result  in  (2.45).  It  is  now  reasonable  to  conjecture  that  the  formula  in 
(2.39)  for  Cn  could  be  simplified  to  Cn  =  y| Am  (where  rm  =  n  +  1)  although  this 
result  has  yet  to  be  shown. 


2.4.3  Example:  The  TD-UTD  Edge  Diffraction  Transition 
Function 

An  example  is  now  examined  to  demonstrate  some  of  the  formulas  presented  in  the 
previous  section.  The  particular  function  which  is  examined  is  the  TD-UTD  edge 
diffraction  transition  function  which  is  used  in  Chapter  3.  This  function  is  well 
understood  and  the  analytic  time  transform  is  known  in  closed  form.  So,  although 
one  would  normally  not  need  the  early  or  late  time  expansions  of  this  function  for 
computations,  it  is  an  interesting  example  since  the  early  and  late  time  expansions 
can  be  verified  by  examining  the  closed  form  time  domain  function. 

The  frequency  domain  function  of  interest  is 


F{ui )  =  y/iFxe’uxfce. rfc 

-  (~) 


ggjV4  [OC+jt  g -Ufi7/c 

- 7=~  /  —7 - —dfl 

y7T  J-oo+jt  fl2  +  JX 


(2.46) 

(2.47) 

(2.48) 


where  e  >  0  and  Ftr(z)  is  the  usual  frequency  domain  UTD  transition  function  defined 
in  (3.33).  For  now,  simply  assume  that  x  and  c  are  positive  constants.  Their  physical 
meaning  is  fully  explained  in  Chapter  3.  The  analytic  time  transform  of  (2.46)  can 
be  obtained  in  closed  form  as 


F(t) 


(2.49) 
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Now,  a  low  frequency  power  series  expansion  of  (2.46)  which  happens  to  be  con¬ 
vergent  for  all  finite  values  of  u>  is 


jp/  \  r~  iux/c  i  2  ^  (  l)n(ja;®/c)n+2 

1  ("J  =  v  1 - 7=  -  ./«,  V\ - 

V*  n=o  n\[2n  -f  1) 


(2.50) 


Using  (2.26),  one  can  write  down  the  corresponding  late  time  expansion  of  the  analytic 
time  transform  of  (2.50) 


F(t) 


x  y 


where 


\/jr  [—  jr  y/c(—  jr)3/2  2 c3/2  (— yr)5/2 


T  -t  +  x/c 


(2.51) 


Notice  that  (2.51)  is  convergent  for  \x/(ct)\  <  1.  Now,  the  late  time  expansion  in 
(2.51)  can  also  be  obtained  directly  from  the  closed  form  result  in  (2.49).  First  rewrite 
(2.49)  as 

>(()=  ^  x^e 


XUJ/C 


^ /n(-jr)  x/r(-jr )3/2 


and  then  use 


(1  +  z)-'»  =  1  -  \z  +  ^ 


-JTC 


for  |z|  <  1 


(2.52) 


to  obtain  (2.51).  Although,  using  the  closed  form  time  domain  function  to  obtain 
the  late  time  expansion  is  trivial  for  this  particular  function  one  should  keep  in  mind 
that  for  some  other  frequency  domain  functions  one  doesn’t  have  a  closed  form  time 
domain  formula  and  may  need  the  late  time  expansion  to  perform  computations. 

An  alternate  late  time  expansion  whose  convergence  is  centered  at  t  =  0  (instead 
of  t  =  —x/c  as  in  (2.51))  can  be  obtained.  To  obtain  this  alternate  expansion,  first 


substitute 


■^=£3(7) 


into  (2.50)  and  after  a  good  deal  of  algebra  obtain 


'<->  -  *=[■-£  (¥)”*(¥) 
4  fjux\3/2  1  /  ju>x\2  1 
3\/7F  \  c  /  2  \  c  / 


(2.53) 
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Now  use  (2.25)  and  (2.26)  to  obtain  the  alternate  late  time  expansion 


-j3'2 


(2.54) 


where  this  series  is  convergent  for  \x/(ct)\  <  1.  The  late  time  series  in  (2.54)  can  also 
be  obtained  directly  form  the  closed  form  time  function  in  (2.49).  Rewrite  (2.49)  as 


F(t) 


(2.55) 


and  then  use 


(1  +  z)  1  =  1  —  2  +  z2  —  z3  +  z4  •  •  •  for  \z\  <  1 


to  derive  the  late  time  series  in  (2.54). 

Now  an  early  time  expansion  is  derived  from  a  high  frequency  asymptotic  expan¬ 
sion  of  F(u)  in  (2.46),  which  can  be  written  as 


I'M 


C_\3/2  ±(±\“’ 

u )  4x2  Ijw ) 


15  c 


8a:3  \  jui 


2x  \juj ) 

7/2  105  /  c 

+  77  7  I  T- 


16a:4  V  ju) 


(2.56) 


for  u  — *  oo,  or  more  generally, 


F(u) 


~  .  .  1  •  3  •  •  •  (2m  -  1)  /  c 

’  2m  \iu>3 


(2.57) 


To  obtain  the  early  time  expansion  of  F  (£)  from  (2.57),  use  (2.30).  Notice  that  if 
the  ATT  is  applied  to  the  asymptotic  expansion  in  (2.57)  term-by-term  then  only 
the  first  summation  in  (2.30)  is  obtained.  But,  the  second  summation  is  obviously 
important  as  will  be  seen  shortly.  Using  (2.30)  on  (2.57),  the  first  summation  is 

S  (^“ ) 

which  is  the  result  which  could  be  obtained  by  a  term-by-term  transformation.  The 
second  summation  in  (2.30)  contains  the  Mellin  transform  of  the  frequency  domain 
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function  F(w).  This  can  be  evaluated  as  follows 


M[F\m  4- 1] 


roo  xe^!A  f<x>+jt  e-uiJ?/c 


—7=—  /  — - —da  umdu) 

V7T  J-oo+jt  flZ  +  JX 


poo  xeJ  1  poo-rjt 

JO  J—oo+jt  1 

xe^4  roo+jt  1 

\FH  J—oo+it  a2  + 


eJ  '  roo+jt  1  roo  , 

— -  /  - - -  /  ume~u^/cdu)  da 

y/ir  J-oo+jt  fi2  +  jx  Jo 

gi71"/4  fOO  -fje  1 

-7=-cm+1r(m  + 1)  /  .  ,  -m 

\/7T  ./-oo+ie  (jU2  + 


=  — m!\/7ra; 


-cm+1r(m  +  1)27T7 


»+ie  (n2  +jx){ii2)m+A 

1 _ 1 

2e-i«'Ay/x  (-jiz)m+1 


(2.58) 

(2.59) 
dp  (2.60) 

(2.61) 

(2.62) 


and  then  the  second  summation  in  (2.30)  is 


Finally,  the  total  early  time  expansion  can  now  be  written  as 


F(t) 


1  p  “  -jet 

V*y/=jiV  h>\  j* 


+j-r=E(-L 

V™  ^0  \3X 


H*r 


(2.63) 


Notice  that  this  series  is  convergent  for  \ct/x\  <  1.  If  the  time  variable  is  allowed  to 
be  real  Im  (<)  =  0  then  it  is  more  convenient  to  write  (2.63)  as 


+  ...  yfc  ^  /  ct\m  .  c  ^  /  ct\’ 

F  ~  S  ( * )  +J  v™  J?0  ( * ) 


(2.64) 


Recall  that  a  term-by-term  transformation  would  provide  the  first  summation  in 
(2.64)  and  therefore  for  this  example  a  term-by-term  transformation  correctly  predicts 
the  early  time  behavior  of  the  real  part  of  F  ( t )  when  evaluated  for  real  time  ( Im  (<)  = 
0)  but  the  imaginary  part  of  F  (t)  can  not  be  correctly  obtained  from  a  term-by-term 
transformation.  Now,  the  early  time  expansion  in  (2.63)  can  also  be  derived  directly 
from  the  closed  form  result  in  (2.49).  First,  write  the  closed  form  formula  in  (2.49) 


1  1  je 


(2.65) 


and  use 


(1  +  z)-1  =  5?(-l)mzm  for  |z|  <  1 
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to  obtain  the  early  time  expansion  in  (2.63). 

The  important  observation  to  be  made  from  this  example  is  that  one  can  not 
apply  the  analytic  time  transformation  to  a  high  frequency  asymptotic  series  term  by 
term.  Instead  one  must  use  a  more  careful  asymptotic  analysis  such  as  in  Bleistein 
and  Handelsman’s  book  [31]. 


2.5  Efficient  Time  Convolution 


Consider  the  problem  of  finding  the  ATT  of  the  product  of  two  frequency  domain 
functions  F(u)H(w).  The  time  convolution  property  can  be  used  here 


F{u)H(u)  ^\f  (0*  H  (t) 


(2.66) 


where  the  time  convolution  integration  path  is  placed  above  the  real  time  axis  in  the 
complex  time  plane.  If  the  frequency  domain  function  F(u)  can  be  modeled  using 


exponential  functions 


or  equivalently,  if 


F(uj)  =  ^2  Ane  a”u  for  u  >  0 


+  i  N  A 

ht)  =  Lyj-^— 

K)  *^t  +  jan 


=  YjAn  5  (t  +  jan) 


(2.67) 


(2.68) 


(2.69) 


then  the  convolution  can  be  evaluated  in  closed  form  as 

1  +  +  N  + 
r  F  (<)*  H  (0  =  £  An  H  (t  +  jan) 

n=l 

where  the  shifting  property  of  the  analytic  delta  function 

(t  +  j/3)  =H  (t+j/3) 


was  used. 


(2.70) 


(2.71) 


A  useful  approximate  formula  for  the  convolution  operation  may  be  obtained 
when  the  frequency  response  of  the  excitation  is  bandlimited  (i.e.  .Fq(u;)u(u;)  =  0 
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when  u;  >  u>m  for  some  constant  u>m).  This  useful  approximate  formula  is  derived 

+ 

with  the  knowledge  that  Fq  (t)  is  analytic  in  the  entire  t-plane. 

\  h  (0*  S(f)  =  \  r+”  ^  (*-*)&  (“)*»  (2-72) 

Z  Z  — oo+je 

=  \r  3  -  J>)  B  (0  +  }i)dff  (2.73) 

Z  J —oc 

1  r°°  +.  + 

=  2  J-oc  F°  ^  ~  ^  n  ^  "  7  +  ;e^7  (2‘74) 

d  M  + .  + 

~  -  X]  -Po  (md-je)  H  (t  -  md  +  je)  (2.75) 

2  m=-M 

+  . 

The  above  approximation  is  obviously  the  same  as  (2.70)  with  ,An  =  |  Fp  (md  —  je), 
an  —  e  +  jmd  and  m  =  n  —  N  +  1. 
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Chapter  3 

TD-UTD  for  a  Curved  Wedge 


This  chapter  presents  the  development  of  the  TD-UTD  solution  for  the  scattering 
from  a  perfectly  electrically  conducting  (PEC)  arbitrary  curved  wedge  when  it  is 
illuminated  by  an  analytic  time  impulsive  astigmatic  wavefront.  This  solution  is 
obtained  by  an  analytic  time  transform  (ATT),  as  defined  in  the  previous  chapter,  of 
the  corresponding  frequency  domain  UTD  solution  for  a  perfectly  conducting  curved 
wedge  [2,  4,  38].  Thus,  the  present  solution  constitutes  a  time  domain  UTD  (TD- 
UTD)  for  a  curved  wedge  in  which  the  transient  fields  propagate  along  the  ray  paths 
of  the  UTD.  The  TD  ray  solution  (or  TD-UTD)  is  only  valid  for  “early  to  intermediate 
times”  since  the  UTD  is  an  asymptotic  solution  in  the  frequency  domain  that  remains 
accurate  for  moderate  to  high  frequencies.  Thus,  each  TD-UTD  ray  contribution  is 
most  accurate  in  the  neighborhood  of  the  time  required  to  traverse  its  geometric 
ray  path  length  from  the  source  to  the  observer  (i.e.  the  ray  arrival  time).  It  is 
noted  that  the  transient  response  in  intervals  close  to  ray  arrival  times  of  the  TD- 
UTD  are  generally  the  strongest  and  are  directly  associated  with  the  waves  arriving 
from  the  local  features  of  the  radiating  (scattering)  object.  The  TD-UTD' has  the 
same  geometrical  ray  paths  as  the  frequency  domain  UTD  except  that  each  incident, 
reflected  or  diffracted  TD  ray  is  associated  with  its  corresponding  progressing  wave 
in  space-time.  The  TD-UTD  fields  are  polarized  transverse  to  their  ray  paths  as  in 
the  conventional  frequency  domain  UTD. 

There  have  been  many  investigations  on  the  TD  diffraction  by  a  straight  wedge. 
Keller  and  Blank  solve  the  problem  of  electromagnetic  diffraction  and  reflection  of 
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plane  wave  pulses  by  perfectly  conducting  wedges  using  the  conical  flow  method  [39] . 
Friedlander  solves  the  corresponding  acoustic  problem  of  wedge  diffraction  [7].  Felsen 
derives  the  scattering  from  a  straight  wedge  (and  some  other  canonical  geometries) 
with  various  types  of  illumination  [34,  9,  22].  Nikoskinen  et  al.  use  a  transient 
image  theory  to  obtain  the  transient  scattering  from  a  PEC  wedge  [40].  Ianconescu 
and  Heyman  solve  the  pulsed  electromagnetic  field  diffraction  by  a  straight  wedge 
using  the  Spectral  Theory  of  Transients  (STT)  [19,  14,  15,  16].  Using  the  STT, 
they  also  extend  their  results  for  an  incident  field  which  is  an  astigmatic  collimated 
wavepacket  [18,  17].  The  work  by  Ianconescu  and  Heyman  is  relevant  to  the  work 
in  this  report  since  they  use  an  analytic  time  function  representation  similar  to  the 
one  used  here.  Veruttipong  and  Kouyoumjian  obtain  a  time  domain  version  of  the 
UTD  for  a  straight  perfectly  conducting  wedge  [12,  11,  10]  by  applying  an  inverse 
Laplace  transform  to  the  corresponding  frequency  domain  UTD  wedge  diffraction 
coefficient  [2],  Unfortunately,  Veruttipong’s  solution  is  not  easily  generalized  to  the 
curved  wedge  case-where  the  reflected  or  diffracted  rays  may  pass  through  smooth 
or  line  caustics. 

The  TD-UTD  solution  for  a  perfectly  conducting  wedge  presented  in  this  chapter 
differs  from  previous  solutions  in  that  the  present  solution  is  applicable  to  the  more 
general  problem  of  diffraction  from  an  arbitrary  curved  wedge  which  may  contain 
curved  faces  and/or  a  curved  edge,  while  the  incident  time  dependent  wavefront  is 
allowed  to  be  astigmatic.  It  is  noted  that  plane,  cylindrical,  conical  and  spherical 
wavefronts  are  all  special  cases  of  the  general  astigmatic  wavefront  which  is  char¬ 
acterized  by  two  distinct  radii  of  curvature.  In  fact,  Veruttipong’s  TD  diffraction 
coefficient  for  the  straight  wedge  [12]  can  be  found  as  a  special  case  of  the  solution  in 
this  chapter.  In  addition,  this  work  employs  an  analytic  time  function  representation 
for  the  transient  fields  where  the  actual  time  response  is  simply  the  real  part  of  the 
analytic  function,  as  discussed  in  the  previous  chapter.  It  is  shown  in  this  TD-UTD 
development  that  the  use  of  analytic  time  functions  allows  one  to  circumvent  in  a 
simple  manner  some  complexities  that  arise  when  inverting  into  the  TD  the  UTD  ray 
fields  which  pass  through  caustics.  These  caustics  may  be  produced  by  reflection  and 
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edge  diffraction  from  general  curved  geometries.  The  analytic  time  functions  used 
here  are  closely  related  to  the  analytic  time  functions  used  in  the  STT  [14,  15,  16]. 
The  STT  can  in  general  provide  an  analytical  TD  impulse  response  which  is  valid  at 
all  times  for  some  special  canonical  geometries,  but  the  TD-UTD  which  is  valid  for 
early  to  intermediate  times  is  immediately  applicable  to  more  general  and  relatively 
complex  geometries  due  to  its  ray  interpretation. 

The  format  of  this  chapter  is  as  follows.  First  the  well  known  frequency  domain 
UTD  solution  for  the  curved  wedge  is  briefly  presented.  In  the  UTD,  the  dominant 
response  from  the  PEC  curved  wedge  is  decomposed  into  the  incident  field,  the  re¬ 
flected  field  and  the  edge  diffracted  field.  The  edge  diffracted  field  may  contain  two 
terms,  one  is  the  first  order  diffracted  field  while  the  other  is  a  higher  order  slope 
diffracted  field.  The  slope  diffracted  field  is  significant  when  the  field  incident  at  the 
point  of  diffraction  is  rapidly  varying  as  a  function  of  space.  Then  the  time  domain 
analytic  signal  representation  for  each  of  these  mechanisms  is  derived  to  obtain  a 
TD-UTD  for  the  curved  wedge.  This  chapter  does  not  include  an  analysis  of  the 
surface  rays  excited  and  diffracted  by  the  edge  of  a  wedge  with  convex  faces,  nor  does 
it  include  the  diffraction  and  excitation  of  whispering  gallery  type  fields  by  the  edge 
of  a  wedge  with  concave  faces.  Nonetheless,  some  curved  surface  diffraction  effects 
for  smooth  surfaces  are  discussed  in  Chapter  IV.  The  simplification  of  the  TD-UTD 
wedge  diffraction  coefficient  for  some  important  special  cases  is  also  presented.  The 
convolution  of  the  TD-UTD  curved  wedge  impulse  response  with  a  more  general  tran¬ 
sient  excitation  is  discussed.  Some  numerical  examples  are  presented  to  illustrate  the 
utility  of  these  developments. 

3.1  The  Frequency  Domain  UTD  Solution  for 
Scattering  from  a  Curved  Wedge 

The  general  forms  in  the  frequency  domain  of  the  incident,  reflected  and  edge 
diffracted  UTD  rays,  illustrated  in  Figures  3.1,  3.2  and  3.3  respectively,  are  now 
presented  [2,  4,  38]. 
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Figure  3.1:  An  incident  geometrical  optics  (GO)  ray  tube. 

3.1.1  The  Frequency  Domain  Geometrical  Optics  (GO) 
Fields 

The  incident  ray  shown  in  Figure  3.1  is  a  diverging  (i.e.  convex)  wavefront  but  one 
could  also  allow  the  wavefront  to  be  convergent  (i.e.  concave)  or  even  saddle  shaped. 
The  incident  field  at  the  point  of  reflection  or  point  of  diffraction  is  assumed  to  be 
ray  optical  and  may  be  written  as  a  geometrical  optics  (GO)  ray  field  such  that 

E’M  =  E‘H  |A,(s‘)|  jnie-jk9<  (3.1) 

where  n'  is  the  number  of  line  or  smooth  caustics  that  the  incident  GO  ray  has 
traversed  as  in  Figure  3.1  and  A, -(a*)  is  the  spreading  (or  spatial  divergence)  factor 
which  may  be  written  as 


M*)  = 


P\  \  (  P2 

Pl+SJ  \p*2  +  5* 


Note  that  p\  and  p'2  are  the  principal  radii  of  curvature  of  the  corresponding  wavefront 
at  the  reference  point  0.  A  GO  field  reflected  from  the  point  Qr  as  in  Figure  3.2 
may  be  written  as 

E»  =  .  it  \AT(sr)\jnre~ik"T  (3.3) 

where  the  dyadic  reflection  coefficient  R  is  assumed  to  be  independent  of  frequency 
and  can  be  written  as 


R  —  efienifo  +  cj_ej_jRs 


(3.4) 


where  Ra  =  —  1  and  Rh  —  1  for  a  perfect  electric  conductor.  The  .Ar(sr)  is  the 
spreading  factor  for  the  reflected  ray  which  is  given  by  (3.2)  with  the  superscript  i 
replaced  by  r,  and  nr  is  the  number  of  smooth  or  line  caustics  that  the  reflected  ray 
has  traversed  after  leaving  the  reflection  point  Qr.  It  is  assumed  in  the  UTD  that 
(3.3)  represents  the  field  reflected  from  the  unperturbed  surface  formed  by  a  smooth 
extension  of  the  original  surface  past  the  edge  (so  that  there  is  no  edge  effect  present 
in  the  evaluation  of  (3.3)  as  in  Figure  3.2). 


3.1.2  The  Frequency  Domain  UTD  Edge  Diffracted  Field 

The  edge  diffracted  field  may  be  modeled  by  an  edge-diffracted  ray  as 


EdH  =  E‘H|^-f>(u;)|^d(sd)|i 


;nde~jksd 


(3.5) 


where  Ad(sd)  is  the  spreading  factor  for  the  edge  diffracted  ray,  D  is  the  UTD  dyadic 
diffraction  coefficient  [4,  2]  and  nd  is  the  number  of  line  or  smooth  caustics  that  the 
edge  diffracted  ray  has  traversed  after  leaving  the  diffraction  point  Qe •  Note  that 
Equation  (3.5)  does  not  include  slope  diffraction  [38,  4]. 


3.1.3  The  Frequency  Domain  UTD  Slope  Edge  Diffraction 

The  UTD  Slope  diffraction  is  a  higher  order  term  which  must  be  added  to  the  edge 
diffraction  term  described  above  when  the  incident  field  at  the  point  of  diffraction  is 
rapidly  varying  in  a  direction  transverse  to  the  direction  of  incidence.  In  the  frequency 
domain,  the  slope  diffraction  is  proportional  to  a  spatial  derivative  of  the  incident  field 
at  the  point  of  edge  diffraction.  Two  versions  of  the  slope  diffraction  are  used  here. 
The  first  one  is  by  Y.  M.  Hwang  [21,  4].  While  the  other  one  is  by  Veruttipong  [20]. 
The  slope  diffraction  by  Hwang  is  more  widely  used  than  Veruttipong’s  formulation 
since  Hwang’s  result  has  a  more  compact  form.  Nonetheless,  it  is  interesting  to 
compare  these  two  solutions  in  the  time  domain  as  is  done  later  in  this  chapter. 

The  UTD  slope  diffracted  field  for  both  versions  has  the  same  form,  only  the 
slope  diffraction  coefficient  is  different.  In  general,  one  can  write  the  total  UTD  edge 


37 


X/(2-n)7TV\ 

0  =  n7i  0  =  0 
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AT  Qe 


Figure  3.3:  Coordinate  systems  for  the  TD-UTD  dyadic  diffraction  coefficient 


diffracted  field  as  the  sum  of  two  terms 


EdV)  =  Ed(u)  +  ESV) 


(3.6) 


where  the  first  order  term  Ed(u;)  is  given  in  (3.5)  and  Esd(u>)  is  the  higher  order 
slope  diffraction  term.  The  slope  diffracted  field  for  a  curved  wedge  can  be  written 
as  Ead  =  j30Ef0  +  <\>Ef  with 


.  dE\ 


dn • 


-  BEZ, 


ir9Jk 
Vh  Unr 


M*d)  jnde-jk°d. 


The  extra  superscript  s  on  the  diffraction  coefficients  indicates  slope  diffraction.  The 
result  in  (3.7)  assumes  that  the  incident  and  reflected  field  vary  rapidly  at  the  edge 
in  the  n*  and  the  nT  directions,  respectively.  On  the  other  hand,  if  the  diffracted  field 
spreading  facter  ^(.s'*)  is  a  rapidly  varying  function  then  the  result  in  (3.7)  is  not 
accurate. 

For  a  straight  wedge,  this  can  be  simplified  to 


.  dEL 


h  dn * 


Msd)\j 


nde~jksd 


where  Dss  —  DSI  +  Dssr  and  Dsh  =  Dsl  -f  .  The  expressions  for  the  slope  diffraction 
coefficient  for  both  Hwang’s  and  Veruttipong’s  versions  are  given  in  Appendix  C 
where  the  time  domain  version  of  these  quantities  are  derived. 

For  more  information  about  the  parameters  in  the  UTD  ray  expressions,  one 
may  consult  the  literature  [2,  38,  4].  It  should  be  noted  that  there  are  also  edge 
excited  surface  diffracted  rays  or  whispering  gallery  effects  present  for  a  curved  wedge 
geometry;  however,  the  present  chapter  does  not  deal  with  these  surface  diffracted 
and  whispering  gallery  type  ray  fields.  Surface  diffraction  by  smooth  convex  surfaces 
is  addressed  in  Chapter  4. 
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3.2  The  TD-UTD  Impulse  Response  for  a  PEC 
Curved  Wedge 

The  TD-UTD  analytic  impulse  response  for  a  PEC  Curved  Wedge  (not  including 
surface  diffracted  rays  and  whispering  gallery  effects)  may  be  written  as 

E 7°  (t)  =E/  (t)Pi+  E';  (t)Ur+  E?  (1)+  E?  (()  (3.9) 

where  E/  (t),  E/  (t) ,  E/  (f)  and  E/d  ( t )  are  the  analytic  signal  representations  for 
the  incident  field,  the  reflected  field,  the  first  order  edge  diffracted  field  and  the 
slope  edge  diffracted  field,  respectively.  The  field  in  (3.9)  results  from  an  impulsive 
transient  illumination  of  the  wedge  as  denoted  by  the  subscript  I.  The  analytic 
impulse  response  is  defined  as  the  response  which  results  when  the  incident  field 
E/  ( t )  has  an  analytic  delta  function,  8  (t)5  time  behavior.  The  spatial  unit  step 
functions  £/,  and  UT  are  1  on  the  lit  side  of  the  incident  shadow  boundary  and  the 
reflection  shadow  boundary,  respectively,  and  0  otherwise. 

Specifically,  it  is  noted  that  the  TD-UTD  analytic  impulse  response  of  the  curved 
wedge  is  the  response  due  to  an  excitation  of  an  astigmatic  wavefront  with  an  analytic 
delta  function  time  behavior.  This  analytic  impulse  response  is  found  by  setting 
E*,(w)  =  Ej,  in  (3.1),  (3.3),  and  (3.5),  where  Ej,  is  a  complex  constant  with  respect  to 
frequency,  and  then  applying  the  analytic  time  transform  in  (2.6)  to  these  equations. 
Finally,  these  results  are  combined  as  shown  in  (3.9).  This  procedure  provides  analytic 
functions  of  time  where  the  actual  real  time  responses  are  found  from  the  real  part 
of  these  functions  evaluated  for  real  time,  as  discussed  in  Chapter  2.  The  response  to 
any  desired  finite  energy  pulsed  excitation  is  obtained  by  performing  a  convolution 
of  the  finite  pulsed  field  excitation  with  the  TD-UTD  analytic  impulse  response. 


3.2.1  The  TD-UTD  Geometrical  Optics  Field 


The  incident  GO  ray  field  with  an  analytic  impulse  time  behavior  is  obtained  from 
(3.1)  and  (2.6) 

E/  (t)  =  E U,(s*')|  8  (t  -  s*/c)  (3.10) 
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with  the  understanding  that  the  actual  real  time  incident  field  E}(f)  is  given  by 

E )(t)  =  Re  [E/  (<)]  for  Im  t  =  0  (3.11) 

It  is  important  to  note  that  in  general  Eq  is  a  complex  vector  which  is  constant  with 
respect  to  frequency  for  the  analytic  impulse  response  and  n,  may  not  be  zero;  it 
therefore  becomes  necessary  to  define  the  incident  field  excitation  more  precisely  as 
a  generalized  impulse  (or  analytic  delta  function),  8  (t).  As  long  as  the  analytic 
signal  representation  of  the  incident  field  is  a  complex  constant  times  the  analytic 
delta  function,  then  the  transient  response  is  called  the  “impulse  response”.  This  is 
consistent  when  using  the  analytic  time  signals,  since  any  transient  response  due  to 
any  other  excitation  can  be  found  by  convolution.  The  analytic  delta  function  8  ( t ) 
was  discussed  earlier  and  is  given  in  (2.12).  One  of  the  benefits  of  using  the  analytic 
time  representation,  which  is  now  clearly  evident  from  equations  (3.10)  and  (2.12),  is 
that  the  analytic  time  function  automatically  represents  the  TD  ray  field  before  and 
after  the  ray  has  traversed  a  caustic  or  any  number  of  caustics.  For  example,  if  the 
incident  ray  has  not  traversed  any  caustics  then  n,  =  0  and,  assuming  Ef,  is  real,  the 
real  time  incident  field  behavior  is 

E  )(t)  =  Ej  \Mj)\ S(t  -  a'/c)  (3.12) 

but  if  the  incident  ray  has  traversed  one  smooth  or  line  caustic  then  n,  =  1  and  the 
real  time  behavior,  assuming  EJ,  is  real,  is 

Ej(t)  =  — l (3.13) 

When  the  analytic  signal  representation  in  (3.10)  is  used,  one  does  not  have  to  be 
concerned  about  the  number  of  caustics  that  the  ray  has  traversed,  and  it  avoids  the 
need  to  do  a  separate  analysis  for  each  case;  this  is  especially  important  in  inverting 
the  diffracted  field  Ed(u;)  when  a  ray  caustic  is  traversed  making  the  analysis  of 
this  case  rather  complicated  to  get  E j(t)  if  the  analytic  signal  representation  is  not 
used.  Note  that  the  ray  field  description  (3.10)  breaks  down  if  the  observation  point 
is  located  near  a  caustic.  The  reflected  ray  contribution  to  the  TD-UTD  impulse 
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response  is 


-  E;  (t)  =  E‘  •  Rf  <+"r  |A(S’K(^r)|  6  (Tr)  (3.14) 

where  Tr  —  t  —  s'/c  —  sr /c.  Note  that  the  variable  5*  in  (3.14)  is  the  distance  from  the 
source  to  the  point  of  reflection  Qr\  whereas,  s'  in  (3.10)  is  the  straight  line  distance 
from  the  source  to  the  observer.  Also,  note  that  the  convolution 

\  £  (f  -  s'  j c)*  6(t~  *7c)  =6  (t  ~  j/c  ~  ^/c)  (3.15) 

is  used  to  obtain  (3.14). 

3.2.2  The  TD-UTD  Edge  Diffracted  Field 

The  first  order  edge  diffracted  ray  contribution  to  the  TD-UTD  impulse  response  is 

Ei  ( t )  =  E*  •  D  0 rd)jni+n d  |a,(5*')^(^)|  (3-16) 

where  tj  —  t  —  s'/c  —  sd/c,  and  s'  in  (3.16)  denotes  the  distance  from  the  source  to 
the  point  of  edge  diffraction  Q /■; .  The  TD-UTD  dyadic  diffraction  coefficient  for  a 
PEC  curved  wedge  is 

D  (t)  =  -/3o/30  Ds  (0  -  44>  Dh  (i t )  (3.17) 

where  the  unit  vectors  in  (3.17)  are  depicted  in  Figure  3.3. 


Derivation  of  the  TD-UTD  Wedge  Diffraction  Coefficient 

The  derivation  of  the  TD-UTD  dyadic  diffraction  coefficient  for  a  PEC  curved  wedge 
is  now  presented.  The  corresponding  time  harmonic  dyadic  diffraction  coefficient  [2] 
can  be  written  in  the  following  form  [11] 

D(u;)  =  -  $4>DkW  (3-18) 
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and  a±(/3)  =  2  cos2[(2n7rTV±  —  /3)/2]  where  jV*  is  the  nearest  integer  solution  of 
2irnN±  —  (3  =  ±7r.  For  a  thorough  explanation  of  the  parameters  used  in  the  UTD 
edge  diffraction  coefficient,  one  is  referred  to  the  literature  [2,  38,  4].„  All  of  the 
frequency  dependence  has  been  collected  in  the  function 

Or1)  <3-32> 

where  Ftr(z)  is  the  conventional  UTD  transition  function  defined  by 

Ftr(z)  =  2 r  e~**dT  (3.33) 

•V* 
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The  F(xm,u>)  function  can  also  be  written  as 

F(xm,  u)  =  v/’f*meJU':rm/cerfc 


where  the  complement  of  the  error  function  is  defined  as 

2  r°°  ,1 

erfc(z)  =  — ■=  /  e  1  dt 

\/TT  Jz 


(3.34) 


(3.35) 


Note  that  the  function  F(xm,w)  is  not  the  same  as  the  transition  function  defined  in 

[2]. 

The  analytic  time  transform  of  the  frequency  domain  function  F(xm,u ;)  is  now 
derived.  The  following  integral  (which  may  be  found  in  a  table  of  Laplace  transforms 
[28])  is  used 

/“  ettTeitc(y/ar)e~FTdT  =  —  1  (3.36) 

Jo  #^  +  V“) 

where  —tt  <  arg(p)  <  7r  and  —ir  <  arg(a)  <  7r.  Note  that  — 7r  <  arg (jxm)  <  ir  and 
therefore  -Zir/2  <  arg(a;m)  <  7t/2.  Now,  using  the  integral  in  (3.36),  the  analytic 
signal  corresponding  to  the  frequency  domain  function  F(xm,w)  is 


F  (®m)  ^)  — 


yf-ftW  7^  "b  \J 3xm! c) 


(3.37) 


for  Im  t  >  0  where  —  n  <  arg(-jt)  <  7r.  Notice  that  the  singularity  on  the  real  time 
axis  ( Im  t  =  0)  is  only  a  branch  point  singularity  and  not  a  pole,  and  that  there  is 
no  difficulty  in  using  the  expression  for  F  (x, m,t)  in  (3.37)  for  real  time  (Im  t  =  0). 
Paying  close  attention  to  the  Riemann  sheet  on  which  the  parameter  xm  is  defined 
(i.e.  — 37t/2  <  arg(zm)  <  7t/2),  (3.37)  may  be  recast  into  a  more  convenient  form 


F  (®m  j  t'j  — 


_ -iy-gm/g _ 

+  e-W4v/-Xm/c) 


(3.38) 


for  Im  t  >  0  where 


v—=  I  *">0  (3.39) 

yj  ®m  5  ^  0 

and  Re  (\/—jt)  >  0.  The  form  in  (3.38)  is  more  convenient  for  numerical  computa¬ 
tions  since  most  calculators  and  computer  programming  languages  choose  the  branch 
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cut  of  the  square  root  function  to  be  the  negative  real  axis  and  then  (3.38)  may  be 
directly  implemented  without  any  ”if-then”  conditions.  Finally,  the  TD-UTD  dyadic 
diffraction  coefficient  is  given  by 


D  (f)  =  — /3q/30  Ds  ( t )  -  $4  Dh  ( t ) 


(3.40) 


where 


Ds,H  (t)  =  -  W  ,  £  F 

2n V27T  sm  p0  m=i 


(3.41) 


The  dyadic  diffraction  coefficient  in  (3.17)  can  be  applied  to  the  general  geometry 
of  a  PEC  curved  wedge  which  may  have  curved  faces  and  a  curved  edge.  Nevertheless, 
it  is  interesting  to  examine  some  special  cases  such  as  a  straight  wedge  and  a  straight 
half  plane. 


Simplification  of  TD-UTD  Diffraction  Coefficient  for  Some  Important  Spe¬ 
cial  Cases 

First,  the  behavior  of  the  diffraction  coefficient  near  a  shadow  boundary  (SB)  is 
investigated,  then  the  special  case  geometries  of  a  straight  wedge  and  flat  but  possibly 
curved  screen  are  examined. 

There  are  four  possible  SB  locations  corresponding  to  the  four  terms  in  (3.41). 
Define  four  corresponding  angles 


*  +  (<t>  -  <t>) 

<t>') 

7T  +  (^  -f  <f>)  —  2n7T 
*-{<!>  + 


(3.42) 

(3.43) 

(3.44) 

(3.45) 


When 


0+  the  observer  is  approaching  the  mth  SB  from  the  lit  side.  When 


0  then  xm  — ►  0  and  it  can  be  shown 


KmhV^ r  ~  sgn(em) 


(3.46) 


where  C^h  =  1  if  m  =  1,2  and  C^h  =  q=l  if  m  =  3,4.  Using  (3.46),  it  can  be 
shown  that  the  TD-UTD  wedge  solution  has  the  same  behavior  near  a  SB  as  the 
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exact  solution  in  [19]  (where  A,  in  Equation  (27)  in  [19]  is  related  to  X2  in  this  report 
by  A ,  «  yfx 2/  c  near  the  SB).  By  using 

F  ( Xm,t )  ~  y/xmv  6  ( t )  for  xm  *  0  (3.47) 


with  (3.46),  it  can  be  shown  that  the  TD-UTD  diffracted  field  has  the  proper  dis¬ 
continuity  which  when  added  to  the  GO  field  produces  a  continuous  total  field  across 
the  SB.  This  is  expected  since  the  same  behavior  across  a  SB  occurs  in  the  frequency 
domain  UTD. 

Assume  the  two  faces  of  the  wedge  are  flat  and  the  incident  GO  field  is  a  real 
impulse  in  time  (i.e.  5(f))  so  that  all  of  the  L  parameters  are  positive  and  therefore 
xm  >  0.  Now,  the  analytic  time  function  in  (3.38)  may  be  simplified  to 


*)  = 


-/t(t  +  xm/c) 


(3.48) 


for  Im  t  >  0,  and  the  corresponding  real  time  function  is 


,*-■*)-  [*(-•)]  =  (3'49) 
at  Im  t  =  0  where  u(t)  is  the  Heaviside  unit  step  function.  Equation  (3.49)  along  with 
(3.41)  is  exactly  the  TD-UTD  diffraction  coefficient  derived  by  Veruttipong  [12],  which 
therefore  is  a  special  case  of  the  more  general  TD-UTD  diffraction  coefficient  presented 
in  this  report.  The  solution  may  be  further  simplified  by  letting  |xm/(ct)|  — *  00,  which 
is  justified  when  the  observation  point  is  far  away  from  any  SB.  When  this  is  done 
the  result  in  (3.49)  may  be  further  simplified  to 

*X»)  =  .  (3.50) 

y/irt 

Equation  (3.50)  together  with  (3.41)  is  what  one  obtains  by  taking  the  inverse  Fourier 
transform  of  Keller’s  frequency  domain  GTD  diffraction  coefficient  [12]. 

The  TD-UTD  diffraction  coefficient  for  a  flat  screen  (which  may  be  curved)  is 


found  by  setting  n  =  2  in  (3.41)  to  obtain 

+  /,x  —  1  F  (*i4>0  F  (xBjt) 

Da'h  ^  ~  2^  sin  &  [cos  (*=£)  ^  cos  (**), 


(3.51) 
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where 


xA  =  2  V  cos2  {^2^^ 

XB  =  2L'cos>(tt£\ 

and 

r  =  -^-dA  »*A 

/9*  +  3d 

rr  _  &  P  *  2  a 

L  -  d  sin  (Jo 

pr  +  sd 

3.2.3  The  TD-UTD  Slope  Diffraction 


(3.52) 

(3.53) 

(3.54) 

(3.55) 


The  analytic  time  transform  (ATT)  can  be  applied  to  the  higher  order  frequency 
domain  UTD  slope  diffraction  term  in  (3.7)  to  obtain  a  TD-UTD  slope  diffraction  in 
closed  form. 


The  TD  response  corresponding  to  slope  diffraction  is 


o/f  A  ■(■  o (J  a  -j-  „J 

with  E/  (t)  =  /30  EpQ  (tj)  +  (f>  E$  (tj)-  The  *  denotes  a  TD  convolution.  The 
analytic  impulse  response  is  obtained  by  letting  the  incident  field  (and  therefore  also 
the  reflected  field)  in  (3.56)  be  impulsive  in  time. 

The  TD-UTD  slope  diffraction  coefficient  obtained  from  Hwang’s  formulation  for 
a  straight  wedge  is  derived  in  Appendix  C.  Here  is  the  same  coefficient  •■generalized 
for  a  curved  wedge 


Dsi  (t)  = 


4n2\/27r  sin2  /3C 


[Cl  Fa  (*ljf)  ~  C*2  Fa  (*2j0] 


&  W  =  +F‘  M ' 

'  If  <  / —  •»#  4-  f>~ 3*  ^ .  / — m  / r 


V-ft  +  e~i*/4J-xm/c 


(3.57) 

(3.58) 

(3.59) 


=  “‘’fir1)  (3-60) 

C2  =  CSC2  )  (3-61) 

°3  =  csc’p^f~)  (3(i2) 

Ci  =  c,c2(1^")  (3-63) 

The  TD-UTD  slope  diffraction  coefficient  for  a  straight  wedge  obtained  from 

Veruttipong’s  formulation  is  derived  in  Appendix  C.  Here  is  the  same  coefficient 
generalized  for  a  curved  wedge 

A“  (0  =  — E  A™  (*«.,<)  +  Bm  Fs  (xm,t)  (3.64) 

4n* V27T  sin4  Po  mrrl 


A::,  (o = 


■=Fl 


4ti2\/2tc  sin2  j30 


T,A 


m  F vs  (®tnjf)  d”  Bm  Fa  (®rai^) 


m=3 


where  JPS  (a:m,<)  is  defined  in  (3.59), 


and 


with 


0^(13)  —  2  cos2 


(3.65) 

(3.66) 

(3.67) 

(3.68) 

(3.69) 

(3.70) 

(3.71) 

(3.72) 

(3.73) 

(3.74) 
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d^/?)  =  sm(2nirN±  —  /3) 


(3.75) 


The  other  analytic  time  function  is 


+  ,  2  C-y/  2Jm  . 

{Xm,t)  =  - -= —  In 

v71"  L 


e-7’r/4v/-jct  +  \/-*r 


■\/  ®r 


+  jarg  (eW4^ct  + 


(3.76) 


where  arg(z)  is  the  argument  (phase  angle)  of  the  complex  number  z,  y/—xm  is  defined 
in  (3.39)  and  Re  yf—jt  >  0. 

Although  the  slope  diffraction  coefficients  in  (3.57)  and  (3.58)  (or  (3.64)  and 
(3.65))  have  been  generalized  for  a  curved  wedge,  one  must  be  careful  when  the 
illuminated  wedge  face  is  curved.  In  fact,  the  slope  diffracted  field  presented  here 
will  not  fully  compensate  for  the  discontinuous  spatial  derivative  in  the  reflected  field 
when  the  reflecting  surface  is  curved.  Nonetheless,  the  slope  diffracted  field  presented 
here  will  provide  a  good  approximation  for  many  pratical  situations.  For  a  more 
detailed  discussion  on  the  this  problem  one  is  referred  to  the  thesis  by  Zheng  [41]. 


3.3  Some  Numerical  Examples 

First,  this  section  presents  a  convenient  expansion  of  a  general  transient  wave  exci¬ 
tation  which  can  be  convolved  in  closed  form  with  the  TD-UTD  impulse  response  of 
the  curved  wedge,  next  some  numerical  examples  are  included  to  illustrate  the  utility 
of  the  concepts  developed  in  this  chapter.  In  this  section,  Hwang’s  version  of  slope 
diffraction  (instead  of  Veruttipong’s  version)  is  used  unless  stated  otherwise. 

Convolution  of  the  TD-UTD  Impulse  Response  with  a  Gen¬ 
eral  Pulsed  Astigmatic  Wave  Excitation 

It  is  easy  to  examine  the  impulse  response  provided  by  the  TD-UTD  formulas  derived 
in  the  last  section,  but  there  are  two  important  reasons  why  one  should  be  interested 
in  different  types  of  excitations  besides  an  impulsive  wavefront.  The  first  reason  is 
because  the  TD-UTD  is  based  on  the  asymptotic  (high  frequency)  formulas  of  the 
time  harmonic  UTD  and  therefore  the  “late  time”  behavior  of  the  TD-UTD  impulse 
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response  is  expected  to  be  erroneous.  On  the  other  hand,  if  the  TD-UTD  impulse 
response  is  convolved  with  an  excitation  waveform  which  has  a  frequency  spectrum, 
Eq(u;)  as  in  (3.1),  dominated  by  frequencies  for  which  the  time  harmonic  UTD  solution 
is  expected  to  be  accurate  then  the  transient  waveforms  predicted  by  the  TD-UTD 
should  be  very  accurate  also.  The  second  reason  is  because  in  real  physical  situations 
the  frequency  spectrum  of  pulsed  fields  which  can  propagate  without  dispersion  is 
limited.  Also,  there  is  typically  a  lower  limit  on  the  frequency  which  a  source  can 
radiate.  Therefore  the  response  due  to  an  excitation  waveform  which  has  an  effectively 
band  limited  frequency  spectrum,  Eq(u;),  and  of  finite  time  duration  is  more  realistic. 
Furthermore,  one  must  convolve  the  TD-UTD  impulse  response  with  a  more  realistic 
excitation  in  order  to  predict  or  compare  with  measurements  or  other  numerical 
calculations. 

Numerical  convolution  of  any  given  excitation  time  waveform  with  the  impulse 
response  formulas  in  the  previous  section  is  always  possible,  but  it  may  not  be  the  most 
efficient  approach.  A  different  approach  which  is  typically  more  efficient  for  narrow 
pulse  (i.e.  wide  band)  excitations  is  to  express  the  excitation  waveform  as  a  sum  of 
a  small  number  of  simple  expansion  functions  which  can  be  convolved  with  the  TD- 
UTD  impulse  response  in  closed  form.  One  choice  of  a  very  simple  expansion  function 
is  a  waveform  whose  analytic  signal  representation  is  a  simple  pole  in  the  complex 
time  plane.  First,  assume  that  the  frequency  response  of  each  vector  component  in 
the  excitation  Eq(u>)  is  the  same,  so  that 

®iW  =  pHM  (3-77) 


where  p  is  a  polarization  unit  vector,  and  Fq(u)  is  the  Fourier  transform  of  the 
excitation  time  waveform  F^t).  Let 


Fi(t)  =  Re  [f  ‘  (()]  =  Re 


\i  N  A  1 

J_  An 


(3.78) 


with  a  frequency  response  of 

^(u,)  =  f:A„e— ,  w  >  0  (3.79) 

n— 1 
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The  parameters  {A„}  and  {an}  are  in  general  complex.  Now,  suppose  the  analytic 
impulse  response  is  E/™  ( t ),  then  using  the  convolution  property  in  Table  2.1  and 
the  excitation  in  (3.78),  the  response  of  the  curved  wedge  to  a  pulsed  astigmatic  wave 
is 

E UTD(<)  =  \  H  (*)*  eTD  (t)  (3.80) 

=  £^„ErD(i  +  i“n)  (3.81) 

n=l 

and  the  real  response  EUTD(t)  is  the  real  part  of  (3.80)  with  Im  t  =  0.  Notice  that  as 
long  as  Re  an  >  0  for  all  n,  the  analytic  function  resulting  from  the  convolution  in 
(3.80)  is  analytic  on  the  real  time  axis  (Im  t  =  0)  and  therefore  EUTD(i)  is  bounded 
and  well  behaved. 

A  convenient  frequency  domain  window  function,  which  satisfies  the  more  general 


form  in  (3.79),  is 

r?i( t  ,\  /O  / 1  gn—u:T\N„—ujMT 

7*0(0;)  =  G0(l  -  e  )  e 

(3.82) 

where 

T  -,1MutN) 

2i r/c  \  M  J 

(3.83) 

and  fc  is  the  center  frequency.  The  peak  of  the  window  Pq(u>)  is 

normalized  to  one 

by  choosing 

_  (M  +  N\m  /M  +  N\n 

Co=\  M  )  {  N  ) 

(3.84) 

An  Example:  Scattering  from  a  Two-dimensional  Curved 
Strip 

Figure  3.4  shows  the  geometry  of  a  two  dimensional  curved  parabolic  strip.  The  strip 
is  excited  by  an  incident  plane  wave  with  the  electric  field  polarized  perpendicular  to 
the  plane  of  the  paper.  The  observation  point  P,  as  shown  in  Figure  3.4,  is  near  the 
reflection  shadow  boundary  (RSB)  of  edge  (1).  Figure  3.5  shows  the  impulse  response 
calculated  using  the  TD-UTD  when  the  dominant  three  scattering  mechanisms  are 
included.  Figure  3.6  shows  the  impulse  response  corresponding  to  each  scattering 
mechanism  separately,  namely  the  reflected  field,  the  diffracted  field  from  edge  (1) 
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and  the  diffracted  field  from  edge  (2).  Notice  that  the  contribution  to  the  impulse 
response  from  the  edge  diffracted  ray  from  edge  (1)  appears  to  be  very  close  to  -0.5 
times  the  reflected  ray  contribution  to  the  impulse  response.  This  is  to  be  expected 
since  the  observation  point  is  very  near  to  the  RSB  of  edge  (1).  For  verification,  the 
impulse  response  predicted  by  the  TD-UTD  is  convolved  with  a  finite  (band  limited) 
excitation  and  then  compared  with  a  reference  solution.  For  the  frequency  domain 
window  shown  in  the  inset  in  Figure  3.7,  the  parameters  in  the  window  function  of 
(3.82)  are  fc  =  20  GHz,  M  =  1  and  N  =  2. 

The  reference  solution  is  based  on  a  frequency  domain  method  of  moments  (MM) 
analysis  which  has  been  weighted  by  the  appropriate  frequency  window  and  trans¬ 
formed  numerically  into  the  time  domain  using  an  inverse  fast  Fourier  transform 
(FFT).  As  can  be  seen  in  Figure  3.7,  the  reference  solution  MM+FFT  and  the  TD- 
UTD  result  are  almost  indistinguishable  from  each  other. 

Electric  Dipole  Illuminating  a  Wedge 

Figure  3.8  shows  the  space- time  behavior  of  the  total  far  zone  field  as  predicted  by 
the  TD-UTD  when  a  PEC  wedge  is  illuminated  by  an  infinitesimal  electric  current 
element.  The  curves  are  calculated  by  a  convolution  of  the  TD-UTD  impulse  response 
with  a  smooth  finite  energy  time  pulse  with  a  frequency  response  described  by 

F(u>)  =  Ce~au{  1  -  e~au)  for  u;  >  0  (3.85) 

where  C  and  a  are  constants.  The  dipole  is  pointed  just  5  degrees  off  the  direction  to 
the  edge,  so  that  the  incident  field  is  non-zero  and  rapidly  varying  at  the  edge.  This  is 
done  so  that  both  the  first  order  and  slope  diffraction  are  significant.  The  cfiffracted 
field  is  much  smaller  than  the  incident  or  reflected  field  and  is  difficult  to  see  in  Figure 
3.8.  Nonetheless,  notice  that  the  total  field  is  a  smooth  and  continuous  function  of 
both  space  and  time.  This  is  expected  since  the  TD-UTD  diffracted  field  corrects  the 
discontinuities  of  the  geometrical  optics  (GO)  field  at  the  shadow  boundaries. 

Figure  3.9  is  similar  to  Figure  3.8  except  it  shows  the  first  order  diffracted  field. 
The  field  is  also  multiplied  by  a  factor  of  10  to  make  it  visible  on  this  plot  scale. 
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(7)  '  RSB 


Figure  3.4:  Two  dimensional  parabolic  strip  geometry.  The  focus  is  located  at 
F(x=0,z=4),  the  width  of  the  aperture  is  8,  and  the  observation  point  is  P(x=4,z=8). 
All  linear  dimensions  are  in  centimeters. 


Time  (t-p/c)  (ns) 


Figure  3.5:  The  TD-UTD  scattered  field  at  the  observation  point  P  when  the  2-D 
parabolic  strip  in  Figure  5  is  excited  by  an  impulsive  plane  wave.  The  time  axis  is 
shifted  by  p/c  for  convenience  where  p  is  the  distance  of  P  from  the  origin  and  c  is 
the  speed  of  light. 
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Time  (t-p/c)  (ns) 


Figure  3.6:  The  impulse  response  for  each  scattering  mechanism  is  plotted  here  sep¬ 
arately,  where  Er(t ),  Edl(t )  and  Ed2(t )  correspond  to  the  impulse  response  of  the 
reflected  ray,  the  edge-diffracted  ray  from  edge  (1)  and  the  edge-diffracted  ray  from 
edge  (2),  respectively.  The  time  axis  is  shifted  by  p/c  for  convenience. 


Figure  3.7:  Same  as  in  Figure  6  except  that  the  incident  temporal  excitation  is  a 
finite  energy  time  waveform. 
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Notice  that  the  first  order  diffracted  field  is  discontinuous  as  a  function  of  angle  at 
the  shadow  boundaries.  This  discontinuity  corrects  the  discontinuous  GO  field  at 
these  boundaries. 

Figure  3.10  is  similar  to  Figures  3.9  and  3.8,  except  it  shows  the  slope  diffracted 
field.  Notice  that  the  slope  diffracted  field  is  a  smooth  function  of  both  space  (angle) 
and  time,  but  it  does  have  a  discontinuous  derivative  as  a  function  of  angle  at  the 
shadow  boundaries,  which  corrects  the  discontinuous  derivative  of  the  GO  plus  first 
order  diffracted  field. 

Figures  3.11,  3.12  and  3.13  compare  the  TD-UTD  impulse  response  with  an  exact 
solution.  The  exact  solution  is  presented  in  Appendix  B  and  is  based  on  Felsen’s 
work  [22].  The  PEC  wedge  is  illuminated  by  a  infinitesimal  current  element  where 
the  current  has  a  unit  step  function  time  dependence.  This  illuminates  the  wedge  with 
a  spherical  wave,  and  the  first  order  term  in  (1  /r)  has  an  impulsive  time  dependence. 
The  TD-UTD  solution  assumes  that  the  incident  field  is  a  time  impulsive  spherical 
wave  with  a  (1  /r)  distance  dependence.  Figure  3.14  is  a  closeup  view  of  the  plots  in 
Figure  3.13.  The  TD-UTD  based  on  Veruttipong’s  work  is  labeled  “TD-UTD  (Ver)” 
and  the  one  based  on  Hwang’s  work  is  labeled  “TD-UTD  (Hwa)”.  Notice  that  the 
TD-UTD  slope  diffraction  solution  based  on  Veruttipong’s  work  is  more  accurate  for 
early  time,  whereas  the  TD-UTD  slope  diffraction  based  on  Hwang’s  work  is  more 
accurate  for  intermediate  to  late  time.  This  is  anticipated  from  the  equations  since 
the  formulation  based  on  Veruttipong’s  work  contains  an  approximation  for  a  higher 
order  term  in  the  frequency  domain  asymptotic  expansion.  Nonetheless,  it  appears 
from  these  simple  examples  that  the  solution  based  on  Hwang’s  work  is  adequate  for 
many  applications  and  is  still  a  simple  and  more  compact  formulation.  * 

Magnetic  Dipole  Illuminating  a  Wedge 

The  comparison  in  Figures  3.15  and  3.16  is  similar  to  Figures  3.13  and  3.14  except  the 
wedge  is  excited  by  an  infinitesimal  magnetic  current  element.  Figure  3.16  is  simply  a 
close  up  view  of  the  bottom  plot  in  Figure  3.15.  The  TD-UTD  slope  diffraction  based 
on  Veruttipong’s  work  is  labeled  “TD-UTD  (Ver)”  and  the  one  based  on  Hwang’s 
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Electric  Dipole  and  PEC  Wedge,  Total  Field 


Time  (ns) 

Figure  3.8:  A  waterfall  plot  showing  the  space-time  behavior  of  the  total  field  when 
an  electric  dipole  illuminates  a  PEC  wedge.  The  observer  is  in  the  far  zone  and  the 
incident  field  at  the  edge  is  rapidly  varying  so  slope  diffraction  is  significant. 


Time  (ns) 

Figure  3.9:  A  waterfall  plot  showing  the  space-time  behavior  of  the  first  order 
diffracted  field  when  an  electric  dipole  illuminates  a  PEC  wedge  and  the  observer 
is  in  the  far  zone.  The  field  here  is  multiplied  by  a  factor  of  10  as  compared  with 
Figure  3.8. 


Electric  Dipole  and  PEC  Wedge,  Slope  Diffraction  (lOx) 


Time  (ns) 

Figure  3.10:  A  waterfall  plot  showing  the  space-time  behavior  of  the  slope  diffracted 
field  (Hwang’s  version)  when  an  electric  dipole  illuminates  a  PEC  wedge.  The  observer 
is  in  the  far  zone.  The  field  here  is  multiplied  by  a  factor  of  10  as  compared  with 
Figure  3.8. 


58 


Electric  Dipole  and  PEC  Wedge 


2  4  6  8  10  12  14 
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Figure  3.13:  The  two  versions  of  the  TD-UTD  slope  diffraction  are  compared  with  an 
exact  result.  The  infinitesimal  electric  current  has  a  unit  step  time  behavior  which 
approximately  illuminates  the  wedge  with  an  impulsive  spherical  wave. 
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Electric  Dipole  and  PEC  Wedge 


TIME  (ns) 


Figure  3.14:  The  two  versions  of  the  TD-UTD  slope  diffraction  are  compared  with  an 
exact  result.  The  infinitesimal  electric  current  has  a  unit  step  time  behavior  which 
approximately  illuminates  the  wedge  with  an  impulsive  spherical  wave. 
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work  is  labeled  “TD-UTD  (Hwa)”.  The  magnetic  current  element  is  placed  directly 
on  one  of  the  wedge  faces,  since  this  could  be  used  to  model  a  small  slot  antenna. 
Notice  that  there  is  excellent  agreement  for  an  observer  near  the  shadow  boundary 
and  far  from  the  shadow  boundary.  Also,  notice  that  once  again  the  TD-UTD  based 
on  Veruttipong’s  work  is  more  accurate  for  early  time;  whereas,  the  formulation  based 
on  Hwang’s  work  is  more  accurate  for  intermediate  to  late  time  as  expected. 

3.4  Conclusions 

The  TD-UTD  solution  for  the  scattered  field  from  a  general  curved  wedge  is  sum¬ 
marized  as  follows.  The  TD-UTD  impulse  response  for  the  field  in  the  presense  of  a 
general  curved  wedge  may  be  written  as 

E/TD  (0  =E/  (t)Ui+  E;  (t)Ur+  E/  (0+  E/d  (f)  (3.86) 

where  the  spatial  unit  step  functions  f/,  and  Ur  are  1  on  the  lit  side  of  the  incident 
shadow  boundary  and  the  reflection  shadow  boundary,  respectively,  and  0  otherwise. 
The  TD-UTD  incident  GO  field  E /  ( t )  is  in  (3.10)  and  the  TD-UTD  reflected  field 
E/  (t)  is  in  (3.14).  The  first  order  TD-UTD  edge  diffracted  field  E/  ( t )  is  in  (3.16) 
where  the  TD-UTD  dyadic  diffraction  coefficient  is  defined  in  (3.40)  and  (3.41).  The 

-f  grf 

higher  order  edge  diffracted  field  called  the  slope  diffracted  field  E j  (t)  is  in  (3.56) 
where  the  TD-UTD  dyadic  slope  diffraction  coefficient  based  on  Hwang’s  frequency 
domain  UTD  slope  result  is  given  in  (3.57)  and  (3.58)  while  the  slope  diffraction 
coefficient  based  on  Veruttipong’s  frequency  domain  UTD  result  is  given  in  (3.64) 
and  (3.65).  Finally,  the  impulse  response  formulas,  denoted  by  the  subscript  7,  may 
be  convolved  with  an  excitation  in  closed  form  as  in  (3.80)  and  (3.81)  when  the 
excitation  has  the  form  in  (3.77)  or  (3.78). 

This  chapter  presents  a  time  domain  version  of  the  UTD  (TD-UTD)  for  the  general 
case  of  a  curved  PEC  wedge,  where  the  wedge  faces  and  the  edge  may  be  curved.  This 
time  domain  solution  can  handle  more  general  geometries  than  what  has  been  done 
in  the  past  and  supplies  a  good  stepping  stone  for  the  development  of  a  complete 
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Magnetic  Dipole  on  PEC  Wedge 
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Figure  3.15:  The  two  versions  of  the  TD-UTD  slope  diffraction  are  compared  with  an 
exact  result.  The  infinitesimal  magnetic  current  has  a  unit  step  time  behavior  which 
approximately  illuminates  the  wedge  with  an  impulsive  spherical  wave. 
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Figure  3.16:  The  two  versions  of  the  TD-UTD  slope  diffraction  are  compared  with  an 
exact  result.  The  infinitesimal  magnetic  current  has  a  unit  step  time  behavior  which 
approximately  illuminates  the  wedge  with  an  impulsive  spherical  wave. 


TD-UTD  which  includes  multiple  diffractions,  caustic  corrections,  surface  diffraction, 
etc. 
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Chapter  4 

TD-UTD  for  a  Smooth  Convex 
Surface 


The  study  of  electromagnetic  surface  diffraction  has  been  an  active  area  of  research  for 
many  years  both  in  the  frequency  domain  (i.e.  the  time  harmonic  case)  and  also  in  the 
time  domain.  Smooth  surface  diffraction  as  described  here  assumes  that  the  surface 
is  a  perfect  electric  conductor.  Even  with  this  assumption,  smooth  surface  diffraction 
is  a  very  important  phenomena  since  many  geometries  encountered  in  engineering 
contain  conducting  smooth  surfaces,  such  as  airplanes,  automobiles,  etc.  The  work 
in  this  report  is  based  on  the  frequency  domain  UTD  solution  and  Pathak  provides 
an  excellent  summary  of  this  work  in  [4]. 

There  is  an  abundance  of  past  research  on  the  topic  of  diffraction  from  smooth 
surfaces  which  are  excited  by  short  pulses,  either  acoustic  pulses  or  electromagnetic 
pulses.  One  of  the  classic  works  on  acoustic  pulse  scattering  is  a  book  by  Friedlander 
[7].  In  chapter  6  of  [7],  Friedlander  formulates  the  time  domain  acoustic  Green’s  func¬ 
tion  (i.e.  impulse  response)  for  a  circular  cylinder  and  also  for  a  sphere.  He  obtains 
an  approximate  “early  time”  formula  for  the  diffraction  from  a  circular  ^cylinder  by 
implementing  an  asymptotic  evaluation  of  an  inverse  Laplace  transform  of  the  fre¬ 
quency  domain  result.  This  early  time  result  is  a  summation  of  time  domain  creeping 
wave  modes  and  therefore  is  not  valid  when  the  observer  is  near  the  surface  shadow 
boundary  (SSB),  but  it  is  still  very  useful  for  “deep  shadow”  observer  locations  or 
backscatter  calculations.  Friedlander  solves  the  acoustic  scattering  from  a  sphere  sim¬ 
ilarly.  Also,  he  finds  an  explicit  early  time  formula  for  when  the  observer  is  located  at 
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the  caustic  in  the  shadow  region.  But,  Friedlander  does  not  obtain  the  solution  to  the 
sphere  problem  for  the  case  when  the  diffracted  ray  field  has  passed  through  a  caus¬ 
tic.  In  [42],  Weston  obtains  the  back-scattered  electromagnetic  field  from  a  perfectly 
conducting  sphere  when  the  incident  field  is  a  plane  wave  and  the  temporal  excitation 
is  a  modulated  square  pulse.  Weston  finds  a  power  series  approximation  for  small 
(or  early)  time  from  the  inverse  Laplace  transform  of  the  Luneberg-Kline  expansion. 
For  large  (or  late)  time,  Weston  uses  an  eigenfunction  representation  and  when  the 
sphere  radius  is  large  in  terms  of  the  wavelength  of  the  modulation  frequency  he  used 
a  creeping  wave  representation.  Wait  and  Conda  discussed  diffraction  of  electromag¬ 
netic  pulses  by  curved  surfaces  in  [43].  The  first  section  of  Wait  and  Conda’s  paper 
is  devoted  to  finding  the  electric  currents  induced  on  a  circular  cylinder  when  it  is 
excited  by  a  plane  wave.  The  temporal  excitation  is  a  step  function.  The  currents  in 
the  frequency  domain  are  expressed  as  an  approximate  integral  expression  containing 
an  Airy  integral.  This  approximate  frequency  domain  result  is  valid  when  the  radius 
of  the  cylinder  is  large  in  terms  of  wavelength.  They  then  expand  this  frequency 
domain  expression  into  a  power  series  and  use  an  inverse  Laplace  transform  to  obtain 
an  asymptotic  time  series  in  terms  of  inverse  powers  of  time.  Although  it  appears 
that  they  use  both  high  frequency  and  low  frequency  approximations  here,  one  can 
interpret  their  solution  to  be  valid  for  intermediate  time  values  and  observations  near 
the  shadow  boundary.  In  the  second  half  of  the  paper,  Wait  and  Conda  [43]  dis¬ 
cuss  the  diffraction  from  a  smooth  surface  when  the  source  and  observer  are  removed 
from  the  surface.  They  solve  the  specific  case  of  a  circular  cylinder,  although  they 
mention  at  the  end  of  the  paper  that  these  results  could  be  easily  generalized.  This 
solution  is  valid  for  observations  near  the  shadow  boundary.  In  [44],  Chen  dbmputes 
the  large-time  transient  behavior  of  the  diffracted  field  from  a  circular  cylinder  ex¬ 
cited  by  a  plane  wave  with  a  step  function  temporal  behavior.  He  solves  the  problem 
for  various  types  of  boundary  conditions  at  the  surface  and  also  pays  close  attention 
to  causality.  Chen  bases  his  solution  on  the  exact  frequency  domain  eigenfunction 
summation.  He  uses  an  inverse  Laplace  transform  on  the  eigenfunction  summation, 
applies  the  transform  term  by  term,  and  then  evaluates  each  term  asymptotically  for 
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large  time.  Uberall,  Doolittle  and  McNicholas  solve  the  acoustic  scattering  problem 
of  a  cylinder  excited  by  an  impulsive  plane  wave  (i.e.  the  impulse  response)  [45]. 
They  first  solve  the  frequency  domain  problem  and  then  obtain  the  impulse  response 
by  using  an  inverse  Laplace  transform.  Since  their  frequency  domain  solution  is  an 
asymptotic  high  frequency  result,  their  impulse  response  is  valid  for  early  (or  small) 
times.  Also,  in  their  solution  the  scattered  field  is  decomposed  into  a  reflected  field 
and  diffracted  field  which  is  written  as  a  sum  of  creeping  wave  modes.  Therefore, 
their  solution  is  not  valid  near  the  shadow  boundary.  In  [46]  and  Appendix  II  of  [1], 
Moffatt  explores  that  nature  of  the  creeping  wave.  He  does  this  by  examining  the 
early  time  backscattered  electromagnetic  field  from  a  circular  cylinder  and  a  sphere 
which  he  obtains  by  using  an  inverse  Laplace  transform  on  the  asymptotic  high  fre¬ 
quency  results.  He  discovers  that  the  early  time  backscattered  diffracted  field  from 
the  sphere  is  the  Hilbert  transform  of  a  causal  time  function  and  then  correctly  pre¬ 
dicts  that  the  early  time  approximate  backscattered  diffracted  field  from  a  sphere 
is  not  causal.  This  is  easily  understood  now,  since  we  know  that  the  arrival  of  the 
reflected  field  represents  the  “turn  on”  time  of  the  total  scattered  field,  and  therefore 
the  “small  time”  solution  for  the  diffracted  field  contribution  may  contain  information 
before  and  after  the  time  of  arrival  of  the  diffracted  ray.  This  is  consistent  with  the 
asymptotic  nature  of  the  “small  time”  solution.  Wait  investigates  the  currents  on  a 
cylinder  when  the  cylinder  is  excited  by  an  impulsive  plane  wave  in  [47].  Wait  first 
formulates  the  solution  in  the  frequency  domain  and  approximates  the  currents  in  the 
penumbral  (near  the  shadow  boundary)  region  as  a  contour  integral  containing  Airy 
functions.  He  expands  this  representation  into  a  creeping  wave  mode  series  and  uses 
an  inverse  Laplace  transform  term  by  term  to  obtain  a  time  domain  creeping  wave 
mode  series.  He  is  able  to  transform  each  creeping  wave  mode  (for  the  currents  in  the 
shadow  region)  in  closed  form  in  terms  of  Airy  functions.  In  [48,  49],  Schafer  calcu¬ 
lates  the  currents  on  a  circular  cylinder  when  the  cylinder  is  excited  by  an  impulsive 
plane  wave.  He  finds  an  early  (small)  time  representation  in  the  lit  region  by  using  an 
inverse  Laplace  transform  on  the  frequency  domain  Luneberg-Kline  series  which  gives 
a  power  series  time  representation.  He  obtains  an  early  (small)  time  representation 
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in  the  shadow  region  by  applying  the  inverse  Laplace  transform  to  the  high  frequency 
creeping  wave  summation.  He  obtains  a  large  time  representation  by  applying  the 
inverse  Laplace  transform  to  the  frequency  domain  eigenfunction  summation  term 
by  term.  Finally,  for  intermediate  time  values  he  uses  an  interpolation  scheme  to 
combine  the  small  and  large  time  representations.  Lee,  Jamnejad  and  Mittra  exam¬ 
ine  the  early  time  response  for  the  scattered  field  and  current  on  a  circular  cylinder 
by  using  an  inverse  Laplace  transform  on  the  high  frequency  creeping  wave  repre¬ 
sentation  [33].  In  [50],  Heyman  and  Felsen  develop  a  solution  for  the  currents  on  a 
circular  cylinder  (or  a  smooth  surface  in  general)  when  the  excitation  is  located  on  the 
cylinder.  Their  solution  combines  a  low  frequency  eigenfunction  representation  with 
a  high  frequency  creeping  wave  representation,  then  they  do  a  similar  combination  in 
the  time  domain.  Ma  and  Ciric  in  [51],  obtain  an  early  time  solution  for  the  currents 
on  a  circular  cylinder  when  it  is  excited  by  a  line  current.  Their  solution  is  similar  to 
Schafer’s  for  the  plane  wave  case.  In  [52],  Naishadham  and  Yao  present  the  transient 
scattering  by  a  circular  cylinder  when  it  is  excited  by  magnetic  line  source  when  the 
source  and  observation  point  are  on  the  cylinder.  They  obtain  their  solution  using 
an  efficient  numerical  evaluation  of  a  contour  integral  which  approximates  the  exact 
contour  integral  solution. 

The  time  domain  curved  surface  diffraction  presented  in  this  chapter  is  different 
from  the  solutions  discussed  above  since  the  present  work  is  based  on  the  frequency 
domain  UTD  [4],  The  frequency  domain  UTD  formulation  of  surface  diffraction  can 
be  applied  to  a  very  general  class  of  smooth  surfaces,  whereas  most  of  the  previous 
time  domain  solutions  discussed  above  involve  canonical  geometries  such  as  cylinders 
or  spheres.  More  importantly,  the  previous  time  domain  solutions  discussed  above 
assume  simple  excitations  such  as  a  plane  wave  or  line  source,  whereas  in  the  TD- 
UTD  solution  presented  here,  the  excitation  is  an  astigmatic  ray  field  where  plane, 
cylindrical  and  spherical  wavefronts  are  special  cases.  Also,  with  the  analytic  time 
representation  used  here,  one  could  easily  incorporate  incident  ray  fields  which  have 
passed  through  a  smooth  caustic,  incident  ray  fields  that  are  circularly  polarized  and 
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diffracted  ray  fields  which  have  passed  through  a  smooth  caustic.  This  greatly  in¬ 
creases  the  usefulness  of  the  present  TD-UTD  solution  for  smooth  surface  diffraction. 

The  content  of  this  chapter  is  as  follows.  First  there  is  a  brief  review  of  the 
frequency  domain  UTD  formulation  for  smooth  surface  diffraction.  Then  the  TD- 
UTD  is  obtained  for  observers  in  the  deep  lit  and  deep  shadow  regions  by  applying 
the  analytic  time  transform  to  the  frequency  domain  UTD  formulas.  Also,  a  TD- 
UTD  solution  is  obtained  for  observers  near  the  shadow  boundary.  Finally,  all  of 
these  solutions  are  combined  to  obtain  a  uniform  TD-UTD  surface  diffraction  solution 
which  remains  valid  in  the  deep  shadow  region,  near  the  shadow  boundary  and  deep 
in  the  lit  region.  This  uniform  TD-UTD  solution  is  the  time  domain  version  of  the 
UTD  solution  in  Pathak  [4].  This  solution  is  restricted  to  the  “scattering  case”  where 
the  observer  and  source  are  located  off  the  smooth  surface,  although  the  radiation 
and  coupling  problems  could  be  dealt  with  in  a  similar  fashion.  Finally,  some  simple 
numerical  results  are  presented  to  demonstrate  the  TD-UTD  and  compare  it  with  a 
reference  solution. 

4.1  General  Formulation  in  the  Frequency  Do¬ 
main 

This  chapter  discusses  only  the  smooth  surface  scattering  situation  where  the  source 
and  observer  are  off  of  the  surface.  Although  the  radiation  and  coupling  situations 
where  the  source  or /and  the  observer  are  located  on  the  smooth  surface  are  of  con¬ 
siderable  interest  only  the  scattering  situation  is  considered  here.  The  radiation  and 
coupling  situations  could  be  analyzed  in  a  similar  fashion  as  the  scattering  problem 
studied  here. 

The  frequency  domain  UTD  solution  for  the  smooth  surface  diffraction  geometry 
shown  in  Figure  4.1  is  as  follows.  The  UTD  formulation  for  smooth  surface  diffraction 
briefly  presented  here  is  based  on  the  summary  by  Pathak  [4] .  For  the  scattering  case 
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Figure  4.1:  Geometry  of  the  smooth  surface  scattering  situation  where  the  source  and 
observer  are  off  the  surface. 


shown  in  Figure  4.1,  the  total  field  as  approximated  by  the  UTD  may  be  written  as 

E^P)  =  |  E'tft)  +  E ”(Pl)  P  =  Pi  in  lit  region  ^ 

E d(Ps)  P  =  Ps  in  shadow 

The  surface  diffracted  field  Erf(Ps)  may  also  appear  in  the  lit  region  when  the  surface 
is  closed  and  one  includes  multiple  encirclements  around  the  surface.  The  incident 
field  E^P^)  is  the  usual  GO  incident  field  which  exists  in  the  lit  region.  The  reflected 
field  E 9T{Pl)  is  a  generalized  reflected  field  which  behaves  as  the  usual  GO  reflected 
field  deep  in  the  lit  region,  but  behaves  quite  differently  near  the  surface  shadow 
boundary  (SSB),  it  can  be  expressed  in  terms  of  generalized  reflection  coefficients 


Wr{PL)  ~  E*(Q*)  •  [R.ae±e±  +  the jejj]  Ar(sr)e~^r  (4.2) 

where  E*(<5/i)  is  the  incident  field  at  the  point  of  reflection  Qr.  The  unit  vectors  in 
(4.2)  are  the  same  as  in  the  usual  GO  reflection  formula  discussed  in  the  previous 
chapter.  The  spreading  factor  for  the  reflected  field  is  given  by 

A^S  )  =  /(*>!+ Sr)(/>2  +^)  (4‘3) 
where  p\  and  pr2  are  the  reflected  ray  caustic  distances  and  sr  is  the  distance  from  the 
reflection  point  Qr  to  the  observer  at  P.  The  surface  diffracted  field  can  be  written 
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Figure  4.2:  Ray  vectors  for  surface  diffraction. 


as 

E d{Ps)  ~  E'(Qi)  •  [V.bibi  +  %AiA2]  (4.4) 

where  E’(Qi)  is  the  incident  field  at  the  attachment  point  as  shown  in  Figure  4.1. 
The  unit  ray  vectors  in  (4.4)  are  shown  in  Figure  4.1.  The  unit  vectors  ni  and  n2  are 
normal  to  the  surface  at  points  Q\  and  Q2  respectively.  The  unit  tangent  vectors  ti 
and  t2  point  in  the  direction  in  which  the  “surface  ray”  is  moving  at  points  Q\  and 

A  A  A 

Q2  respectively.  The  other  tangent  vectors  bli2  are  defined  by  bii2  =  ti)2  x  ni>2.  The 
spreading  factor  for  the  diffracted  ray  is 


M>d) 


Sd(p,  +  3d) 


■it' 

(4.5) 


where  pa  is  the  caustic  distance  of  the  diffracted  ray  and  sd  is  the  distance  from 
the  point  of  diffraction  to  the  observer.  The  frequency  domain  generalized  reflection 
coefficient  can  be  written  as 


fcs,h  = 


IHr-MZ1)*/™ /. 


-4 


z:L 


\  2EL^/iry/juj 


- L M  '  —  O 


[1  -  Ftr(  u>xL)}  +  —P3,h(  u>-< 


o;e 


(4.6) 
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where 

M(-)  =  j  ‘  (4.7) 

Eh  =  -2M{QR)cos9i  (4.8) 

xL  =  2—  cos2  6'  (4.9) 

c 

and  6'  is  the  angle  between  the  surface  normal  vector  and  incident  ray  direction  at 
the  point  of  reflection.  The  function  Ftr(z)  is  the  UTD  transition  function  given  in 
(3.33)  and  the  Fock  type  function  Pa,h{z )  is  described  shortly.  The  P9(Qr)  in  the 
above  equations  is  the  radius  of  curvature  of  the  surface  at  the  point  of  reflection  Qr 
in  the  plane  of  incidence.  The  frequency  domain  UTD  surface  diffraction  coefficient 
is  given  by 

V,„  =  -,/2cM(g1)M(<?2)  -  £.(“'*■')] 


where 


(4.10) 

jQl  Pg(l  ) 

(4.11) 

x  2  cM(Q1)M(Q2) 

(4.12) 

tQi 

1=  dl’ 

(4.13) 

The  path  from  Q\  to  Q2  on  the  surface  is  a  geodesic  where  the  path  length  is  a 
minimum.  The  pg(V )  in  (4.11)  is  the  radius  of  curvature  of  the  surface  at  V  on  the 
geodesic  path  in  the  1'  (or  t)  direction.  The  distance  parameter  in  (4.9)  a^d  (4.12) 
when  evaluated  at  the  shadow  boundary  is  given  by 

J  __  _ Pl(Ql)p2(Ql) _ S(Pb(Q1)  f)  (A 

{p\{Qi)  +  s){p2(Qi)  + s)  Pb(Qi) 

where  p\  2{Q\)  are  the  principle  radii  of  curvature  of  the  incident  ray  field  at  Qi  and 
pl(Qi)  is  the  radius  of  curvature  of  the  incident  ray  field  in  the  bi  direction  at  Q\. 
Note  that  the  distance  parameter  L  in  (4.9)  and  (4.12)  would  in  general  be  different 
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for  the  lit  and  shadow  regions  but  it  is  usually  more  convenient  to  approximate  L  by 
its  value  at  the  shadow  boundary  given  by  (4.14)  [3]. 

For  the  special  case  of  an  incident  spherical  wave 

L  =  (4-15) 


*' = = Am = fim 


(4.16) 


where  s'  is  the  distance  from  the  source  to  the  point  Q\  on  the  surface. 

The  Fock  type  functions  P„AX)  are  related  to  the  soft  and  hard  Pekeris  functions 
p*(x)  and  q*( x)  by 


p  /  \  P  W  I  -w 4  e~W4 

P*AX)  =  S  \e  3  ,A-  —j=- 

q*(x)  I  2V*X 


(4.17) 


where  p*(x)  and  q*(x)  are  bounded  and  well  behaved.  The  Fock  type  function  PSi/,(*) 


can  also  be  written  as 


e-w  r° o  QV(z)  ■ 

p-h(x)  =  -jrLam7)e  dz 


where  the  operator  Q  is 


The  Fock  type  Airy  functions  are 


Vtt  J-°oQW2{z ) 


1  soft  case 

d/dz  hard  case 


(4.18) 


(4.19) 


2 jV(z)  =  W,(z)  -  W2(z) 


1  /■  ■ 


(  7T  Joo  exp(— j27r/3) 

1  r°° 


'y3/3dy 


W2(z)  =  -r  ^  e^-"3/3 dy 

y/TZ  Joo  exp(+j27r/3) 


(4.20) 

(4.21) 


(4.22) 


4.2  Deep  Lit  Region 

The  frequency  domain  UTD  formulation  for  the  fields  in  the  lit  region  can  be  simplified 
when  the  observer  is  far  from  the  shadow  boundary.  When  the  observer  is  in  the  deep 
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lit  region  (i.e.  far  from  the  shadow  boundary),  xL  is  large  and  —  ZL  is  large  (or  ZL  «C  0) 
and  the  following  approximations  can  be  made 

Ftr(  WXL  )  ~  1  for  xL  —>  oo  (4.23) 

Pa,h(u;*ZL)~±w*J^^ejulEL)3/12  for  EL  — >  — oo  (4.24) 

Using  the  above  approximations  for  the  deep  lit  region,  the  generalized  reflection 
coefficient  becomes 

Ha,h  ~  1  deep  lit  region  (4.25) 

where  the  ~  above  the  reflection  coefficient  77  has  been  dropped  since  it  is  no  longer  a 
function  of  frequency  with  this  approximation.  It  is  now  obvious  that  the  generalized 
reflection  coefficient  'R.a<h  simply  becomes  the  ordinary  GO  reflection  coefficient  Ra,h 
discussed  in  Chapter  3  when  the  observer  is  in  the  deep  lit  region.  Now  the  TD-UTD 
field  in  the  deep  lit  region  can  be  written  as 

E9  (-Pl;0~E  (Qr]  t  —  5VC)  *  +  ”^/»e||eJj]ylr(5r)  (4.26) 

where  7 Zs,h  is  given  in  (4.25).  Notice  that  this  is  an  early  time  deep  lit  region  result 
which  directly  corresponds  to  the  high  frequency  deep  lit  region  approximations.  The 
early  time  deep  lit  region  result  in  (4.26)  is  not  expected  to  be  accurate  for  late 
times,  since  the  deep  lit  region  approximations  made  in  the  frequency  domain  are 
increasingly  accurate  as  the  frequency  increases.  Nonetheless,  the  result  in  (4.26)  is 
a  good  approximation  for  early  time  and  for  an  observer  deep  in  the  lit  region. 


4.3  Deep  Shadow  Region 


The  frequency  domain  UTD  formulation  can  be  simplified  when  the  observer  is  deep 
in  the  shadow  region.  When  the  observer  is  in  the  deep  shadow  region,  xd  in  (4.12) 
and  2  in  (4.11)  are  large  and  positive  and  the  following  approximations  can  be  made 


Ftr{ wxd)  ~  1  for  xd  — >  oo 


A,h(w®  2)  ~  < 


-e  J7r/4  A  eJs  exp(  — (yu;)3  2qn 
v/ir  ^  2[Ai,(-9n)]2 

_e-W4  ^  ej f  exp(— (yu>)3  2qn 


“i  2g„  [Ai  (-gn)]: 


(4.27) 

(4.28) 
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for  S  — >  <X).  In  (4.28),  Ai(x)  =  V(x)ItJic  is  a  Miller  type  Airy  function  and  Ai^x)  = 
dA.\(x)/dx  is  its  derivative  [29].  The  roots  of  Ai(— qn)  =  0  and  Ai'(  — qn)  =  0, 
denoted  by  qn  and  qn  respectively,  are  also  used  in  (4.28).  Typically  two  terms  (N  =  2) 
are  used  for  the  summation  in  (4.28)  [4].  The  surface  diffraction  coefficient  in  (4.10) 
can  now  be  simplified  for  an  observer  in  the  deep  shadow  region  by 
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s,h 


1  N 

1  V' 

2c 

-M(Q1)M{Qa\ 
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0«)*  n=l 
1  N 
ttVe 

UW)6  n=l 

ex-p(-(ju)*Zqn) 


exp(-(ju)*Zqn) 


(4.29) 


Finally,  the  analytic  time  transform  can  be  applied  to  (4.29)  to  obtain  the  analytic 
time  representations  of  the  surface  diffraction  coefficient 


Vs,h  (: t ) 


where  T-p  =  t  —  l/c  and 


2c 


7T 


M(Q1)M(Q2) 


MQi) 

dT](Q2) 


Few  (“(7n}TZ>) 

h  2[Ai'(-fc)]» 
l  h\  2?n[Ai {-qn)}2 


(4.30) 


Few  (a,t)  =  -  f  ( ju) )  « e  “0“)’  eJU)tdu> 
7 r  Jo 


(4.31) 


for  Im  t  >  0.  The  analytic  time  function  in  (4.31)  is  a  new  special  function  which 
must  be  calculated.  Appendix  D  explains  in  detail  how  to  calculate  Few  ( a ,£)  in 
(4.31). 


4.4  Observer  Near  the  Shadow  Boundary 


This  section  presents  the  TD-UTD  smooth  surface  reflection  and  diffraction  coeffi¬ 
cients  for  when  the  observer  is  near  the  shadow  boundary.  First  rewrite  the  frequency 
domain  generalized  reflection  coefficient  in  (4.6)  as  the  following 


JzlP-MzL)3n2  < 

r 

Ftr(  loxl)  e-j*'6 

f(w'/3SL)  ' 

L 

V  -Le  1 

V  t-j 

2SLv/7rv/yo;  (yw)1/6 

k 

r(w1/32L)  J 

(4.32) 


and  the  surface  diffraction  coefficient  of  (4.10) 

Vs,h  =  -yji cM(Q1)M(Q2) 


as 

Ftrj  t*KCd) 
2Zy/lCy/JuJ 
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g-iTT/e 

GH1/6 


p*(  u1'3-) 
q*(^1/3E) 


MQi)  -jui/c 

^  dr)(Q2) 


(4.33) 


Now  apply  the  ATT  to  (4.32)  and  (4.33)  to  obtain 


where 


and 


+  J^iF{cx\Tn)_+  (Z}L 

Us'h  VHL1  2^Tc=L  *p’q[~ 


TK  =  t-  (El)3/12 
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Vs,h=  ft cM{Qx)M{Qt)  F^£)~  ita  (S,ty>) 
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jWi) 

^(Q2) 


where 


Tp  =  t  —  l/c 
+ 


(4.34) 


(4.35) 


(4.36) 


(4.37) 


See  Appendix  E  for  a  series  expansion  for  Fp,q  (2,  t)  which  is  appropriate  when  |E| 
is  small  (i.e.  when  the  observer  is  close  to  the  shadow  boundary).  The  other  time 
function  used  in  (4.34)  and  (4.36)  is  related  to  the  TD-UTD  wedge  transition  function 
in  (3.37)  by 


Ftr(uix) 


[  VJw  J 


1  + , 

=  ScF{CX,t) 


■Jxfi 


+  VJx) 

At  the  shadow  boundary,  xd  — >  0  and  xL  —*  0,  therefore 

-j-  _  -j- 

F  ( cx,t )  ~  \/itcx  8  ( t ) 

Also,  EL  — >  0“  and  E  — »  0+  so  that  the  following  could  be  used 


(4.38) 

(4.39) 


(4.40) 


£  (o  t)  -  r(5/6)  e—  < 
}t  (-jf)s/6 


Po 

9o 


(4.41) 


where  po  —  0.354064  and  qo  ==  —0.307177.  Equations  (4.40)  and  (4.41)  could  be  used 
to  demonstrate  that  the  TD-UTD  is  continuous  across  the  surface  shadow  boundary, 
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just  as  is  done  in  [3]  for  the  frequency  domain  case.  The  term  containing  8  (4)  function 
will  produce  minus  one  half  of  the  incident  field  when  used  in  the  reflection  coefficient 
and  plus  one  half  of  the  incident  field  when  used  in  the  diffraction  coefficient.  But, 
the  term  containing  (— jt)~ 5/6  produces  a  diffracted  field  effect  even  on  the  shadow 
boundary.  This  is  in  contrast  to  the  first  order  diffraction  from  the  wedge  which 
becomes  one  half  of  the  incident  field  at  the  shadow  boundary. 

4.5  Uniform  TD-UTD  Surface  Diffraction 


This  section  presents  the  uniform  TD-UTD  formulation  for  the  scattering  from  a 
smooth  convex  surface.  This  requires  the  analytic  time  transform  (ATT)  of  the 
corresponding  frequency  domain  UTD  generalized  reflection  coefficient  in  (4.6)  and 
surface  diffraction  coefficient  in  (4.10).  The  following  transform  pairs  are  used 
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g-ix/4 


(4.42) 


and 


1  - 


Ps,k(u1/3Z)\  =Fs<h  (2,4) 


U J 


1/6 


(4.43) 


along  with  (4.38).  The  transform  in  (4.43)  can  be  used  when  E  <  0  and  when  E  >  0, 
so  it  can  be  used  for  the  generalized  reflection  coefficient  and  the  surface  diffraction 
coefficient.  The  special  time  function  Fa%h  (2,4)  is  not  known  in  closed  form,  but 
instead  it  can  be  calculated  using  various  expansions  for  early  time,  late  time  and  so 
forth.  Appendix  E  fully  explains  the  properties  of  the  Fah  (2,  t)  function  and  explains 
the  algorithm  used  for  its  calculation.  Nonetheless,  some  basic  results  concerning  the 
F^h  (2,4)  are  now  summarized.  .. 

In  the  shadow  region,  S  >  0  and  the  following  expansions  can  be  used.  When 
the  observer  is  in  the  deep  shadow  region  (or  for  very  early  time)  a  creeping  wave 
expansion  can  be  used 


i^(2,4)~  < 


1  Few  (,=,9n5^) 

2|,Ai'(-«„)]2 

1  ^  Few  (l— '9n,4) 


(4.44) 
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for  \t/E3\  — >  0,  where  Few  (<*,£)  are  the  time  domain  creeping  wave  mode  functions 
described  in  Appendix  D.  When  the  observer  is  near  the  shadow  boundary  (or  for 
late  time)  an  inverse  time  power  series  can  be  used 


where 


+p  ,  1  e"-^4 

Fs,h  (->*)  -Fp, g  (=,«)  -  2lrE  ^yt 


+  N  e-j7r/12  ^  J  pn  r(n/3  +  5/6)e^/6  En 


(4.45) 


(4.46) 


for  |t/E3|  +  0. 

For  the  lit  region,  EL  <  0,  and  (4.46)  can  be  used  when  \t/(EL)3\  — ►  oo  by  simply 
replacing  E  in  (4.46)  with  EL.  But,  for  \t/(EL)3\  ->  0  the  creeping  wave  mode  series 
in  (4.44)  is  not  valid  in  the  lit  region,  instead  one  can  use 


Fe,h  (s>  0  ~  ± Y  6  (*  +  (2L)3/12) 


(4.47) 


for  [£  +  (EL)3/12]/(EL)3  0.  See  Appendix  E  for  more  details  about  the  computa- 

-|-  p 

tion  of  these  expansions  and  Fa<h  (—>0  in  general. 

Now,  using  the  transform  relationships  in  (4.38),  (4.42)  and  (4.43),  transform  (4.6) 


+  /— 4  e-J*74  +p 

%s,h=  -\/=r  { - 7 - +  F**  ( 3  ’r*) 

V=L  2s'-x(v^ji+7^) 


(4.48) 


and  transform  (4.10)  into 


Va,H=  -\]2cM(Qi)M(Q2)  - 7 - t=v+  Fa,h  (H,**)  *  (4.49) 

2Ett  J 

It  is  interesting  to  notice  how 

e-jw/4 

\/~jt  +  \/Jx 

behaves  for  x  >  0  and  real  time  Im  (£)  =  0.  The  above  can  be  rewritten  as 

\Zt~  jy/x 

t  +  * 
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and  notice  that  this  gives 


Re 


) 


y/tu(t) 
t  +  x 


for  x  >  0  and  Im  (t)  =  0.  This  is  the  same  type  of  transient  behavior  observed  for 
the  slope  diffraction  transition  function  (C.17). 


4.6  Numerical  Example 

Figures  4.4  through  4.8  show  a  comparison  between  a  TD-UTD  solution  and  an 
eigenfunction  reference  solution  for  the  TE*  (soft)  far  zone  scattering  from  a  two 
dimensional  (2-D)  cylinder  at  various  bistatic  angles.  The  reference  eigenfunction 
solution  is  a  frequency  domain  solution  [53]  which  is  multiplied  by  the  frequency 
domain  window  in  Figure  4.3  and  transformed  into  the  time  domain  using  the  IFFT. 
The  time  pulse  excitation,  in  the  top  plot  of  Figure  4.3,  corresponding  to  the  frequency 
window  in  the  bottom  plot  of  Figure  4.3,  is  convolved  with  the  TD-UTD  analytic 
impulse  response  in  closed  form  where  only  two  frequency  domain  exponentials  are 
used  (see  (2.67)  in  Chapter  2).  In  the  far  zone,  xd  — »  oo  and  xL  — >  oo  so  that  the 
coefficients  ils,h  and  T>s<h  are  proportional  to  the  j p  h  (E,£)  function.  The  source  and 
observer  locations  are  chosen  to  be  in  the  far  zone  for  this  comparison  because  the 

-| ~  p 

Fa>h  (H,  t)  function  is  more  difficult  to  calculate  than  the  other  terms  in  reflection 
coefficient  j Zs,h  and  diffraction  coefficient  T>a,h  and  it  is  therefore  desirable  to  check 
the  accuracy  of  this  Fs<h  (H,t)  function.  The  radius  of  the  cylinder  is  r  =  1  meter 
which  makes  the  radius  one  wavelength  at  the  peak  frequency  of  the  window  function 
in  Figure  4.3.  In  Figures  4.4  through  4.13,  the  top  plot  shows  the  magnitude  of  the 
analytic  time  function  while  the  bottom  plot  shows  the  real  part  of  the  analytic  time 
function.  Although  the  bottom  plots  (i.e.  the  real  part)  are  the  actual  time  signals 
arising  from  the  excitation  in  Figure  4.3,  the  top  plots  (i.e.  the  magnitude)  may  be 
thought  of  qualitatively  as  the  signal  from  an  envelope  detector.  Notice  that  there  is 
excellent  agreement  between  the  TD-UTD  solution  and  the  eigenfunction  reference 
solution,  which  confirms  that  the  numerical  algorithm  for  computing  Fah  (H5t)  i*1 
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Frequency  (GHz) 

Figure  4.3:  Excitation  pulse  w(t)  and  corresponding  frequency  response  (or  frequency 
window)  W (/)  used  for  the  scattering  from  a  2-D  cylinder  example. 

Appendix  E  is  working  well  for  this  case  and  that  the  TD-UTD  is  very  accurate  for 
this  case. 

Figures  4.9  through  4.13  are  essentially  the  same  as  the  Figures  4.4  through  4.8  just 
discussed  except  that  Figures  4.9  through  4.13  are  for  the  TM2  (hard)  polarization 
case. 


4.7  Conclusions 


This  chapter  presents  a  TD-UTD  formulation  for  the  scattering  from  a  smooth  per¬ 
fectly  conducting  convex  surface.  The  TD-UTD  impulse  response  for  the  field  in  the 
presence  of  a  convex  smooth  surface  may  be  written  as 
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Figure  4.4:  Comparison  of  the  TD-UTD  solution  and  an  eigenfunction  reference 
solution  (eigen)  which  is  transformed  into  time  domain  with  the  IFFT  after  the 
window  in  Figure  4.3  is  applied.  The  far  zone  backscatter  from  a  2-D  circular  cylinder, 
for  TE*  (hard)  polarization  and  the  radius  is  r  =  1  meter. 


Time  (ns) 

Figure  4.5:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (at  bistatic  angle  ip  =  90°) 
from  2-D  cylinder,  for  TEZ  (hard)  polarization  and  radius  r  =  1  meter. 
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Figure  4.6:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (at  bistatic  angle  ip  —  135°) 
from  2-D  cylinder,  for  TE2  (hard)  polarization  and  radius  r  =  1  meter. 


Time  (ns) 


Time  (ns) 

Figure  4.7:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (at  bistatic  angle  ip  =  165°) 
from  2-D  cylinder,  for  TE2  (hard)  polarization  and  radius  r  —  1  meter. 


86 


Figure  4.8:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (at  bistatic  angle  if)  =  175°) 
from  2-D  cylinder,  for  TEZ  (hard)  polarization  and  radius  r  =  1  meter. 
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Figure  4.9:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  backscatter  from  2-D  cylinder,  for  TMZ  (soft) 
polarization  and  radius  r  =  1  meter. 


88 


Time  (ns) 

Figure  4.10:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (bistatic  angle  ^  =  90°)  from 
2-D  cylinder,  for  TMZ  (soft)  polarization  and  radius  r  =  1  meter. 
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Figuje  4.11:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (bistatic  angle  ip  =  135°) 
from  2-D  cylinder,  for  TM2  (soft)  polarization  and  radius  r  =  1  meter. 
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Figure  4.12:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (bistatic  angle  ip  =  165°) 
from  2-D  cylinder,  for  TM*  (soft)  polarization  and  radius  r  =  1  meter. 
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Figure  4.13:  Comparison  of  the  TD-UTD  versus  a  reference  eigenfunction  solution 
(“eigen”)  which  is  transformed  into  the  time  domain  with  the  IFFT  after  the  window 
in  Figure  4.3  is  applied.  The  far  zone  bistatic  scattering  (bistatic  angle  rj;  =  175°) 
from  2-D  cylinder,  for  TM*  (soft)  polarization  and  radius  r  =  1  meter. 
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where  the  impulse  response  of  the  incident  GO  field  is 

E’;  ( t )  =  E*  A(s’)  s  (t  -  s{/c)  (4.51) 

with  Eg  being  an  initial  value  vector  which  is  a  constant  with  respect  to  time  (or 
frequency)  and  contains  the  polarization  information  of  the  incident  field.  Note  that 
the  field  in  (4.51)  is  from  the  source  and  travels  to  the  observer  along  a  straight  line 
of  length  s'.  The  fields  incident  at  the  point  of  reflection  Qr  and  the  attachment 
point  Qi  have  the  same  form  as  in  (4.51).  The  TD-UTD  generalized  reflected  field  is 

Ef  (()  =  Ej  ■  [ft,  (rr)eLe±+  %h  W^ejj]  Ar(S')A(»')  (4.52) 

where  rr  =  t  —  s' /c  —  sT/c,  and  the  TD-UTD  surface  diffracted  field  is 

E;  (<)  =  E‘  •  [is  (rd)b1b2+  Vh  (r(,)n1n2]  A^A^)  (4.53) 

where  =  t  —  s'  jc  —  sd/c.  Note  that  the  variable  s'  in  (4.52)  is  the  distance  from  the 
source  to  the  reflection  point  Qr  but  the  variable  s'  in  (4.53)  is  the  distance  from  the 
source  to  the  point  of  attanchment  Q\.  The  TD-UTD  coefficients  7 ZSfh  and  Vath  can  be 

-f  p 

found  in  (4.48)  and  (4.49),  respectively.  The  special  function  Fah  (E,t),  used  in  the 
uniform  coefficients  7ta,h  and  can  be  calculated  using  the  numerical  algorithm 

in  Appendix  E  for  all  values  of  E  and  Im  ( t )  >  0.002|E3|  (or  Im  ( t )  >  0.002|(EI/)3|). 
Finally,  the  impulse  response  formulas,  denoted  by  the  subscript  /,  may  be  convolved 
with  an  excitation  in  closed  form  as  in  (3.80)  and  (3.81)  when  the  excitation  has  the 
form  in  (3.77)  and  (3.78). 
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Chapter  5 
Conclusions 


This  report  presents  a  development  of  a  time  domain  version  of  the  uniform  geometri¬ 
cal  theory  of  diffraction  (TD-UTD).  The  analytic  time  transform  (ATT)  is  formulated 
and  is  shown  to  be  very  useful  in  deriving  the  TD-UTD  solution  to  the  scattering  from 
an  arbitrary  curved  wedge  and  also  from  an  arbitrary  smooth  convex  surface,  both 
of  which  configurations  are  perfectly  conducting.  These  two  types  of  configurations 
constitute  a  very  important  set  because  many  complex  geometries  can  be  built  up 
from  these  configurations  as  has  already  been  demonstrated  in  the  frequency  domain 
UTD  based  ray  approach.  This  same  ray  picture  is  valid  in  the  TD-UTD  developed 
here,  and  therefore  the  TD-UTD  retains  all  the  advantages  of  the  frequency  domain 
UTD. 

The  ATT  is  defined  and  its  properties  are  discussed  in  Chapter  2.  An  important 
observation  is  that  the  ATT  may  be  expressed  as  a  one-sided  forward  Laplace  trans¬ 
form  when  the  roles  of  time  and  frequency  are  interchanged  as  expressed  in  (2.13). 
This  is  an  important  observation  since  it  implies  that  the  abundant  math  texts  and 
tables  on  the  one-sided  Laplace  transform  may  be  applied  in  deriving  results  for  the 
ATT.  Another  important  result  presented  in  Chapter  2  is  the  ATT  applied  to  an 
asymptotic  high  frequency  power  series  expansion,  which  is  useful  when  deriving  ex¬ 
pressions  for  the  numerical  computation  of  special  functions.  It  is  noted  that  the 
ATT  can  not  simply  be  applied  to  an  asymptotic  frequency  domain  power  series 
term-by-term,  nonetheless  more  rigorous  results  are  presented  in  Chapter  2.  When 
an  impulse  response  is  represented  by  an  analytic  time  function  (obtained  by  the 
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ATT),  the  convolution  with  an  excitation  pulse  function  can  be  performed  in  closed 
form  as  in  (2.70)  when  the  frequency  response  of  the  excitation  is  modeled  as  a  series 
of  exponential  functions  as  in  (2.67). 

The  TD-UTD  formulation  for  the  scattered  field  from  a  general  curved  wedge  when 
the  wedge  is  illuminated  by  an  astigmatic  wavefront  is  presented  in  Chapter  3.  The 
TD-UTD  impulse  response  for  the  curved  wedge  is  obtained  by  an  application  of  the 
ATT  to  the  frequency  domain  UTD  formulation.  This  TD-UTD  impulse  response  may 
then  be  convolved  with  a  more  general  time  pulse  excitation  where  the  convolution 
may  be  performed  efficiently  as  described  in  Chapter  2.  In  general,  the  TD-UTD 
impulse  response  for  the  curved  wedge  may  be  expressed  as 

E T°  ( t )  =E /  (t)Ui+  E/  (t)Ur+  E/  (t)+  Ef  ( t )  (5.1) 

where  the  spatial  unit  step  functions  [/,  and  Ur  are  1  on  the  lit  side  of  the  incident 
shadow  boundary  and  the  reflection  shadow  boundary,  respectively,  and  0  otherwise. 
The  TD-UTD  incident  GO  field  E/  (t)  is  in  (3.10)  and  the  TD-UTD  reflected  field 
E/  (t)  is  in  (3.14).  The  first  order  TD-UTD  edge  diffracted  field  E/  (t)  is  in  (3.16) 
where  the  TD-UTD  dyadic  diffraction  coefficient  is  defined  in  (3.40)  and  (3.41).  The 
higher  order  edge  diffracted  field  called  the  slope  diffracted  field  E/*  (t)  is  in  (3.56) 
where  the  TD-UTD  dyadic  slope  diffraction  coefficient  based  on  Hwang’s  frequency 
domain  UTD  slope  result  is  given  in  (3.57)  and  (3.58)  while  the  slope  diffraction 
coefficient  based  on  Veruttipong’s  frequency  domain  UTD  result  is  given  in  (3.64) 
and  (3.65).  The  two  versions  of  the  slope  diffraction  coefficient  are  derived  in  detail 
in  Appendix  C. 

Chapter  4  presents  a  TD-UTD  formulation  for  the  scattering  from  smooth  con¬ 
ducting  convex  surfaces.  Only  the  scattering  case  where  the  source  and  observer  are 
removed  from  the  surface  is  analyzed  in  Chapter  4,  but  the  radiation  and  coupling 
cases  could  be  handled  in  a  similar  fashion.  The  TD-UTD  impulse  response  for  the 
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field  in  the  presence  of  a  convex  smooth  surface  may  be  written  as 

E/  (t)+  E/r  ( t )  lit  region 

E r  (*)  =  (5-2) 

E /  (t)  shadow  region 

where  the  impulse  response  of  the  incident  GO  field  is 

E/ (t)  =  KMs')  6(t- sih)  (5.3) 

with  Eq  being  an  initial  value  vector  which  is  a  constant  with  respect  to  time  (or 
frequency)  and  contains  the  polarization  information  of  the  incident  field.  Note  that 
the  field  in  (5.3)  is  from  the  source  and  travels  to  the  observer  along  a  straight  line 
of  length  s'.  The  TD-UTD  generalized  reflected  field  E/r  (t)  is  given  in  (4.52)  and 
the  TD-UTD  surface  diffracted  field  E/  ( t )  is  in  (4.53).  The  surface  diffracted  field 
E/  (t)  may  also  exist  in  the  lit  region  for  closed  surfaces  when  multiple  encirclements 
are  included.  In  the  same  manner  as  the  TD-UTD  wedge  scattering  formulation,  the 
impulse  response  formulas  for  the  smooth  surface  scattering,  denoted  by  the  subscript 
/,  may  be  convolved  with  an  excitation  in  closed  form  as  in  (3.80)  and  (3.81)  when 
the  excitation  has  the  form  in  (3.77)  and  (3.78). 

The  development  of  the  TD-UTD  is  important  for  many  reasons.  There  are  many 
problems  of  engineering  interest  for  which  the  more  robust  time  domain  numerical 
methods  (such  as  the  finite  difference  time  domain  (FDTD),  or  the  time  domain 
integral  equations  (TD-IE)  method)  are  not  well  suited,  because  the  geometries  are 
too  large  in  terms  of  pulse  width.  Also,  the  TD-UTD  can  provide  an  approximate 
impulse  response  which  may  provide  more  interesting  information  about  a  scatterer 
or  an  antenna,  since  the  impulse  response  is  essentially  independent  of  excitation 
waveform  shape.  In  addition,  the  development  of  the  TD-UTD  has  requires  the  use 
of  some  novel  techniques  such  as  the  analytic  time  transform  (ATT)  and  so  forth, 
there  are  many  other  research  topics  which  the  development  of  the  TD-UTD  and  the 
ATT  may  have  a  positive  benefit. 

An  important  observation  from  the  development  of  the  TD-UTD  is  that  the  early 
time  asymptotic  results  obtained  from  the  high  frequency  asymptotic  frequency  do- 
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main  solutions  are  not  necessarily  causal.  This  is  understood  now,  because  the  early 
time  results  correspond  to  scattering  from  localized  portions  of  an  object,  and  the 
shifted  times  about  which  they  are  asymptotically  valid  is  NOT  the  so  called  “turn 
on  time”  for  the  general  geometry.  This  is  an  important  observation,  because  it  rules 
out  using  an  inverse  Laplace  transform  for  solving  general  problems  since  the  inverse 
Laplace  transform  only  provides  one-sided  time  functions.  Although,  in  principle  a 
two-sided  Laplace  transform  or  a  Fourier  transform  could  be  used,  the  ATT  developed 
in  Chapter  2  is  much  more  convenient.  The  analytic  time  signals  obtained  from  the 
ATT  provide  a  very  compact  general  solution  which  can  be  convolved  with  a  broad 
class  of  excitation  functions  in  closed  form. 

The  TD-UTD  formulations  developed  in  this  report  for  a  general  curved  wedge 
including  slope  diffraction  and  for  the  diffraction  from  smooth  convex  curved  surfaces 
(scattering  case)  represent  a  substantial  progress  towards  achieving  a  full  TD-UTD 
framework  which  is  just  as  useful  as  the  frequency  domain  UTD.  Nonetheless,  there 
is  a  lot  of  research  which  still  must  be  done  to  extend  the  usefulness  of  the  TD- 
UTD  even  farther.  First  of  all,  work  should  be  done  on  improving  the  efficiency  of 
computing  the  TD-UTD  transition  functions  (related  to  the  Fock  type  function)  in  the 
smooth  surface  diffraction  coefficient.  Although  the  numerical  algorithm  described  in 
Appendix  E  is  adequate  for  most  purposes,  it  is  likely  that  there  is  a  more  efficient 
method  for  the  calculation  since  accuracy  not  efficiency  was  the  primary  concern  in 
the  development  of  the  numerical  algorithm  in  Appendix  E.  Next,  there  are  still  some 
canonical  (localized)  geometries  to  be  analyzed  such  as  the  coupling  and  radiation 
problems  for  antennas  located  on  smooth  convex  surfaces,  diffraction  from  a  vertex, 
a  general  equivalent  current  type  of  formulation  for  observations  near  fau sties,  an 
impedance  wedge,  and  so  on.  In  particular,  the  development  of  the  TD-UTD  solution 
to  the  coupling  and  radiation  on/from  smooth  convex  surfaces  should  essentially 
follow  the  approach  used  for  the  scattering  case  considered  in  this  report. 

In  addition  to  the  obvious  potential  research  topics  concerning  the  TD-UTD  men¬ 
tioned  above,  there  are  some  other  potential  research  problems  related  to  the  TD- 
UTD.  One  of  these  topics  is  using  the  TD-UTD  as  the  guiding  principle  in  the  de- 
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velopment  of  new  parametric  models  used  in  the  analysis  of  radar  scattering  data 
[5,  6].  There  are  some  inherent  benefits  of  formulating  parametric  scattering  models 
in  the  time  domain  (or  more  generally  the  image  domain)  instead  of  the  frequency 
domain  since  for  scatterers  large  in  terms  of  wavelength  (or  pulse  width)  the  terms 
in  the  parametric  model  naturally  become  isolated  in  time.  Another  potential  use  of 
the  TD-UTD  is  in  the  development  of  a  hybrid  numerical  method  combining  a  robust 
numerical  approach  such  as  the  finite  difference  time  domain  (FDTD)  or  time  domain 
integral  equations  (TD-IE)  with  an  asymptotic  ray  based  methodology.  Although, 
it  is  not  currently  known  how  a  hybrid  time  domain  numerical  method  such  as  this 
would  perform,  it  is  felt  that  the  TD-UTD  formulation  as  presented  in  this  report 
would  be  helpful  in  its  development. 

Recently,  there  has  been  a  great  interest  in  a  new  type  of  transform  called  the 
wavelet  transform  [54,  55].  There  are  fundamentally  two  types  of  wavelet  transforms; 
the  discrete  wavelet  transform  (DWT)  and  continuous  wavelet  transform  (CWT).  It 
appears  that  the  DWT  is  more  popular  because  there  have  been  recent  developments 
in  fast  algorithms  for  the  computation  of  the  DWT.  Nonetheless,  it  turns  out  that  the 
second  time  derivative  of  the  ATT  defined  in  (2.11)  is  a  continuous  wavelet  expansion. 
See  Kaiser  [56]  to  understand  this  relationship.  In  particular,  Chapter  9  of  [56] 
“Introduction  to  Wavelet  Electromagnetics”  is  very  helpful.  To  utilize  this  CWT, 
there  may  be  further  modifications  necessary  to  compensate  for  the  poor  behavior  of 
the  asymptotic  high  frequency  results  at  the  low  frequencies,  but  these  modifications 
should  be  trivial.  Because  of  the  observation  that  the  ATT  is  actually  a  CWT,  the 
TD-UTD  which  uses  the  ATT  may  be  applied  to  compute  time-scale  or  time-frequency 
distributions  of  various  diffraction  phenomena  directly  or  essentially  in  closed  form. 
It  is  not  known  if  looking  at  the  TD-UTD  formulation  from  a  wavelet  transform 
perspective  would  be  helpful,  but  there  may  be  future  research  opportunities  here. 

In  summary,  the  current  development  of  the  TD-UTD  will  likely  impact  many 
areas  of  research  beyond  just  the  analysis  of  electromagnetic  radiation  and  scattering 
phenomena,  because  of  the  novel  use  of  the  ATT  and  also  the  increased  insight  into 
transient  electromagnetics. 
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Appendix  A 

Synthesis  of  an  Arbitrary  Finite 
Energy  Time  Pulse  for  Efficient 
Convolution 


This  appendix  describes  one  practical  method  of  obtaining  the  parameters  in  the 
model  of  the  excitation  waveform  in  (2.68)  or  (3.78).  As  discussed  earlier,  the  advan¬ 
tage  of  using  the  model  in  (2.68)  or  (3.78)  is  that  the  convolution  of  this  model  with 
the  TD-UTD  impulse  response  can  be  performed  in  closed  form  as  shown  in  (2.70) 
or  (3.80). 

Assume  that  the  constants  {An}  are  real  and  an  =  a!n  +  jet"  where  a'n  and  a"  are 
real  for  n  —  1, . . . ,  N.  Now,  the  real  time  waveform  constructed  by  the  model  (using 
the  series  in  (2.68)  or  (3.78))  is 


N 


m=  ^  f(o  =  - E 


AnQ!n 


*  ^1  (t  -  a'n)2  +  <2 


for  Im  t  =  0 


(A.l) 


Assume  that  the  entire  excitation  pulse  F^t)  is  contained  in  the  time  interval  0  < 
t  <  a  (or  at  least  approximately  contained  within  this  interval).  Notice  that  the 
parameter  a"  controls  the  position  of  the  expansion  term  on  the  t  axis  while  a'n 
controls  the  “width”  of  the  term.  It  would  be  convenient  to  select  explicit  values 
(which  may  depend  on  the  index  n)  for  the  parameters  a!n  and  a"  so  that  the  solution 
for  the  An  parameters  will  be  a  "system  of  linear  equations”  type  of  problem.  Let 
the  expansion  functions  be  equally  spaced  along  the  t  axis  and  choose 

a"  =  (n_1)A^l  (A-2) 
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for  n  =  1  Let  the  “width”  of  the  expansion  functions  be  constant 

for  all  n)  and  enforce 


a'  > 


1  a 
21V  —  1 


=  a 


to  make  sure  the  expansion  functions  are  wide  enough.  Typically,  one  can  choose 
a'  =  0.75 (a/(N  —  1))  or  a!  =  1.0 (a/(JV  —  1))  for  satisfactory  results.  Note  that 
choosing  an  a'  value  that  is  too  large  may  make  the  expansion  approximation  have  a 
large  “tail”  for  times  where  the  excitation  waveform  is  actually  zero,  on  the  other  hand 
choosing  a1  too  small  may  cause  too  much  “ripple”  in  the  expansion  approximation. 
Now  sample  the  excitation  waveform  F^t)  at  M  points  defined  by 


tm  =  (m-  1) 


a 

M-  1 


for  m  =  1,2 ,M  and  M  >  N.  The  expansion  in  (A.l)  is  set  equal  to  the  excitation 
waveform  F^t)  at  these  sample  points  to  form  an  overdetermined  system  of  linear 


equations 

[C][A\=[B] 

(A.3) 

where 

[C)mn 

1  oj, 

*  (tm  -  O2  +  a'n 

(A.4) 

\B]m 

=  Fo(tm) 

(A.5) 

=  An 

(A.6) 

One  can  solve  the  overdetermined  system  of  equations  in  (A. 3)  in  a  least  squares  sense 
by  using  the  singular  value  decomposition  (SVD)  [57].  Or,  a  simpler  but  less  robust 
approach  to  solve  (A. 3)  in  a  least  squares  sense  is  to  form  the  moment  equations 


\C]T[C]\A]  =  [CHS] 


(A.7) 


and  then  use  Cholesky  decomposition  on  (A.7)  to  solve  for  the  unknown  vector  [A] 
[57].  The  solution  of  these  equations  provides  the  {An}  parameters  and  all  the  pa¬ 
rameters  in  the  expansion  are  now  specified. 
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Appendix  B 

Exact  Time  Domain  Field  of  an 
Arbitrarily  Oriented  Electric  or 
Magnetic  Dipole  in  the  Presence 
of  a  Perfectly  Conducting  Wedge 


This  appendix  summarizes  an  exact  time  domain  solution  for  the  electromagnetic  field 
from  an  arbitrarily  oriented  electric  or  magnetic  dipole  in  the  presence  of  a  perfectly 
conducting  wedge.  The  solution  in  this  appendix  is  based  on  Felsen’s  work  [22].  In 
this  appendix,  a  on  top  of  a  variable  indicates  that  it  is  a  frequency  domain 
quantity.  An  eJUJf  time  convention  is  used  for  the  frequency  domain  quantities  in  this 
appendix. 

The  coordinate  system  used  here  is  shown  in  Figure  B.l.  Without  loss  of  generality 
the  source  is  assumed  to  lie  in  the  z  —  0  plane  and  the  observer  may  be  out  of  this 
plane.  The  wedge  angle  a  is  assumed  to  be  less  than  180  degrees,  but  the  solution 
for  the  general  case  is  given  by  Felsen  [22].  The  more  general  case  simply  includes 
more  image  contributions,  whereas  for  the  a  <  180  case  there  is  at  mosj  one  image 
contribution  to  the  geometrical  optics  field  corresponding  to  the  reflected  ray. 

B.l  Electric  Dipole  Case 

Electromagnetic  fields  can  be  expressed  in  terms  of  vector  potentials.  Felsen  chooses 
to  use  Hertz  vector  potentials  in  [22]  so  they  are  also  used  here.  The  electric  field  E 
and  magnetic  field  H  can  be  expressed  in  terms  of  the  electric  Hertz  vector  potential 
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Figure  B.l:  An  electric  (or  magnetic)  dipole  radiating  in  the  presence  of  a  PEC  wedge. 
The  dipole  is  in  the  z  =  0  plane  but  the  observer  may  be  out  of  this  plane. 


n  by  [58] 


d2 

E  =  W-Tl-fie— n 

(B.l) 

K  =  eVx^n 

(B.2) 

in  the  time  domain,  or 

e  =  v  v  *  n  -|-  u)2  fi€  n 

(P-3) 

H  =  jtveV  x  ft 

(B.4) 

in  the  frequency  domain.  The  space  surronding  the  wedge  is  assumed  to  be  a  linear  ho- 

mogeneous  isotropic  nonconducting  medium,  and  the  permittivity 

e  and  permeability 

fi  are  constants  with  respect  to  space  and  time  (or  frequency).  The  electromagnetic 

fields  satisfy  Maxwell’s  Equations 

8E  dP 

(B.5) 

VxH  =  £aT+9T 

_  _  aH 

VxE  =  ~'iHi 

(B.6) 

V-E  =  -iv-P 

e 

(B.7) 

V-H  =  0 

(B-8) 
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where  P  is  an  electric  dipole  moment  density  which  is  related  to  the  electric  current 
density  J  by 

J  -  <B-9> 


The  Hertz  electric  vector  II  satisfies  the  wave  equation 

v2n-„4n  =  -ip. 


(B.XO) 


Notice  that  the  Hertz  electric  vector  II  is  related  to  the  conventional  vector  potential 


A  by  [59] 


or  in  the  frequency  domain 


A  =  fie- 


A  =  jujfiell. 


(B-ll) 


(B.12) 


In  the  far  zone,  defined  here  as  the  limit  r  — »  oo,  the  Hertz  vector  can  be  written 


e-ifcr _ 

n  nsj  —  - f p(6,  <j>)  for  r  — ►  oo 

r 


and  the  electric  field  becomes 


u>2  e~jkr 

E  ~  — - r  x  r  x  f p(9,<f>)  for  r  — ►  oo 

Cr  T 


where  c2  =  l/(/xe).  In  the  time  domain,  these  far  zone  relationships  become 

II  ~  <j),t  —  r/c)  for  r  — >  cxd 

r 

1  d 2 

E  ~  x  ^  x  “  r/c)  for  t  >  oo. 

clr  ot i 


(B.13) 


(B.14) 


(B.15) 

(B.16) 


Usually,  the  1/r  factor  and  r/c  time  delay  in  (B.15)  or  (B.16)  are  not  included  in  far 
zone  calculations. 

The  “impulse  response”  in  this  appendix  is  defined  in  terms  of  the  dominant  term 
of  the  electric  field  incident  on  the  wedge,  therefore  to  obtain  the  impulse  response 
let  the  electric  dipole  moment  density  have  a  ramp  time  dependence  as  in 


P(f)  =  p£(r  —  v')tu[t) 

or  equivalently  let  the  electric  current  element  be 

J(t)  =  p£(r  —  r')tt(t) 


(B.1T) 


(B.18) 


where  u(t)  is  the  unit  step  function. 
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The  Geometrical  Optics  Field 

* 

It  is  assumed  here  that  the  orientation  of  the  dipole  p  is  contained  in  the  z  =  0  plane 
(i.e.  a  transverse  source).  The  solution  for  a  longitudinal  source  can  be  found  in 
Felsen’s  book  [34]. 


The  Hertz  vector  can  be  written  as 


n  =  ns°  +  nd 


(B.19) 


where  11®°  corresponds  to  the  geometrical  optics  (GO)  field,  while  IId  corresponds 
to  the  diffracted  field.  When  the  wedge  angle  a  is  less  than  180  degrees,  the  GO 
Hertz  vector  can  be  written  as  two  terms 

ng0  _  1 

47re  |r  —  r  | 

+P.7j—(‘ - |r -  |:l)“(t,~  |r - r:l  W  -\t- #1)  (B.20) 

47t  e  |r  —  r,| 

where  the  orientation  of  the  dipole  is 


p  =  x  cos  u  +  y  sin  v 


(B.21) 


and  the  orientation  of  the  image  is 


p,  =  —x  cos  v  +  y  sin  v. 


(B.22) 


The  position  of  the  dipole  in  cylindrical  coordinates  is  r'  =  0),  while  the 

position  of  the  image  source  is  r|  =  (p',  $,0)  where  <$>\  —  2iv  —  <f>' . 

To  obtain  the  far  zone  (r  —>  o o)  expressions  use 


r-rol 


r  in  magnitude 

r  —  f  •  ro  in  phase  (or  time  delay) 


(B.23) 


where  ro  may  be  r'  or  r'.  Using  the  above  far  zone  approximation,  the  Hertz  vector 


becomes 


,  1  (r  +  r  •  t'/c)u(t  +  r  •  r'/c)  ,  ,  A,iS 

p  — i - '-U{T  -\(f)-(f>\) 

47re  r 

1  (x  +  f.r;/cKr  +  f.r;/c)  _  _  ^  (B24) 

47re  r 
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where  t  =  t  —  r/c.  Using  the  above  result  together  with  (B.15)  and  (B.16),  the  far 
zone  geometrical  optics  electric  field  can  be  written  down  from  inspection  as 

ZQ 


Eg° 


rxrx  p- - 6(t  +  r •  t'/c)u(tc  —\<j>  —  <j>'\) 

47 xrc 

Irxrx  p  +  r  •  r'/c)u(7r  -\<t>-  $|) 

47rrc 


(B.25) 


where  Zq  =  \J p/e- 


The  Diffracted  Field 

The  Hertz  vector  corresponding  to  the  diffracted  field  11^  as  presented  by  Felsen  is 
repeated  here  for  convenience  [22].  Also,  the  simplification  of  11^  for  an  observer  in 
the  far  zone  and  the  electric  field  in  the  far  zone  is  derived. 

When  the  orientation  of  the  dipole  p  is  transverse  to  z  the  solution  is  obtained 
from  only  the  transverse  components  of  the  Hertz  vector  II  =  xllx  +  yHy.  Felsen 
gives  [22] 

TTrf  C  Re  — ' i(3)]  I  V  ('R 

aPn-  =  -  5S - - “(‘  - r/c)  (B-26) 

where  the  double  time  derivative  on  n!*  „  occurs  because  Felsen  uses  the  source  P  = 
pS(r  —  r')6(t )  which  is  the  double  time  derivative  of  the  source  used  here.  The  other 
parameters  in  (B.26)  are 


=  cosh 


-i 


cH 2  —  pl  —  pf*  —  z 


,/2 


.21 


2  PP' 

T  =  yf{p  +  p')2  +  ~Z2 


(B.27) 

(B.28) 


and 


7T 


Ax(<t>,<j>',w)  =  {[Q i{<l> -<f>',w)-Q ^(f) +  (/>', w)]  cosu 

~  [Qi{(j>  -  <j>',  w)  +  Q2(<I>  +  ™)]  sin  v } 


(B.29) 


7 T 


Ay{(f>,(t>' ,w)  =  {[Q2((/>-<t>',w)-Q2(<l>  +  <t>,,w)}cosv 

+  ~  <f>' ,™)  +  Qi{<f>  +  (f>',w)}  sin  v} 


(B.30) 
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with 


Qi(ip}w)  —  cos^  —  w  —  7r)cot 


—  cos  (i/)  +  w  +  7r)  cot 


Ip  —  W  —  7r\ 

2a/n  ) 

( Ip  +  W  +  7T 


2ol/tt 


(B.31) 


Q2(ip,w)  —  sin(-0  —  w  —  7r)  cot 


sin  (ip  +  w  +  tt)  cot 


l/j  —  W  —  7r\ 

2a/7r  J 

/  if)  +  W  +  7T 


2a/7r 


(B.32) 


To  obtain  the  far  zone  expressions,  the  following  far  zone  approximation  is  used 


C2t2  _p2_  p'2  _  z2  CT 

-  rsj  - 

2  pp1  p1  sin  6 


for  r  — >  oo 


(B.33) 


where 


r  =  t-r/c 


(B.34) 


and  notice  that 


sinh  cosh 


p  sin  i 


/c2r2  _  ^/2  sjn2  0. 


So,  in  the  far  zone,  the  Hertz  vector  y  becomes 


&_T\d  „  ~c  Re  [A*M  £± zj£& 

dt 2  *,tf  47 r2er  ^.2  _  ^  sin2  # 


u.  I  r  —  —  sin  6 


(B.35) 


(B.36) 


where 


/ 3  ~  cosh 


[p'  sin0j 

Now  the  far  zone  diffracted  electric  field  can  be  written  down 
Ed  ~  -f^-r  x  r  x  {x  Re  [Ax((j>,<f>\  -j/3)\ 

47 TlT 


+  y  Re  [Ay(<f>,  4>',-jP)}} 


u(t  —  ^  sin  8^j 
Jc2t2  —  p'2  sin2  6 


(B.37) 


(B.38) 
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B.2  Magnetic  Dipole  Case 


Once  again,  the  electromagnetic  fields  can  be  expressed  in  terms  of  vector  potentials. 
Felsen  chooses  to  use  Hertz  vector  potentials  in  [22]  so  they  are  also  used  here.  The 
electric  field  E  can  be  expressed  in  terms  of  the  magnetic  Hertz  vector  potential  IIm 
by  [58] 

E  =  -/vxnm  (B.39) 

ot 

in  the  time  domain,  or 


E  =  -ju/fiV  x  nni 


(B.40) 


in  the  frequency  domain.  The  electromagnetic  fields  satisfy  Maxwell’s  Equations 


V  x  H 

V  x  E 
V-E 


dE 
'  dt 


-IL 


SH 

dt 


VH  =  -VM 


(B.41) 

(B.42) 

(B.43) 

(B.44) 


where  M  is  a  magnetic  dipole  moment  density  which  is  related  to  the  magnetic  current 
Kby 

K  =  /.  jjM.  (B.45) 

The  magnetic  Hertz  vector  IIm  satisfies  the  wave  equation 

Q 2 

V2nm  -  nei-nm  =  -M.  (B.46) 

ot 1 

Notice  that  the  Hertz  magnetic  vector  IIm  is  related  to  the  conventional  vector  po¬ 
tential  F  by  [59] 

dUm 

P  =  (B.47) 

or  in  the  frequency  domain 

F  =  jidfieir71.  (B.48) 
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In  the  far  zone,  defined  here  as  the  limit  r  — »  oo,  the  Hertz  vector  can  be  written 
as 

e-ikr 

IIm  ~ - gm(0,  (f>)  for  r  — >  oo  (B.49) 

r 

and  the  electric  field  becomes 

2  ~  jkr 

E  ~  —  - r  xgm(^)  for  r  — *  oo.  (B.50) 

c  r 

In  the  time  domain,  these  far  zone  relationships  become 

nm  ~  -gm(^,^,t-r/c)  for  r  — >  oo  (B.51) 

r 

q2 

E  ~  —  —  r  x  —  r/c)  for  r  — >  oo  (B.52) 

cr  ot 1 

Usually,  the  1/r  factor  and  r/c  time  delay  are  not  included  in  the  far  zone  calculations. 

The  “impulse  response”  in  this  appendix  is  defined  in  terms  of  the  dominant  term 
of  the  electric  field  incident  on  the  wedge,  therefore  to  obtain  the  impulse  response 
let  the  magnetic  dipole  moment  density  have  a  ramp  time  dependence  as  in 

M(f)  =  —  m£(r  —  r')tu(t)  (B.53) 

t1 

or  equivalently  let  the  magnetic  current  element  have  a  step  time  dependence  as  in 

K(t)  =  m£(r  —  r')u(t)  (B.54) 

where  u(t)  is  the  unit  step  function. 

The  Geometrical  Optics  Field 

It  is  assumed  here  that  the  orientation  of  the  magnetic  dipole  rii  is  contained  in  the 
z  —  0  plane  (i.e.  a  transverse  source).  The  longitudinal  source  case  can  be  found  in 
Felsen’s  book  [34]. 

The  Hertz  vector  can  be  written  as 

nm  =  nmso  +  it"’11  (b.55) 
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corre- 


where  IIm’§0  corresponds  to  the  geometrical  optics  (GO)  field,  while  IT71’^ 
sponds  to  the  diffracted  field.  When  the  wedge  angle  a  is  less  than  180  degrees,  the 
GO  Hertz  vector  can  be  written  as  two  terms 

1  (t  -  |r  -  -  |r  -  r'|) 


nm,6°  = 


m 


47 Tfi 


r  —  r 


-u(rr  -  \<f>-<f>'\) 


izfl  lr_ril 


where  the  orientation  of  the  magnetic  dipole  is 


m  =  x  cos  v  +  y  sin  v 


(B.57) 


and  the  orientation  of  the  image  is 


m,  =  x  cos  v  —  y  sin  v. 


(B.58) 


The  position  of  the  magnetic  dipole  in  cylindrical  coordinates  is  r'  =  0),  while 

the  position  of  the  image  source  is  r(  =  (p' ,  </>',  0)  where  <$>\  =  27r  —  <f>' . 

To  obtain  the  far  zone  (r  — >  oo)  expressions  use 


r-  r0  ~  < 


r  in  magnitude 

r  —  r  •  ro  in  phase  (or  time  delay) 


(B.59) 


where  ro  may  be  r'  or  r'.  Using  the  above  far  zone  approximation,  the  magnetic 
Hertz  vector  becomes 

■IT-*0  ~  AJ_(^  +  fT7cMr  +  f-r-/c) 

47 r/x  r 

tf.|)  (B.60) 

47 Tfi  r 

where  t  —  t  —  r/c.  Using  the  above  result  together  with  (B.51)  and  (B.52),  the  far 
zone  geometrical  optics  electric  field  can  be  written  down  from  inspection  as 


E®°  ~  r  x  m- - S(r  +  r  •  r'/c)u(7r  —  \<f>  —  <j>'\) 

4tt  rc 

+r  x  m,— —  8(t  +  r  •  r^/c)u(7r  -  \<f>  -  #|). 
47T  rc 


(B.61) 
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The  Diffracted  Field 


The  magnetic  Hertz  vector  corresponding  to  the  diffracted  field  IIm’^  as  presented 
by  Felsen  is  repeated  here  for  convenience  [22].  Also,  the  simplification  of  IIm’^  for 
an  observer  in  the  far  zone  and  the  electric  field  in  the  far  zone  are  presented. 

When  the  orientation  of  the  magnetic  dipole  m  is  transverse  to  z  the  solution  is 
obtained  from  only  the  transverse  components  of  the  Hertz  vector  IIm  =  idl™ +yn™. 
Felsen  gives  [22] 


«V,d  = 

dt1  xu 


c  Re 


u(t  —  r/c) 


(B.62) 


47T2  pp'sinh/3 

where  the  double  time  derivative  on  H!?’^  occurs  because  Felsen  uses  the  source 
M  =  im£(r  —  r ')S(t)  which  is  the  double  time  derivative  of  the  source  used  here. 
The  parameters  in  (B.62)  are  the  same  as  in  the  electric  dipole  case,  except 


=  —  —  {[Qi(<f>  —  +  Qi{4>  +  (f>',w)]cosi/ 

-  [Q2{(f>  -  <t>,  w)  -  Q2(<t>  +  <£',  w)]  sin  I/}  (B.63) 

A^{(f>,(t>\w)  =  -^{[Q2(4>-4>',w)  +  QM+(t>',w)\™SV 

+  [Qi{<f>-  <f>',w)  -  Q\((f)  +  <t>  ,w)\  sin  v}  (B.64) 


The  far  zone  expressions  are  obtained  in  the  same  manner  as  in  the  electric  dipole 
case.  In  the  far  zone,  the  Hertz  vector  H™^  becomes 


92  nm,d  -c  Re  [A^v(^,^,-j/3)] 


-II"'’"  ~ 
dt 2  x'v  4t r2r 


Ic2t2 


where 


yf- 

/3  ~  cosh-1 


p'2  sin2  6 


u  I  t  —  —  sin  6 
c 


CT 


[p1  sin  9\ 


as  r  — >  oo,  and 

t  -t-  r/c. 

Now  the  far  zone  diffracted  electric  field  can  be  written  down 


Ec 


1 


47 r2r 


r  x  {x  Re  [A™(</>,  <j)' , -j/3)} 


u  (t  —  ^  sin  6 ) 
^/cV2”—" /9'2sin2~0 


(B.65) 

(B.66) 

(B.67) 


(B.68) 


as  r  — ►  oo. 
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Appendix  C 

Derivation  of  Two  TD-UTD  Slope 
Diffraction  Coefficients,  One 
Based  on  Veruttipong’s  Work  and 
the  Other  on  Hwang’s  Work 


This  appendix  derives  two  versions  of  the  TD-UTD  slope  diffraction  coefficient;  one 
version  is  based  on  a  frequency  domain  slope  diffraction  coefficient  derived  by  Hwang 
and  Kouyoumjian  [21,  4]  and  the  other  is  based  on  a  frequency  domain  slope  diffrac¬ 
tion  coefficient  derived  by  Veruttipong  and  Kouyoumjian  [20].  For  simplicity,  these 
two  results  will  be  referred  to  as  Hwang’s  version  and  Veruttipong’s  version  of  slope 
diffraction,  respectively.  The  TD-UTD  slope  diffraction  coefficient  for  either  case  is 
obtained  by  an  application  of  the  analytic  time  transform  (ATT)  described  in  Chapter 
2. 

Hwang  obtains  the  slope  diffraction  coefficient  by  taking  the  derivative  with  re¬ 
spect  to  <f>'  of  the  exact  straight  wedge  diffraction  coefficient,  where  this  exact  coeffi¬ 
cient  is  written  as  a  spectral  integral.  The  operations  of  integration  and  differentiation 
are  interchanged  and  the  new  spectral  integral  is  evaluated  asymptotically.  On  the 
other  hand,  Veruttipong  finds  the  slope  diffraction  coefficient  by  taking  the  derivative 
with  respect  to  (f>'  of  the  UTD  wedge  diffraction  coefficient,  where  the  UTD  coeffi¬ 
cient  is  obtained  by  an  asymptotic  evaluation  of  the  exact  diffraction  coefficient.  So, 
Veruttipong  takes  the  derivative  of  an  asymptotic  expression  whereas  Hwang  inter- 
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changes  the  derivative  and  integration  operations  first  and  then  does  the  asymptotic 
evaluation. 

This  appendix  derives  the  time  domain  slope  diffraction  coefficients  for  the  straight 
wedge,  but  these  slope  diffraction  coefficients  can  be  generalized  to  the  curved  wedge 
case  as  is  done  in  Chapter  3. 


C.l  Hwang’s  Slope  Diffraction  Coefficient 

Hwang’s  frequency  domain  slope  diffraction  coefficient  for  a  straight  wedge  can  be 
written  as 

Kh  =  Z  2  /o-  •  2  a  X) 

An V27T  sur  pa  m=1 


where 


S,/l  _ 


C[ 

Cs2'h 
C-'h 


3  “ 


Cs'h 


4  ~ 


+  csc2[(7 r  +  /3  )/(2n)] 
—  csc2[(7t  —  P~)/(2n)] 
±CSC2[(7T  +  /3+)/(2u)] 
-F  csc2[(7t  —  /3+)/(2n)] 


and 


with 


Fa  (®nt)^)  — 


W) 


P^Xm^-Jju,Ptr{  c  ) 


where  Ftr(z )  is  the  UTD  transition  function  defined  by 

too  .  „ 

Ftr(z)  =  2j^e3Z  [  e~3T  dr 

Jy/z 

The  F(xm, w)  function  can  also  be  written  as 


F(xm,u)  =  v'ir*meJ"*m/eerfc 


Ijux, 


as  is  done  in  Chapter  3.  The  xm  for  a  straight  wedge  are  defined  by 


(C.l) 


(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 

(C.7) 

(C.8) 


(C.9) 


=  La+(/3  ) 


(C.10) 
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®2  =  La~(P~)  (C.ll) 

®3  =  La+{(3+)  (C.12) 

®4  =  La~((3+)  (C.13) 

where  a±(/3)  =  2  cos2[(2n7riV±  —  /3)/2]  and  iV±  is  the  nearest  integer  solution  of 
2imN±  —  /?  =  ±7r.  The  L  parameter  for  a  straight  wedge  can  be  found  in  [4,  2]. 

Now  the  TD-UTD  slope  diffraction  coefficient  based  on  Hwang’s  coefficient  is 

i(‘)=TT^nTEtff,M  (C.14) 

4n2V27r  sm  p0  m=1 

where  Fs  (xm,t)  is  the  analytic  time  transform  of  Fs(xm,u>)  given  by 


Fs  (xm,t)  =  2xr 


+ 


which  can  be  easily  simplified  into 


H" 

F t 


2xm  e 


-jV/4 


(C.15) 


(C.16) 


'  v/^jl  + 

It  is  interesting  to  examine  some  special  cases.  Consider  the  straight  wedge  when 
it  is  illuminated  by  a  transient  pulse  which  is  a  real  time  impulse  S(t).  For  this  case, 
xm  >  0  and  we  are  only  interested  in  the  real  part  of  the  slope  transition  function  in 
(C.16).  First  Fs  (®m)f)  can  be  written  as 


Fs  (x, 


( 2Xrn\/c\  I 


)  t  +  xm/c 


(C.17) 


Now  the  real  time  function  is  the  real  part  of  (C.17)  evaluated  at  real  time  ( Im  t  =  0) 


F. s(®m)  ^)  — 


2®, 


R  ) 


Vi 


-u(t) 


(C.18) 


V*  J  t  +  Xm/c 

Notice  that  this  function  is  zero  at  t  =  0  and  has  a  peak  at  t  =  xm/c.  This  signifies 
that  the  impulse  response  of  the  slope  diffracted  field  for  the  straight  wedge  is  a 
continuous  time  function  and  has  a  maximum  after  the  time  of  arrival  of  the  wavefront 
(i.e.  the  time  corresponding  to  the  diffracted  ray  path  length).  If  the  observation  point 
is  far  from  the  shadow  boundary,  then  ®m  can  be  considered  much  greater  than 
ct  and  we  find 


2 r3/2  x 

Fs(xm,t)  ~  —=Vt  u(t)  for  -^—>00 

Y  7T  Ct 


(C.19) 
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C.2  Veruttipong’s  Slope  Diffraction  Coefficient 

Veruttipong’s  frequency  domain  slope  diffraction  coefficient  for  a  straight  wedge  can 
be  written  as 

jy  =  -1 _ ' 

a'h  4n2\/2Tr  sin2/30  . 

+A1Fvs(x  i,w)  +  BiFa(xuu) 

-A2Fvs(x2,u>)  -  B2Fa{x 2,w) 

±^43F„s(x3,o;)  ±  B3Fs(x3,uj) 


(C.20) 

where 

Ai  =  csc2((7t  +  f3~)/2n)  —  Bx 

(C.21) 

A2  =  csc2((tt  —  (3~)/2n)  —  B2 

(C.22) 

A3  =  csc2((7T  +  /3+)/2n)  -  I?3 

(C.23) 

A4  =  CSC2((7T  —  /3+)/2  n )  —  B4 

(C.24) 

and 

Bi  =  +  a+^_^cot((7r  +  ^  )/(2n)) 

(C.25) 

#2  =  .A?  cot((ir  /?  )/(2n)) 

a  ) 

(C.26) 

B3  =  +  a+^+^cot((7r  +  /5+)/(2n)) 

(C.27) 

B4  =  a-(^+^cot((7r  ^+)/(271)) 

(C.28) 

+/_.  «  2/27t  nN*  — /3\ 
a±(P)  =  2  cos  1  2  PJ 

(C.29) 

a±(/ 3 )  =  sin(27i7riV^:  —  /3)  (C.30) 


All  the  other  parameters  in  (C.20)  are  defined  in  the  previous  section  of  this  appendix. 
The  transition  function  Fa(xm, u>)  is  defined  in  (C.6)  and  its  analytic  time  transform  is 
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(C.16).  The  other  transition  function,  which  only  arises  in  Veruttipong’s  formulation, 


c3/2  . 

Fv j(®ni)W)  :  T  p. — F [r 


(“) 


(C.31) 


where  FtT(z)  is  the  UTD  transition  function  defined  in  (C.8).  The  Fvs{xm^u})  function 


can  also  be  written  as 


F(xm,u)  =  —  y/vx^ 

JijJ 


gjWm/cerfc  (  JtH 


(C.32) 


where  the  complement  of  the  error  function  is  defined  as 


2  r°° 

"{c{z)=j;L 


(C.33) 


The  analytic  time  transform  of  Fvs(xm,uj)  is  derived  as  follows.  First  notice  that 


-  c  ~ 

Fvs(xm,u)  =  —F(xm,u) 

ju> 


(C.34) 


where  F(xm,u)  is  defined  in  Chapter  3  for  use  in  the  first  order  diffraction  and  is  given 
in  (3.32)  or  (3.34).  The  analytic  time  transform  of  F(xm,u)  denoted  by  F  (*m,t)  in 
(3.38)  is  derived  in  Chapter  3.  It  is  repeated  here  for  convenience 


F  (xm,t)  = 


__2£xj1 


— - — - j  — ; - 

V? TCy/t(t  +  Xm/c)  t  +  Xm/c 


for  Im  t  >  0  where 


•E  rn  — 


*m  >  0 

s/  xm ,  xm  <  0 


(C.35) 


(C.36) 


and  Re  >  0.  The  following  property  of  the  analytic  time  transform  is  used 


_Ltf  („,)  A  H('l)  ( t )  +  Co 

JU3 


(C.37) 


where  H ^  ^  (t)  is  the  anti- derivative  (or  primitive)  of  H  (t)  and  Cg  is  a  constant 
which  must  be  determined.  Using  (C.37),  the  analytic  time  transform  of  Fva(xm,u ) 
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IS 


+  /_  ^  L  l 

FVs(xm,t)  -  jc^  \n\ 


V  V  i  Wi- jyjxm/cj 

=  j2c\ /^- In  (Vt  +  j y*m/c)  +  Co 
=  ;2cy^  Jin  |\/t  +  j\[x^fc 
+j  arg  +  i  \Jxmlcj  +  Co 


+  ln(t  +  a:m/c)  +  Co 


+J  arg  ^V^  +  iv*m/cJJ  +  C0  (C.38) 

where  arg(z)  is  the  argument  of  the  complex  number  z.  The  constant  Co  are  chosen 
such  that  Fvs  (xm,t)  =  0  when  t  =  0  (for  xm  >  0).  This  is  a  reasonable  requirement 
since  in  general 


+  d  + 

lim  juH(u))  —  lim  —j  H  ( t )  -  H  ( t ) 

(-jf)— ►  0  Ot 


(C.39) 


lim  juFvs(xm,u )  =  0 

W— >oo 

So,  the  analytic  time  transform  of  Fva(xm,w)  can  be  written  as 


L(xm,t)  =  jc^  In 


Wi  + 


J\^Xm/C 


(C.40) 


+  arg  (Vt  +  ~  (C-41) 


To  obtain  a  more  convenient  form  for  calculations,  use  ^/zm  =  —j\/—xm  (since 
-7r/2  <  arg(-*m)  <  37r/2)  and  \/i  =  e^^yj-jt  (since  -7r  <  arg(-jf)  <  i r)  to  arrive 


Fvs  (®mj^)  — 


\/-*m  f,  eiTfAy/-jct  +  y/^Xr 

A  [ 

+J  arg  (eJ7r/4^-jct  +  \Z-xm)  -  . 


(C.42) 
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Appendix  D 

Computation  of  the  Time  Domain 
Creeping  Wave  Function 

This  appendix  discusses  the  analytic  time  function  which  is  used  in  the  TD-UTD 
creeping  wave  analysis  in  Chapter  4.  The  layout  of  this  appendix  is  as  follows.  First 
the  creeping  wave  time  function  is  defined  and  some  of  its  properties  are  examined. 
Then  some  simple  approximate  formulas  for  small  time  values  or  large  time  values 
are  derived.  Finally,  an  efficient  numerical  algorithm  is  presented  to  compute  the 
creeping  wave  time  function  for  all  intermediate  time  values. 

D.l  Definitions  and  Properties 

The  creeping  wave  time  function  is  defined  as 

FcW  (o,0  =  -  r(ju,)-h-a^ejutdw  (D.l) 

77  Jo 

where  Im  t  >  0,  a  is  a  real  constant  and  a  >  0. 

By  taking  the  appropriate  derivatives  and  interchanging  the  differentiation  oper¬ 
ation  with  the  integration,  it  is  easy  to  show  that  Few  satisfies  the  following 

differential  equation 

Q 3  +  9  + 

^  Feu,  (a,  t)  +  —  Feu,  (a,  t)  =  0  (D.2) 

This  differential  equation  with  appropriate  boundary  conditions  could  be  used  to 
derive  a  series  expansion  in  t  or  a  series  in  a  of  Fcw  (a,  t),  although  this  approach  is 
not  used  here. 
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Some  special  values  of  the  creeping  wave  time  function  are 


+  ,  .  r(5/6)e--7,r/12 

Few  (a  =  0,0  =  W  ' 


d  +  .  . 

^  Feu,  (<M) 


o=0 


dn  + 

d^Fcw 


(a,t) 


o=0 


7r(-j05/6 

_  r(7/6)eJ7r/12 

ir(-jt)7/6 

r  eif?(2n_l) 


and 


Feu,  (<M  =  -oo)  =FCU,  (t  =  oo,  a)  =  0 


(D.3) 

(D.4) 

(D.5) 

(D.6) 


Friedlander  derived  the  following  early  time  approximation  by  evaluating  the  in¬ 
verse  Laplace  transform  integral  asymptotically  [7] 


-i  -ces^ 
s  6  e 


ak  exp  [— 2f  2(01/3)2] 
2-y/7T  t 


u(() 


(D.7) 


for  i  — >  0+ ,  and  u{t)  is  the  Heaviside  unit  step  function.  Notice  that  (D.7)  is  only 
valid  for  real  time  (i.e.  Im  t  =  0).  Because  of  the  relationships  between  the  Laplace 
transform  and  the  Fourier  transform,  (D.7)  is  an  approximation  of  the  real  part  of 
Few  (o,  0  in  (D.l).  Using  (D.7)  it  can  be  deduced  that 


[1  3  P 

— 2t"2(o/3)2 

“  - 1 - l 


«(() 


(D.8) 


for  Im  t  =  0  and  t  — >  0+.  Notice  that  (D.8)  goes  to  zero  as  t  -*  0+,  and  in  fact 
any  order  of  time  derivative  of  (D.8)  goes  to  zero  as  t  — >  0+.  Figure  D.l  shows  a 
comparison  between  a  numerical  calculation  of  the  creeping  wave  time  function  in 
(D.l)  and  the  early  time  asymptotic  result  in  (D.8).  Notice  that  only  the  real  part  of 
(D.l)  is  plotted  in  Figure  D.l.  Also  notice  that  the  peak  location  of  the  asymptotic 
result,  which  in  general  occurs  at  t  =  is  very  close  to  the  peak  location  of  the  exact 
result.  The  comparison  in  Figure  D.l  is  for  the  case  of  a  =  1  but  the  comparison 
is  essentially  the  same  for  other  values  of  a  (except  for  a  change  in  height  and  scale 
change  in  t)  because  of  the  scaling  properties  of  Few  (<M)  and  the  asymptotic  result 
in  (D.8). 
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The  creeping  wave  time  function  has  the  following  scale  change  properties. 


F,w  {a,  It)  =  -jTjl  ( a/t1/3,( ) 


(D.9) 


or 

Few  ( ba ,  t)  =  ^  Few  (<*,  t/b 3)  (D.10) 

where  b  >  0. 

An  early  time  power  series  representation  for  Few  (ai  t)  can  be  derived  by  substi¬ 
tuting 

e**  =  f  (D.ll) 

n=0  n! 

into  (D.l)  and  integrating  term  by  term  to  obtain 


+ 

Few 


(a,t) 


-3 j  ~  T(5/2  +  3n) 
7ra5/2  ^  n! 


(D.12) 


for  |t/a3|  — »  0.  This  type  of  early  time  series  is  sometimes  called  an  “expansion  by 
moments”.  Notice  that  this  series  gives  Re  [Few  (tM)]  =  0  f°r  \t/a3\  0.  So,  the 

series  in  (D.12)  is  only  useful  if  |£/a3|  is  “very  small”.  Also  notice  that  this  series 
is  divergent,  but  it  is  asymptotic  for  |t/a3|  — >  0.  This  type  of  asymptotic  series  is 
sometimes  called  “semi-convergent”  [60].  The  method  used  to  obtain  the  asymptotic 
series  in  (D.12)  is  valid  when  all  of  the  integrals  are  convergent  (see  page  134  in 
Bleistein  and  Handelsman’s  book  [31]). 

A  late  time  (or  large  \t/a3\)  series  representation  can  be  derived  by  expanding  the 
frequency  domain  function  in  the  integrand  of  the  transform  integral  of  (D.l)  into  a 
low  frequency  series  expansion  and  integrating  term  by  term.  This  is  done  with  the 
help  of 

e-“(^)1/3  =  V  ju;)"/3  (D.13) 

n—0  n! 

which  is  used  in  (D.l)  to  obtain 


+ 

Fe 


(a,t)  = 


o  3  T 


T(5/6  +  n/3)e,r>6  (-a)” 
n!  (—  it)"/3 


(D.14) 


The  series  in  (D.14)  is  convergent  for  t  ^  0  which  can  be  shown  with  some  effort 
by  using  the  ratio  test.  From  a  practical  stand  point  the  series  in  (D.14)  has  good 
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convergence  when  |t  |  >  a3.  Note  that  the  peak  of  Re  [Few  (a,  /)]  occurs  close  to 
t  =  a3/ 27,  as  seen  in  Figure  D.l,  which  is  significantly  smaller  than  t  =  a3.  So,  the 
series  in  (D.14)  is  only  useful  when  t  is  “very  large”. 


D.2  Numerical  Algorithm 

This  section  describes  an  efficient  algorithm  which  can  be  used  to  evaluate  the  integral 
in  (D.l)  numerically.  This  algorithm  can  theoretically  be  used  for  large  or  small  values 
of  |t|  (provided  that  Im  t  >  0).  But,  a  more  practical  approach  is  to  use  the  early 
time  series  in  (D.12)  when  |i|  <C  a3/ 27,  use  the  late  time  series  in  (D.14)  when 
\t\  a3  and  use  the  numerical  algorithm  for  all  other  |f|  values. 

The  numerical  algorithm  is  derived  by  first  applying  a  “change  of  variables”  on 
the  integral  of  (D.l).  Then  it  is  noticed  that  the  resulting  integrand  contains  a  slowly 
varying  factor  and  a  rapidly  varying  exponential  factor.  So,  the  slowly  varying  factor 
can  be  approximated  by  a  simple  polynomial  or  series  of  exponentials  and  then  the 
integration  is  performed  in  closed  form. 

The  following  change  of  variables 

z  =  (D.15) 


-v  r-*"8’  je-^dz 


is  used  to  change  (D.l)  into 

Few  (a,t)  =  — ^  f 

7T  JO 

Now,  let 

-(!)* 

where  t  =  \t\e^‘  and  0  <  <  7r,  to  change  (D.16)  into 

Few  {a, t)  =  ^/2  Jo  P(J  7  J  7  }  j/3/2 exp  |q  y  +  J  dy  (D.18) 

where 


0j4>t/  2 


ini 


(D.16) 


(D.17) 


(D.19) 


The  integration  contour  of  (D.18),  shown  in  Figure  D.2,  lies  in  the  valley  because  of 
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Figure  D.2:  The  integration  contour  in  the  complex  y  plane.  The  shaded  regions  are 
valleys  where  the  integral  in  (D.18)  is  convergent  and  the  saddle  points  at  y  =  ±1  are 
shown.  The  4>  is  this  figure  is  the  phase  angle  of  the  complex  time  variable  t  =  \t\ej<t>. 


the  following  observation 


6+6“6+2“2+6 


(D.20) 


for  0  <<f>t  <  7 r.  The  integration  contour  in  (D.18)  (shown  in  Figure  D.2  as  a  solid  line) 
is  shifted  to  a  new  contour  (shown  in  Figure  D.2  as  a  dashed  line)  which  goes  straight 
down  the  valley.  This  contour  shift  is  valid,  since  no  singularities  of  the  integrand 
exist  in  this  particular  shaded  valley  region  and  the  integrand  decays  exponentially 
in  the  valley  as  |t|  — >  oo.  So,  the  new  integral  expression  for  Few  (<*,£)  is 


f-  (CM)  =  “"'V’exp  [n  (-y  +  \y’) 


dy  (D.21) 


It  is  now  convenient  to  map  this  contour  integral  such  that  the  contour  is  along  a  real 
axis.  This  can  be  done  by  letting 


x  =  y  exp(— j7r/3  -  j<f>t/ 6) 


(D.22) 


and  noticing  that 


(l(-y  +  p3)  =  \n\e~j4>t/2  (-Xej*/3+j<t,t/6  - 

=  \n\  -  i*3) 

the  integral 

can  now  be  written  as 

Few  (a,t)  -  Tra5/2 

(D.23) 

where 

I(f)  =  H  x3l2e-*Wx*e~WAxdx 

Jo 

(D.24) 

and 

A  —  ei(7r_^‘)/3 

(D.25) 

Notice  that  Re  [A]  >  0  for  Im  t  >  0  (or  0  <  <f>t  <  7r).  Now  the  integral  I  ( t )  is  what 
must  be  evaluated  numerically.  First  I  (t)  is  broken  into  two  integrals 


where 


l,  (t)  =  [BXm  x3/2e-smx3e~M Axdx 

Jo 

I2  (t)  =  f  aj3/2e-®lnl*8e“'nl Axdx 
J  Bxm 


where 


(D.27) 

(D.28) 


(D.29) 


Notice  that  x  =  xm  is  the  location  of  the  maximum  of  x3l2e~*^xi  where  this  is 
considered  the  “slowly  varying”  factor  in  the  integrand,  whereas  e~^Ax  is  considered 
the  “rapidly  varying”  factor  in  the  integrand.  The  constant  B  is  chosen  such  that 
the  integrand  of  (D.24)  is  very  small  at  x  =  Bxm.  This  guarantees  that  Ii  (<)  is  the 
significant  term  in  (D.26)  and  I2  (t)  is  only  a  small  correction  term.  The  constant  B 


is  chosen  to  be 


B  -  2.7 


(D.30) 


A  different  value  for  B  could  also  be  chosen  as  long  as  B  >  1. 

The  I2  ( t )  integral  can  be  approximated  by  using  the  first  term  from  the  integra¬ 
tion  by  parts  procedure 


+  ( Bxm )=  exp  [-||fi|(i?a:TO)3  - 

'2  W  “  |f!|  [(£*„)*  +  4] 


(D.31) 


but  xm  =  (3/2)ly/3 |S7 1  1^3  so 


+  „  (3/2)^lexp[-lB3-(3/2)2|m^B] 

2  *  “  (3/2)>B2|fi|l  +  .4|n|5 


(D.32) 


1 

The  accuracy  of  the  above  equation  is  not  critical,  since  I2  ( t )  is  much  smaller  than 
Ii  (t).  Also,  notice  that  the  approximation  in  (D.32)  improves  when  |fl|  increases 
(i.e.  when  t  decreases). 

Now  all  that  remains  is  the  calculation  of  Ii  (<)  in  (D.27).  It  is  convenient  to 
make  another  change  of  variables  on  Ii  (t)  so  that  the  slowly  varying  factor  in  the 
integrand  is  independent  of  t  and  a.  It  is  very  useful  for  the  slowly  varying  factor  to 
be  independent  of  t  and  a  because  it  can  then  be  approximated  by  a  series  of  simple 
functions  just  once  and  the  integration  can  be  performed  in  closed  form  for  any  values 
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of  t  and  a  without  finding  a  new  approximation  for  the  slowly  varying  factor.  This 
turns  out  be  a  very  efficient  method  for  calculating  the  integral  Ii  (t).  So,  make  the 


following  change  of  variables 


to  arrive  at 


Ii  (f)  =  (Bxm)*  f  u*e  2B3u3e  ^ABxmUdu 

Jo 


(D.33) 


(D.34) 


Notice  that  the  slowly  varying  factor  in  the  integrand  u%e~^Bu7  is  now  independent 
of  t  and  a.  Also,  recall  that  the  constant  B  =  2.7  is  simply  a  fixed  constant.  Let  the 
slowly  varying  factor  in  the  integrand  be  called  g{u)  so  that 


t  \  ?  _tBsu3 

gyu)  =  u?e  2 


(D.35) 


Notice  from  Figure  D.3  that  g(u)  is  a  smooth  and  continuous  function.  The  next 
step  is  to  approximate  g(u)  with  a  series  of  simple  functions  which,  when  substituted 
into  (D.34),  allows  the  integration  to  be  done  in  closed  form.  One  could  achieve  this 
by  using  global  expansion  functions  defined  on  u  €  [0, 1]  or  by  using  subsectional 
expansion  functions.  The  algorithm  in  this  appendix  uses  global  expansion  functions. 
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Two  types  of  approximations  of  g{u)  are  now  explored,  one  uses  a  polynomial  and 
the  other  uses  a  series  of  exponential  functions.  The  polynomial  approximation  can 
be  easily  implemented  using  a  general  least  squares  procedure  [57].  The  approximate 
polynomial  can  be  written  as 

N 

g(u)  fa  Y  gnun  (D.36) 

n=0 

where  the  coefficients  {<7n}^o  are  f°un(i  from  the  general  least  squares  procedure. 
This  is  done  by  sampling  the  equation  in  (D.36)  at  many  values  of  u ,  obtaining  an 
overdetermined  system  of  linear  equations,  and  then  solving  this  system  of  equations 
in  a  least  squares  sense.  After  (D.36)  is  substituted  into  the  integral  of  (D.34),  the 
integration  may  be  done  term  by  term  to  obtain  the  following  approximation 

Ii  (f)  fa  ( Bxm)*  Y  9nK  (D.37) 

n= 0 


where 


ho  — 


(e-|0 \ABxm  _ 


-1  (p-\n\ABxm 

\Vt\ABx 

_  -\fl\ABXm 

hn  =  - + 


n 


\n\ABxm  \q\abxt 


Ln—  1 


for  n  =  1 ...  N 


(D.38) 

(D.39) 


Although  the  polynomial  approximation  for  g(u)  described  above  is  straight  for¬ 
ward,  an  exponential  series  gave  better  performance  (i.e.  the  exponential  series  ap¬ 
proximation  for  g{u)  is  more  accurate  for  a  smaller  number  of  terms).  An  exponential 
series  approximation  for  g(u)  may  be  written  as 


M 

g{u)  ta  Yj  9me~bmU  (D.40) 

m= 1 

where  the  parameters  {5m}m=0  and  {&m}„=o  can  obtained  by  using  the  extended 
Prony  method  (see  page  623  in  [61],  or  page  225  in  [62]).  Notice  that  the  {<7m} 
parameters  here  are  different  from  the  { gn }  described  earlier  for  the  polynomial  ap¬ 
proximation.  The  extended  Prony’s  method  provides  a  very  good  approximation 
although  it  is  not  guaranteed  to  give  a  “least  squares”  approximation.  For  the  prob¬ 
lem  of  concern  here  the  extended  Prony’s  method  gives  an  excellent  approximation. 
Tables  D.l  and  D.2  list  the  parameters  in  the  exponential  series  approximation  of 
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Table  D.l:  The  exponential  factors  for  the  series  of  exponentials  approximation  to 
the  g(u)  function  in  (D.40).  These  parameters  are  the  result  from  using  the  extended 
Prony’s  method  with  100  sample  points. 


m 

Re  (bm) 

Im  (6m) 

i 

2.492857990298293e+00 

-  1.838078138920067e+01 

2 

2.492857990298293e+00 

+  1.838078138920067e+01 

3 

3.699029046816519e+00 

-  1.178361372269887e+01 

4 

3.699029046816519e+00 

-f  1.178361372269887e+01 

5 

4.522095825495970e+00 

0 

6 

4.368725008737067e+00 

-  5.860005489119182e+00 

7 

4.368725008737067e+00 

+  5.860005489119182e+00 

8 

6.547833511389427e+01 

0 

9 

1.501055435645317e+02 

0 

10 

3.732424459349517e+02 

0 

(D.40)  when  M  =  10  and  (D.40)  is  sampled  100  times.  If  one  defines  the  relative 
error  of  the  approximation  as 

^n=l  [9\un)\ 

where  are  the  sampled  u  values,  then  ERR  =  2.937e  — 9  when  the  parameters 

in  Tables  D.l  and  D.2  are  used.  Finally,  substitute  the  approximation  in  (D.40)  into 
the  integral  (D.34)  to  obtain 


ia  ( t ) « {Bxmy 

,1/3 


£  (1 

2_s  9m  I  - 

m= 1  \ 


-  exp(— -  |fl| ABxn)' 
bm  +  \Sl\ABxm 


(D.42) 


where  xm  =  and  A  =  exj>\j(ir  -  <f>t)/3]. 

The  approximation  (D.42)  eventually  breaks  down  if  |t|  is  very  small  or  if  \t\  is 
very  large.  So,  when  |£|  <  0.021(a3/40)  use  the  first  three  terms  of  (D.12) 

3*5-7**- 15  /  t  N  21 


t  -  -3jT(5/2) 

Foi;  ~  va&/2 


,  nnnrrrfi\  3  -  5  -  7  -  *  •  15  /  t  V 
1  + 39.375  (-)+  3,27  (-) 


(D.43) 
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Table  D.2:  Coefficients  for  the  series  of  exponentials  approximation  to  the  g(u)  func¬ 
tion  in  (D.40).  These  parameters  are  the  result  from  using  the  extended  Prony’s 
method  with  100  sample  points. 


m 

Re  (pm) 

Im  (gm) 

i 

-1.309449928924531e-03 

-  4.143903979100756e-04 

2 

-1.309449928924775e-03 

+  4.143903979100027e-04 

3 

3.377694133016045e-02 

+  1.639435160611527e-02 

4 

3.377694133016031e-02 

-  1.639435160611565e-02 

5 

3.798712810824805e-01 

+  5.480557524830053e-16 

6 

-2.228939771564507e-01 

-  4.408366701346665e-02 

7 

-2.228939771564497e-01 

+  4.408366701346653e-02 

8 

6.785853728064312e-04 

-  1.098566438683437e-15 

9 

2.525319404875919e-04 

+  1.430649359954744e-15 

10 

5.057254010578827e-05 

-  5.842608514406963e-16 

When  |t|  >  15a3  use  the  first  three  terms  of  (D.14) 


+ 

Few 


(a,t) 


eJff/3 


r(5/6) 


r(7/6)e^/3a  r(9/6)e^/3a2 


2t2/3 


(D.44) 


Figures  D.4?  D.5  and  D.6  show  the  time  domain  creeping  wave  function  Few 
when  a  =  1  for  various  time  values.  The  function  is  computed  using  the  numerical 
algorithm  described  above  which  uses  a  series  of  exponential  functions  to  approximate 
the  “slowly  varying”  factor  in  the  integral.  Notice  that  when  Few  (<*>  t)  is  plotted  on 
a  large  time  scale,  as  shown  in  Figure  D.4,  it  appears  to  be  vary  sharp  and  almost 
singular  in  nature.  But  when  one  “zooms  in”  on  the  plot  of  Few  (<*,  t),  as  is  done  in 
Figure  D.6,  the  function  is  seen  to  be  smooth  and  continuous.  In  fact,  the  real  part 
of  Few  (a,  t)  is  infinitely  smooth  at  £  — >  0  since  all  of  its  time  derivatives  go  to  zero 
as  explained  earlier. 
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Figure  D.4:  The  time  domain  creeping  wave  function  Few  (ct,  t)  when  a  =  1,  plotted 
for  t  e  [—10, 10]. 


Figure  D.5:  The  time  domain  creeping  wave  function  Few  (<M)  when  a  =  1,  plotted 
for  t  £  [—1, 1]. 
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t 


Figure  D.6:  The  time  domain  creeping  wave  function  Few  (a,  t)  when  a  =  1,  plotted 
for  t  €  [—0.1,0.!]. 
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Appendix  E 

Computation  of  Special  Function 
Used  in  TD-UTD  Surface 
Diffraction 


+  p 

This  appendix  discusses  the  special  function  denoted  by  h  (E,t)  in  Chapter  4  used 
in  the  TD-UTD  surface  diffraction  coefficient.  Some  of  the  properties  of  this  function 

-j-  p 

are  discussed,  but  more  importantly,  a  method  to  compute  the  FSyh  (2,t)  for  complex 
time  values  t  and  both  positive  S  and  negative  S  values  is  presented. 

The  special  function  is  defined  as  the  analytic  time  transform  (ATT)  of  the  fre¬ 
quency  domain  Fock  type  function  (or  Pekeris’s  caret  function)  divided  by  a;1/6 


#*(S>0 


La; 


1  jW«,/3s)' 


1/6 


(E-1) 


where  E  is  a  parameter  which  depends  on  the  source  and  observer  locations  as  defined 
in  Chapter  4.  One  E  parameter  is  used  (instead  of  E  and  EL  as  is  in  Chapter  4).  If 
E  >  0  then  it  is  interpreted  as  the  shadow  region  parameter  in  (4.11),  but  if  E  <  0 
then  it  is  interpreted  as  the  lit  region  parameter  EL  given  in  (4.8).  This  should  not 
cause  any  confusion,  since  this  appendix  is  solely  concerned  with  the  computation  of 
the  function  in  (E.l)  with  little  regard  to  its  actual  physical  significance. 

-j -  p  * 

From  the  definition  in  (E.l)  it  can  be  shown  that  the  Fa<h  (S,  t)  function  satisfies 
the  following  scale  change  property  for  b  >  0 

(E.2) 
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K  (fS,()  =  pL  (H,i/63) 


(E.3) 


Notice  that  the  above  scale  change  property  is  the  same  as  for  the  time  domain 
creeping  wave  function  Few  («,£)  discussed  in  Appendix  D. 

E.l  Shadow  Region 

In  the  shadow  region,  the  E  parameter  is  positive  E  >  0.  When  the  observer  is  deep 
within  the  shadow  region  (or  when  u;E3  — >  oo),  the  frequency  domain  function  of 
(E.l)  may  be  written  as  a  series  of  creeping  wave  modes 


1  1  exp[-(ju’)1/3H5„ 

1  ^  1  exp[-(ju)''3S$„ 

VS^-IA  !(-«,)]’  y  «)>/• 


(E.4) 


This  creeping  wave  mode  series  of  (E.4)  is  asymptotically  valid  as  u;E3  — ►  oo  therefore 
one  can  use  (E.4)  to  find  an  expansion  of  the  analytic  time  transform  of  (D.l)  which 
is  asymptotically  valid  as  |t/E3|  — >  0.  The  corresponding  analytic  time  creeping  wave 


mode  series  is 


#*(2,0 


1  Few  (^ '9m  £) 

2[Ai,(-9n)]2 

1  &  Few  (5g„,t) 

2?„[Ai (~gn)]2 


(E.5) 


where  qn  is  the  zero  of  the  Airy  function,  Ai(— qn)  =  0,  and  qn  is  the  zero 
of  the  derivative  of  the  Airy  function  Ai7(— q„)  =  0.  See  Table  E.l  for  a  list  of  some 
of  these  roots  and  the  associated  values  of  the  Airy  function.  Using  the  table  and 
asymptotic  formulas  in  (E.6)  through  (E.13)  [29]  one  can  continue  the  list  in  Table 
E.l  ad  infinitum.  The  asymptotic  formulas  for  the  roots  of  the  Airy  function  are  [29] 


Ai'(-g„)  = 
Ai(-9„)  = 


/[3x(4n  -  l)/8] 
p[3x(4n  -  l)/3] 
(-l)"-7.[S*(4»-l)/3] 
(-l)”"'si[3ir(4n-  l)/3] 


(E.6) 

(E.7) 

(E.8) 

(E.9) 
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Table  E.l:  Zeros  of  the  Airy  function  and  associated  values  of  the  Airy  function. 


n 

Ai'(-9„) 

9n 

Ai(-g„) 

i 

2.33811 

+0.70121 

1.01879 

+0.53566 

2 

4.08795 

-0.80311 

3.24820 

-0.41902 

3 

5.52056 

+0.86520 

4.82010 

+0.38041 

4 

6.78671 

-0.91085 

6.16331 

-0.35791 

where 


f(  .  \  2/3  A-.  _i_  ®  2  5  4  77125  6 

/W  ~  *'  +  “35*  +  82944* 

108056875  __8  162375596875  __10\ 

6967296  2  +  334430208  Z  ) 


~  ,2/3 

V  48  288  207360 


18683371  „ 

H - z  - 

1244160 


91145884361 

191102976 


,  ,  ,  21/6  /  5  _o 

/l(2)  ~  (l  +  -Z~ 


1525  __4  2397875  __6 

4608 Z  +  663552  2 


*  (,)  ~Z-(l-  -  ***™m-A  (E.13) 

m  1  v^F  V  96  6144  26542080  J  K  ' 

For  observation  locations  near  the  shadow  boundary  (or  if  we  let  u;E3  0)  we 

can  use  a  power  series  representation.  This  is  obtained  by  first  using 


where 


e-W6  f  p*(wi/3H) 


(E.10) 


(E.ll) 

(E.12) 

(E.13) 


(E.14) 


(E.15) 


—  /  .  M/fi  N 

(yu;)1/6  |  J 

Now,  using  the  work  of  Logan  [63]  we  can  derive  the  following  power  series  for 

A.(S.")  =  E  A  {  P"  1 2"U«)(*‘l)  (E.16) 


r“  <r„ 
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Table  E.2:  Constants  used  in  the  power  series  expansion  of  the  Pekeris  functions  p(x) 
and  q(x)  and  therefore  also  in  the  expansion  of  FPig(E,u)  and  FP,q  (E,t) 


0  3.54064e-01  -3.07177e-01 


1 

-1.50139e-01 

2.63755e-01 

2 

-1.91020e-02 

-4.02720e-02 

3 

2.07797e-01 

-2.52283e-01 

4 

-3.04017e-01 

4.17454e-01 

5 

-1.68300e-02 

-3.34820e-02 

6 

1.16557e+00 

-1.37979e+00 

7 

-2.61483e+00 

3.13568e+00 

8 

-5.03520e-02 

-8.66800e-02 

9 

1.77043e+01 

-1.99933e+01 

13  -1.77677e+03  1.93449e+03 


14 

-3.27929e-01 

-4.55469e+00 

15 

2.45097e+04 

-2.62961e+04 

16  -9.71193e+04 


18  1.72330e+06 


19  -7.67284e+06 


20  -1.16244e+03 


1.03833e+05 


-6.78254e+01 


-1.82875e+06 


8.10881e+06 


-1.44881e+03 
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Table  E.2:  (continued) 


n 

Pn 

CTn 

21 

1.68425e+08 

-1.73093e+08 

22 

-8.27602e+08 

8.67552e+08 

23 

0.00000e+00 

0.00000e+00 

24 

2.18280e+10 

-2.23551e+10 

25 

-1.16876e+ll 

1.19558e+ll 

26 

0.00000e+00 

0.00000e+00 

27 

3.62274e+12 

-3.69876e+ll 

28 

-2.09315e+13 

2.13544e+13 

29 

0.00000e+00 

0.00000e+00 

30 

7.49321e+14 

-3.87390e+14 

31 

-4.63552e+15 

4.71980e+15 

32 

0.00000e+00 

0.00000e+00 

33 

1.88985e+17 

-1.69216e+17 

34 

-1.24387e+18 

1.26448e+18 

35 

0.00000e+00 

0.00000e+00 

36 

5.70892e+19 

-5.79813e+19 

37 

-3.97672e+20 

4.03714e+20 

38 

0.00000e+00 

0.00000e+00 

39 

2.03497e+22 

-2.06426e+22 
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Table  E.2:  (continued) 


n 

Pn 

40 

-1.49357e+23 

7.57602e+22 

41 

0.00000e+00 

0.00000e+00 

42 

8.45181e+24 

-8.65857e+24 

43 

-6.51126e+25 

6.66674e+25 

44 

0.00000e+00 

0.00000e+00 

45 

4.04588e+27 

-4.13808e+27 

46 

-3.26101e+28 

3.33365e+28 

47 

0.00000e+00 

0.00000e+00 

48 

2.21139e+30 

-2.25854e+30 

49 

-1.85943e+31 

1.89825e+31 

The  late  time  expansion  corresponding  to  the  power  series  expansion  in  (E.14)  and 
(E.16)  is 

+  p  +  e-W4 

F«h  (3,t)  =FPl9  (S,t)  -  27r5(_i<)1/2  (E.17) 

where 

E.2  Lit  Region 

In  the  lit  region,  the  parameter  E  is  negative  S  <  0.  In  Chapter  4,  the  lit  region 
parameter  is  denoted  by  EL  and  is  defined  in  (4.8).  Here  the  notation  S  is  used  to 
denote  both  S  >  0  for  the  shadow  region  and  E  =  EL  <  0  for  the  lit  region. 

For  an  observer  in  the  lit  region  near  the  shadow  boundary  (or  late  time  where 
|t/E3|  — >  oo),  one  can  use  the  power  series  in  (E.17)  and  (E.18). 

For  an  observer  deep  in  the  lit  region  (or  for  early  time)  where  |t/E3|  —*  0  one  can 
use  an  early  time  asymptotic  series  derived  from  a  high  frequency  asymptotic  series 


-jTr/12  oo 

f)  =  it(-jOS/S  £ 


Pn 


r(n/3  +  5/6)e 


jmr/6  77  ri 


n! 


(E.18) 
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in  Logan  [63].  Note  that  the  high  frequency  creeping  wave  mode  expansion  in  (E.4) 
is  not  valid  in  the  lit  region  where  E  <  0.  The  high  frequency  asymptotic  series  valid 
in  the  lit  region  E  <  0  is 


4A„( -1/3H)  -  ±*Z2e-*-W* 


g-jw(— E)3/12  N  Sn  1 

+  2(— 5)5/2  to)  K  I  [(-“)3iw]n+1 


(E.19) 


where  {sn}  =  {2,-20,560,...}  and  {hn}  =  {2,-28,896,...}.  Assuming  that  the 
ATT  can  be  applied  to  (E.19)  term-by-term  (see  Chapter  2)  then  the  early  time 
expansion  for  the  lit  region  (E  <  0)  is 


(E.20) 


for  |r/E3|  — ►  0,  where 


r  =  t-  (— E)3/12 


(E.21) 


Hn  (H,i) 


2 n\ 

.1  ft 


i  V  (— E)3 


In  -j 


(-E)= 


4>{n  + 1)}  (E.22) 


V’(l)  =  -7  . 

n-l  j 

ifj(n  +  1)  =  -7+^2  —  for  n  >  0 

_  -1  m 

771  —  1 


(E.23) 

(E.24) 


with  7  =  0.5772 . . .  being  Euler’s  constant.  Notice  that  the  expansion  in  (E.20)  is 
centered  at  r  =  0  or  t  =  (— E)3/12.  Because  of  this,  the  expansion  in  (E.20)  is  less 
useful  for  numerical  purposes,  because  it  is  not  centered  at  the  same  t  value  as  the  late 
time  expansion  in  (E.18).  So,  an  alternate  early  time  representation  may  be  needed. 

An  alternate  early  time  expansion  can  be  obtained  by  using  the  following  power 


series  form 


K  (2>o.  *  ±nr  f  [<  -  (-2)3/i2]  +  (ZipiE  Bnh  ((=55; 


(E.25) 
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From  the  results  in  Chapter  2,  it  is  known  that  (E.25)  has  the  proper  form  for  an 
early  time  representation  centered  at  t  =  0,  the  problem  lies  in  obtaining  the  values 
of  the  coefficients.  Notice  that  the  early  time  series  centered  at  t  =  0  in  (E.25)  is  an 
expansion  of  F  h  (E,  t)  minus  the  geometrical  optics  (GO)  term  containing  the  ana¬ 
lytic  delta  function.  This  formulation  is  very  effective  since  the  GO  term  is  dominant 
in  the  deep  lit  region  (where  |[t  —  (— E)3]/E3|  — >  0)  and  the  early  time  series  in  (E.25) 
therefore  represents  a  remainder  which  becomes  important  for  intermediate  time  val¬ 
ues.  Also,  notice  that  the  expansion  in  (E.25)  explicitly  satisfies  the  scale  properties 
in  (E.2)  or  (E.3).  It  is  known  by  numerical  experimentation  that  the  late  time  series 
in  (E.18)  can  be  trusted  (to  within  at  least  3  significant  digits)  for  |t/(E  )3|>  0.15. 
Because  of  this,  an  approximation  of  the  coefficients  in  the  early  time  series 

of  (E.25)  is  obtained  by  setting  the  early  time  expansion  in  (E.25)  equal  to  the  late 
time  expansion  in  (E.17)  with  (E.18)  at  |t/E3|  =  0.15  and  0  <  arg(t)  <  tt  in  a  least 
squares  sense.  This  provides  approximate  values  of  the  coefficients  {B„,h}  which  are 
most  accurate  at  |t/E3|  =  0.15.  This  is  very  desirable,  since  one  can  use  the  late 
time  series  representation  for  |t/E3|  >  0.15  and  then  switch  over  to  the  early  time 
representation  for  |£/E3|  <  0.15  without  any  significant  discontinuities.  Table  E.3 
lists  the  coefficients  {B„,h}  obtained  in  this  least  squares  sense.  Now,  since  the  co¬ 
efficients  {B°t'h}  are  obtained  in  this  approximate  manner  to  minimize  the  error  at 
|t/S3|  =  0.15,  the  early  time  series  in  (E.25)  may  not  work  as  |i/E3|  — »  0.  In  fact,  it 
is  known  from  numerical  experimentation  that  the  early  time  power  series  in  (E.25) 
may  not  be  accurate  when  1 £/E  3|  <  0.15  and  Im  (£/(— E)3)  <  0.002.  Nonetheless, 
when  the  imaginary  part  of  t  is  this  small  (  Im  (</(— E)3)  <  0.002)  the  analytic  delta 
function  (i.e.  the  GO  term)  in  (E.25)  is  the  dominant  term  so  for  practical  applica¬ 
tions  this  representation  works  quite  well.  On  the  other  hand,  if  one  is  interested  in 
computing  an  impulse  response,  where  Im  (£)  =  0  it  is  difficult  to  compute  the  early 
time  behavior  besides  the  dominant  impulse  time  function. 
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Table  E.3:  Coefficients  {£*}  used  in  the  early  time  representation  in  (E.25)  used  in 
the  lit  region  E  <  0  and  for  soft  polarization 


n 

Re  [B*] 

Im  [BJ] 

i 

-2.082816029604068e+09 

-  1.934225071410909e+09 

2 

2.231501364264926e+09 

+  2.248973125688412e-f09 

3 

-1.252865898123653e+09 

-  1.373344664977570e+09 

4 

4.788187068366086e+08 

+  5.744350921 182956e+08 

5 

-1.371893503363450e+08 

-  1.820181556200502e+08 

6 

3.078719411539325e+07 

+4.590574642077959e+07 

7 

-5.503967634776822e+06 

-  9.456164234707374e+06 

8 

7.803327956948597e+05 

+  1.609576880445286e+06 

9 

-8.475429971667759e+04 

-  2.264583143400559e+05 

10 

6.299138967200948e+03 

+  2.60086556591 7071e+04 

11 

-1.674918333321921e+02 

-  2.357302083687360e+03 

12 

-2.947872815842084e+01 

+  1.543029977826218e+02 

13 

4.483644084797442e+00 

-  4.980881275349820e+00 

14 

-6.216333948603559e-02 

-  4.140469301332007e-01 
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Table  E.4:  Coefficients  {J9„}  used  in  the  early  time  representation  in  (E.25)  used  in 
the  lit  region  E  <  0  and  for  hard  polarization 


n 

Re  K1 

Im  K] 

i 

-1.145347064563553e+09 

+  2.321104859253427e+08 

2 

1.247217383801985e+09 

-  6.159421376984853e+07 

3 

-7.075526321707213e+08 

-  7.680045275028466e+07 

4 

2.716051094355760e+08 

+  7.653405928940107e+07 

5 

-7.758226975004484e+07 

-  3.753779340973016e+07 

6 

1.717005587168613e+07 

+  1.266153768031792e+07 

7 

-2.972294299624093e+06 

-  3.237541175734916e+06 

8 

3.933478822263941e+05 

+  6.531900019204784e+05 

9 

-3.619503051989348e+04 

-  1.054006228018756e+05 

10 

1.374538560322166e+03 

+  1.350912057659442e+04 

11 

2.200490845152954e+02 

-  1.327895365123774e+03 

12 

-5.070700867343560e+01 

+  9.003798728735892e+01 

13 

5.004751143977662e+00 

-  2.421287171291755e+00 

14 

-7.860467067425439e-02 

-  3.708472052078492e-01 
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E.3  Numerical  Algorithm 


This  section  presents  an  algorithm  using  the  results  of  this  appendix  for  the  compu¬ 
tation  of  the  Fsth  (S,£)  for  any  real  value  of  E  and  for  any  complex  value  of  t  for 
Im  (t)  >  0. 

Before  the  numerical  algorithm  for  the  computation  of  the  Fa<h  (2,f)  is  presented, 
a  small  correction  term  for  the  creeping  wave  series  of  (E.5)  is  derived.  In  general, 
one  would  like  to  use  the  early  time  representation  in  (E.5)  up  to  some  time,  say 
|t/S3|  =  To,  and  then  switch  over  to  the  late  time  representation  in  (E.17)  (with 
(E.18))  for  |f/E3|  >  T0.  But,  in  order  to  do  this  the  two  representations  (late  time 
and  early  time)  must  overlap  at  |t/E3|  =  To  so  that  the  resulting  function  does 
not  contain  a  discontinuity.  It  is  known  from  numerical  experimentation  that  the 
creeping  wave  mode  series  of  (E.5)  for  soft  polarization  is  not  accurate  enough  for 
large  |t/S3|  values  even  when  50  creeping  wave  modes  are  included  in  the  summation. 
An  estimate  for  the  infinite  number  of  remaining  creeping  wave  modes  (from  the  51st 
on  to  infinity)  can  be  obtained  as  follows.  First,  write  the  following  creeping  wave 
mode  summation  for  soft  polarization  as 


O  _  V->  Few  ('— Q 

"0_J^02[Ai  \-qn))2 


(E.26) 


where  Nq  =  51  for  the  numerical  algorithm  to  be  described  shortly.  Now,  from  the 
asymptotic  formula  in  (E.6)  with  (E.10), 


37r(4n  —  1) 

8 


(E.27) 


and  for  n  >  51  this  is  accurate  to  five  significant  figures.  From  the  asymptotic  formula 
in  (E.8)  with  (E.12), 


Ai'(-9„) 


(— 1)"  1  f  37r(4n  —  1) 


(E.28) 


which  is  accurate  to  six  significant  figures  when  n  >  51.  When  |f/s3|  <  1  and  n  >  51 
the  creeping  wave  mode  function  Few  (<*,  t)  may  be  approximated  with  the  first  term 
in  the  power  series  expansion  of  (D.12)  as 


+,  ,  l  ,,  — 3jT(5/2) 

Few  (^5 1)  ~  5  /2 

7T  na/* 


(E.29) 
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which  is  accurate  to  within  0.1  percent.  Equations  (E.27),  (E.28)  and  (E.29)  can  be 
substituted  into  (E.26)  to  obtain 

— 3jT(5/2)  /  8  \2  1 

""  ~  2H»/2  U  J  (in-  l)2 

which  is  approximate  for  N0  =  51  and  asymptotically  valid  as  No  —■ ►  oo.  The  Euler- 
Maclaurin  summation  formula  can  be  used  to  approximate  the  infinite  summation  in 
(E.30)  [57]. 


<30  J 

So  (4n  -  !)2 


4(41Vo  -  1)  +  2(41Vo  -  l)2  +  12(4JV0  -  l)3 
1536 

_720(41Vo  -  l)5 


(E.31) 


which  gives 

o°  I 

£  rr  _  V,2  «  1.243740 JE?  -  3  (E.32) 

n=51  I4”  11 

Finally,  the  infinite  summation  of  remaining  time  domain  creeping  wave  mode  func¬ 
tions  is 


Sn0  ~ 


-jl.786878E  -  3 

-5/2 


for  No  =  51 


(E.33) 


*4“  p 

Now,  the  numerical  algorithm  which  can  be  used  to  compute  the  Fah  (H,t)  func¬ 
tion  is 


•  For  observers  in  the  shadow  region,  S  >  0: 

—  When  |t/H3|  <  1  use  the  creeping  wave  mode  summation  in  (E.5),  with 
Nh  =  20  and  Ns  =  50.  Also,  add  the  extra  term  in  (E.33)  times  (  —  1/sqrtn ) 
to  the  creeping  wave  mode  sum  in  (E.5)  for  the  soft  polarization  case, 
approximating  an  infinite  creeping  wave  mode  sum. 

—  When  |t/H3|  >  1  use  the  late  time  inverse  power  series  in  (E.17)  with 
(E.18)  including  50  terms  in  the  series. 

•  For  observers  in  the  lit  region,  E  <  0  : 

-  When  |t/E3|  <  0.15  use  the  approximate  early  time  representation  in 
(E.25)  including  14  terms  in  the  power  series. 
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-  When  | t/E  3|>  0.15  use  the  late  time  inverse  power  series  in  (E.17)  with 
(E.18)  including  50  terms  in  the  series. 

The  algorithm  described  above  is  not  the  most  efficient  method  of  using  the  equa¬ 
tions  of  this  Appendix,  but  it  is  robust  and  works  well  as  long  as  Im  ( t )  >  0.002 |H3 1. 
Future  research  can  be  done  on  more  efficient  representations,  such  as  rational  func¬ 
tion  approximations  and  so  forth.  The  challenge  of  approximating  the  analytic  time 
function  Fa<h  (2,£)  hes  in  the  ability  to  compute  this  function  for  all  complex  time 
values  in  the  upper  half  time  plane  ( Im  (£)  >  0).  This  is  a  more  difficult  task  then 
the  more  common  problem  of  approximating  a  given  function  of  a  real  variable. 

Figures  E.l  and  E.2  show  a  comparison  between  the  late  time  representation 
for  F^h  (2,i)  found  in  (E.17)  and  (E.18)  versus  the  early  time  creeping  wave  mode 
representation  for  F^h  (H,£)  found  in  (E.5)  when  the  observer  is  in  the  shadow  region 
(E  >  0).  The  geometry  parameter  E  is  set  equal  to  E  =  1  for  Figures  E.l  and  E.2 
without  loss  of  generality,  since  a  different  S  value  only  causes  a  simple  scale  change  as 
in  (E.3)  but  leaves  the  comparison  valid.  Notice  the  the  two  representations  overlap 
very  closely,  except  when  |t|  is  very  small  such  as  |</S3j  <  0.1.  Recall,  in  the  numerical 
algorithm  outlined  earlier  that,  in  the  shadow  region,  the  early  time  creeping  wave 
mode  representation  of  (E.5)  is  used  for  |f/E3|  <  1.0  while  the  late  time  inverse  power 
series  of  (E.17)  with  (E.18)  is  used  when  |f/S3|  >  1.0. 

Figures  E.3  and  E.4  show  a  comparison  between  the  late  time  representation 
for  Fsh  (2,<)  found  in  (E.17)  and  (E.18)  versus  the  approximate  early  time  power 
series  representation  in  (E.25)  when  the  observer  is  in  the  lit  region  (i  <  0)  and 
the  imaginary  part  of  time  is  held  constant  at  Im  (£/(— S)3)  =  0.01.  The  geometry 
parameter  E  is  set  equal  to  E  =  —1  for  these  figures  without  loss  of  generality,  since 
different  S  values  simply  cause  a  change  in  scale.  Notice  that  the  two  representations 
overlap  very  closely  when  Re  (£)/(— S)3  ~  0.15  and  more  generally  they  overlap  very 
closely  whenever  |£/E3|  =  0.15.  Also,  the  early  time  representation  breaks  down 
violently  as  |£/E3|  — >•  0  and  the  approximate  early  time  representation  breaks  down 
as  |£/E3|  — >  oo. 
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Figure  E.l:  The  special  function  F/,  (H,t)  when 
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